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Chapter 1 

Introduction 

1.1 Overview of study 

My work has dealt with a sequence of basic plasma physics questions that range over 
a wide territory of subject matter. The common aim of these questions is to better understand 
the nature of high speed flow of a magnetized plasma. One important aspect of this work is the 
physics of turbulent flow states. The device that makes this investigation feasible is the CTIX 
plasma accelerator. However, our great success at designing and operating this accelerator has been 
founded on only a very rudimentary model of the plasma physics involved. Like early hominids, 
not understanding combustion chemistry hasn't stopped us from rubbing two sticks together and 
making fire. 

As I have endeavored to do something new and useful with CTIX, I have been repeatedly 
confronted with fundamental questions regarding the plasma physics that makes it all work in the 
first place. I have tried to answer these questions as they have become experimentally accessible, 
and I have attempted to synthesize our results into improved models of what is actually going on 
inside the plasma. The individual results I have produced are interesting in themselves, in that 
they advance our understanding of the inner workings of CTIX-like plasma accelerators; but as a 
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collection I hope that they serve a more general purpose by providing a better understanding of the 
range of possible initial conditions and plasma parameters that give rise to fundamentally interesting 
phenomena such as magnetohydrodynamic (MHD) turbulence. In addition, I have tried to develop 
a useful framework of well posed questions that will serve to guide future progress. 

I have organized my results in a sequence from chapter 2 to chapter 7 in a way that 
logically establishes the key prerequisites for an MHD turbulence experiment. Chapter 2 examines 
the structure and dynamics of the compact toroid's magnetic field. Next, in chapter 3 I develop 
some important results about the electrical resistivity of the plasma. This information is critical 
for properly analyzing the Reconnection/Compression experiments described in chapter 4. Together 
this work provides the initial conditions for the drift section measurement of kink dynamics and the 
wire target turbulence experiments discussed in chapters ?? and 9. The resistivity measurements 
are of key importance in the estimation of the magnetic Reynolds number of the plasma flow, which 
we will use to characterize the flow conditions in our turbulence experiments. 

Chapter 5 describes measurement of the electron temperature of the plasma and looks 
at the dynamics of energy conversion during the process of interaction with a transverse magnetic 
field. The knowledge about the electron temperature helps to construct an accurate model of the 
ionization of helium under our range of conditions. Helium enters the equation because of its ability 
to emit light brightly when bombarded with plasma electrons. This optical emission provides several 
diagnostic opportunities. The first is explored in chapter 6. Using a high resolution spectrometer 
we were able to measure the Doppler shift of a spectral line of He II ions. In doing so we measured 
the velocity of the plasma in previously inaccessible locations in the vessel. 

Several of these preliminary results are combined in Chapter 7 where I calculate values of 
the Reynolds numbers that occur in the CTIX plasma flow. The dimensionless Reynolds numbers 
play a vital role in determining the properties of a fluid flow, and the issues regarding turbulence 
that we can address with CTIX are focused on how experimental observables change with changing 
Reynolds number. 
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From Chapter 8 on we direct our attention to the experiments in the drift section of CTIX 
where conditions are suitable for exploring MHD turbulence. Chapter 8 discusses our primary tool 
for investigating turbulence, a technique called Gas Puff Imaging (GPI) in which a cloud of neutral 
helium can be used to image plasma density fluctuations. Chapter 9 introduces a wire target that 
was used in a set of experiments to act as a magnetic perturbation to the plasma flow, and presents 
the analysis of images that show coherent and incoherent density waves that are created in the 
presence of the wire target perturbation. 

Proposed future work and and a variety of technical details are included as appendices. 

1.2 Structure of CTIX device 

The Compact Toroid Injection experiment, (CTIX), is a coaxial railgun that forms and 
accelerates magnetized plasma rings called compact toroids (CT's). CTIX consists of a pair of cylin- 
drical coaxial electrodes with the region between them kept at high vacuum. The outer conductor 
vessel is 2.5 m in length and 15.24 cm in diameter. The electrodes are connected to high volt- 
age 50 fiF capacitors and metglass passive timing delay inductors that together form the rail gun 
circuit. The electrically conductive CT plasma forms a sliding short between the inner and outer 
electrodes, and completes the path for current to flow down the center conductor, across the plasma, 
and back along the outer conductor. The railgun effect that accelerates the CT can be accounted 
for by the Lorentz j x B force density, where j is the current density driven by the external circuit, 
and B is the magnetic field created by the current flowing in circuit, (primarily down the center 
conductor). There are separate capacitor banks for the formation and accelerator sections that can 
be charged to different voltages. The final velocity of the CT can range from Vct = hem j [is up 
to Vct = 25cm/ /is. Vct can be adjusted within this range by changing the applied accelerator 
voltage. The mass of the injected hydrogen CT is typically in the range of a few micrograms. 

The center electrode has a larger diameter in the formation region, has a step down in 
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Figure 1.1: 3D view of CTIX showing magnetic probe positions and circuit schematic 

diameter, and then a tapered section, and finally it has a long, slightly thinner straight section 
where most of the CT acceleration occurs. The outer electrode has a corresponding tapered section 
at the junction between the formation and acceleration regions. Some of these details are not 
included in the graphical rendering of CTIX displayed in figure 1.1. 

Originally, the center electrode was electroplated with gold to minimize surface effects 
from the raw stainless being exposed to the plasma. However, after years of operation, the sustained 
direct impact of massive high current discharges at gigawatt power levels has resulted in most of 
the gold plate on the middle 10 cm of the formation section center conductor becoming completely 
removed via ablative processes and deposited elsewhere in the vessel. Droplets of gold, typically 
having diameters less than 2 mm, are occasionally found toward the end of the accelerator, and in 
the drift and target chambers. These have apparently been accelerated during a plasma discharge 
and seem to have cooled in mid flight because of their resulting tear-drop shape. It is unknown how 
these droplets affect the electrodynamics of the acceleration discharge. 

Also, some regions near the step down edge on the center conductor (z — 20 cm) have 
accumulated the gold deposits as a highly textured surface that looks like cracked stucco or peeling 
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paint. We don't know exactly what effect this surface has on the formation process, but all evidence 
suggests that either it does not matter at all, or it actually improves the formation of high density 
compact toroids. 

The possibility that this surface texture could be beneficial is interesting, and not incon- 
ceivable. Electric fields will be significantly stronger at the tips of the gold flakes than elsewhere 
in the bulk volume. Positive contributions to the formation process might be accounted for simply 
by an enhanced secondary electron generation at the surface, or some kind of hollow-cathode effect, 
in which the small flat spaces between the flaking gold layer and the vessel wall act to create very 
efficient fast ionization of the plasma in a large number of little pockets distributed over a large 
surface area. Additionally, such an unusual formation boundary condition might partially account 
for the fact that the neutral gas break-down event is held off for such a long time (several ms) while 
the seed magnetic field decays until it is at just the right magnitude and field geometry, and then 
break-down spontaneously occurs. It is also possible that such a rough surface on the boundary of 
the newly forming plasma could modify the fluid velocity field such as by exciting high k turbulence 
for example, which ultimately results in the mysteriously efficient fast dynamo effect that is respon- 
sible for the creation of the large internal magnetic field of every single compact toroid we generate 
on CTIX. 

The observations of extended hold-off of break-down, and unusually high magnetic flux 
amplification factor (~ 200) that occur during the plasma formation are two of the most basic 
features on the CTIX device. They are definitely in need of further explanation, and the effects of 
a non-trivial boundary condition should be an important element to consider in future work on the 
subject. 

CTIX was situated to be able to inject its CT's into the Davis Diverted Torus (DDT), 
which was a small toroidal device that was built in order to investigate basic tokamak physics. It 
was composed of a stainless steel toroidal vacuum vessel, major radius of 44 cm, minor radius 14 cm, 
which is surrounded by 7 toroidal field coils, as well as a variety of poloidal field coils, an Ohmic core, 
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and a set of pre-ionization coils. DDT originally operated in a pulsed mode where a tokamak-like 
plasma was heated and confined. The DDT toroidal field coils had a radius of 40 cm and produced 
about 3 Gauss per Amp at the center of the minor axis of the vacuum vessel. They were water-cooled 
and have achieved a maximum field of 600 Gauss in steady-state mode. 

DavisU Diver tor Torus 



CTIX Injector 
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Figure 1.2: Top view of the CTIX-DDT facility 



However, as the focus of experiments shifted toward CTIX we ceased the full-power pulsed 
mode and only operated DDT with smaller currents in a steady-state mode. In this mode a tokamak- 
like plasma is not possible or desired. Instead, we injected the CT into the vacuum field of DDT 
so that we could observe the interaction of the CT and the external toroidal field. The presence 
of tokamak plasma would have made the diagnostic situation more difficult in certain ways, mostly 
because once reconnection has occurred there is no easy way to distinguish between the tokamak 
plasma and the CT. Without the tokamak plasma getting in the way of our diagnostics, we can look 
at zeroth order stopping effects, as well as post-reconnection dispersion and thermalization since the 
only plasma in the device is that which originated in the CT. For a more detailed description of the 
phenomenon of magnetic reconnection see chapter 4. [ref NF high beta paper] 

Future work at a new experimental facility* will be conducted with CTIX injecting directly 
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into a smaller diagnostic chamber, with DDT removed from the system. This will significantly 
decrease the total volume of the vacuum vessel and allow higher repetition rate because of the 
shorter pump-out time. Also the new design will greatly improve diagnostic access and usability of 
the system. 

1.3 Concept of a compact toroid plasma 

A compact toroid is a magnetized plasma ring formed within a conducting vessel of coaxial 
geometry. The term "compact" is used because traditionally these plasmas tend to be formed 
between coaxial electrodes having only a modest outer radius of about 10 cm. 

It has a strong interior magnetic field with a special geometry that has closed flux surfaces, 
resulting in plasma particles being trapped on the field lines for much longer than the escape time 
due to the particle's thermal velocity. Also, this magnetic field is self-generated by the currents 





Toroidal 



Figure 1.3: Visualization of magnetic flux surfaces within a compact toroid 
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flowing within the plasma ring. In this way the compact toroid is self-confined in a stable geometry 
that persists until the internal currents resistively die down or particles escape as certain non-trivial 
phenomena eventually dominate. 

This special magnetic field geometry has comparable poloidal and toroidal components at 
each point within the plasma. Field lines wrap around toroidal surfaces called flux surfaces, in 
this way field lines can be grouped by what flux surface they lie on, and the entire magnetic field 
constitutes a set of nested tori of flux surfaces. 

The fields are maintained by associated toroidal and poloidal plasma currents. Its magnetic 
geometry is much like a spheromak, the distinction between them is primarily that a compact torus 
is being accelerated by a large railgun current, whereas a spheromak is formed at rest with respect 
to the laboratory. Also, spheromaks are formed in large, roughly spherical vessels, while the coaxial 
vessel allows the compact toroids to have almost any aspect ratio, and can range from doughnut 
shapes to long, cigar shapes. 

Measured surface magnetic fields at the edge of the CT are large B ~ 1 kGauss to 5 kGauss 
and internal magnetic fields are expected to be a factor of 3 to 10 times larger. 

A more detailed account of the theory of the magnetic equilibrium of compact toroids and 
spheromaks will be presented in Chapter 2. The CT has a long lifetime of over 1 ms and is able to 
hold itself together under some rather large forces. The compact toroids formed on CTIX can be 
characterzed according to two parameters, the thermal beta, and the kinetic beta. 



The thermal beta is the ratio of the hydrodynamic thermal pressure to the magnetic pressure, while 
the kinetic beta is the ratio of kinetic energy density to magnetic energy density. Part of what makes 
an accelerated compact toroid plasma special is that it has a low thermal beta, yet a high kinetic 
beta, Pthermai -C 1, yet Akinetic ~ 1- This is important because it means that magnetic effects will 
dominate many pure hydrodynamic effects, yet there is an abundant source of kinetic energy that 
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can drive plasma phenomenon into the non-linear range. 

There is a large toroidal magnetic field behind the CT, initially generated by the railgun 
current flowing in the center conductor. This field becomes spatially uniform and steadily decays as 
the column of magnetized plasma behind the CT expands. During this expansion magnetic energy 
is converted into the kinetic energy of the plasma, thereby driving the acceleration of the compact 
toroid, which rides at the leading edge of the expanding column. 

As it leaves the accelerator, CT plasma is relatively dense with n = 10 14 particles per cm 3 
and of modest temperature T e = T, = btieV ~ 580,000°K. 

The total visible emission of the plasma has been measured to have a peak value of about 
10 17 photons/(cm 3 • sec). This is enough for many amplified or high-sensitivity optical detectors 
to get a good signal, if they have long enough exposure (~ 1 — 10 fis), and integrate over a wide 
enough wavelength range. However, for certain challenging measurements such as high-speed, high- 
resolution spectroscopy, and fast imaging through optical filters, special techniques are required to 
increase either plasma brightness or detector sensitivity. The CT plasma and is optically thin in the 
visible range (i.e., it is transparent to its own light). 

1.4 Timescales within the life of a CTIX plasma 

The innovative aspect of the CTIX design is that it is capable of sustained operation with 
a maximum firing rate of 0.2 Hz, (1 shot every 5 seconds). We will often refer to the firing rate of 
CTIX as its rep rate or repetition rate. When we were taking data the rep rate was typically 0.05 
Hz (1 shot every 20 seconds) so that the vacuum vessel could be fully pumped back down to base 
pressure between shots. 

The standard firing rate is sufficiently fast to be able to take 1000 shots a day. If the 
pumping rate was significantly increased (or injected gas volume was decreased), in principle the 
repetition rate could be increased to 1 kHz if the capacitors could be maintained at high voltage. 
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Figure 1.4: Timescales during a CTIX acceleration discharge 



CTIX's rapid fire ability is due to the fact that it is passively switched. The only thing 
we actively control on CTIX is the gas valve trigger pulse that causes a puff of gas to be let into 
the formation section, already energized to a steady high voltage of about 10 kV. Then at some 
randomly occurring time later, the gas ionizes due to the formation potential, thereby completing 
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the formation circuit. 

After a fixed delay of about 4 /j,s that allows the magnetized plasma to expand into the 
acceleration region and relax into a CT geometry, the accelerator circuit then fires, and accelerates 
the CT to its injection velocity. 

The whole system is passively dependent on the event of the gas break-down, which means 
that CTIX can fire any number of times in quick succession at whatever rate you can successfully 
send in distinct pulses of gas. One limitation is that presently the capacitors loose their charge 
after only one firing cycle. Recharging the capacitor bank takes about 3 seconds when the system 
is running in a fully automated mode. So, in principle we could fire once every 3 second with the 
present machine. However, in practice we run at a much slower firing rate to accommodate several 
other factors. First, the data acquisition system (CAMAC + LabView) typically requires about 
5 seconds to record all the data. And second, as previously mentioned, the pumping rate of our 
present vacuum system is only fast enough to allow a repetition rate of one shot every 20 seconds. 
Any faster than this and the gas pressure will steadily rise over the course of a run. If the electrodes 
a re-energized when the base pressure is too high (above 10 -4 Torr), premature arcing will occur 
that can damage the device. 

The basic operations of CT formation and acceleration have been demonstrated to operate 
very robustly, with a nearly failure-free production of compact toroids over the course of more than 
60,000 firings of CTIX. Although CTIX almost never fails to produce and accelerate viable CT's, 
there is a noticeable amount of random variation of CT properties from one shot to another. With all 
experimentally controllable settings of the device held fixed, important observables such as plasma 
density, CT magnetic field shape and strength, and CT final velocity, will jump around randomly 
within some finite range. These fluctuations are attributed to uncontrollable variations in the initial 
conditions during the CT formation process. 

For example, the formation gas valve sometimes lets in a little more hydrogen then aver- 
age, sometimes less, resulting in a correspondingly greater or lesser CT plasma density. Another 
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uncontrollable event is the ionization of the neutral hydrogen (called the breakdown), which will 
occur at a randomly occurring time, that can vary as much a several milliseconds. 

The time of breakdown is likely to depend on a variety of system parameters in a complex 
way. It is expected that the peak density of neutral hydrogen should play an important role, 
the strength and geometry of the formation magnetic field that acts to bend the trajectories of free 
charged particles, and it most critically depends on the spontaneous creation of enough free electrons 
(possibly via cosmic rays) to trigger the needed runaway cascade of electron impact ionization. 
Overall, there are many more questions than answers regarding the details of the formation process. 



1.5 Typical CTIX plasma parameters 



Typical CTIX Parameters 


Quantity 


Units 


Low 


High 


Density n e = rii 


[cto~ 3 ] 


10 i4 


10 ib 


Electron Temperature 


[eV] 


10 


80 


Peak magnetic field 


[kGauss] 


1 


10 


Electron collision rate 


[MHz] 


4 


900 


Debye Length \n 


[cm] 


7 x 10~ b 


7 x 10~ 4 


Electron gyroradius 


[cm] 


7 x 10~ 4 


0.2 


Fluid Reynolds number VctL/v 


100 


2 x 10 4 


Magnetic Reynolds number VctL/Dm 


10 


10 4 



Table 1.1: Typical range of CTIX plasma paramters 



1.6 Overview of diagnostics on CTIX 



1.6.1 Probe diagnostics 

Along CTIX there are a number of magnetic field probes. Each is capable of resolving 
poloidal and toroidal components of the magnetic field at the outer edges of the CT. These different 
components are recorded on separate digitizer channels. We have three magnetic probes situated 
along the length of the accelerator at 57 cm, 91 cm, and 142 cm from the gas valve. The B z probe 
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signals are used on a daily basis for to see the existence of the CT and to get a quick estimate of 
its average velocity and acceleration using time-of-flight kinematics. Details of the magnetic probes 
will be discussed next, in chapter 2. We also use magnetic probes throughout the drift section and 
DDT. 

A simple measurement of plasma density can be made with Langmuir probes, consisting 
of a thin coaxial line that has its center conductor protruding at the probe tip, allowing the metal 
to be in direct contact with the plasma. By charging the probe tip to a bias voltage relative to 
the plasma, either electron current or ion current will be collected by the wire. Some temperature 
information can be acquired by observing how the current varies when the bias voltage is scanned 
over a range. 

Operated in a similar manner, we have another kind of plasma probe called an energy 
analyzer probe (EAP). It uses a set of fine metal grids to repel charged particles, but allowing a 
high-energy subset of the population to be collected by a plate in the back of the probe assembly It 
is designed to measure the energy distribution function of the plasma, and therefore the temperature, 
by way of scanning the bias voltages applied to the repeller grids. More details of the energy analyzer 
probe will be discussed in chapter 5 

1.6.2 Optical diagnostics 

One of our most complex diagnostics is our set of three laser interferometers that measure 
the electron density of the plasma by detecting changes in the index of refraction. This is done using 
a long coherence-length He Ne laser, whose beam is split into two sub-beams, one that passes through 
the plasma, and a second reference path of the same length as the first that only passes through air. 
These two beams are mixed back together, and the result is that phase of the combined light will 
undergo constructive and destructive interference by an amount that depends on the difference in 
refractive index along the two paths. 

A heterodyne mixing technique is used to efficiently extract the phase information from 
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Phase shifted 
Signal beam 



Frequency shifted 
Reference beam 



Bragg Cell AOM 
40 MHz Modulation 




Plasma in CTIX 



HeNe 5raW laser 



Fiber out to heterodyne 
electronics 



Figure 1.5: Schematic of He Ne laser heterodyne interferometer 



fast time-domain measurements of the laser intensity. A key feature in our system is that the 
plasma passing through the interferometer laser chord results in a very rapid change in phase, 
which can easily be distinguished from much slower modulations do to mechanical oscillations of the 
interferometer support structure (shown in fig. 1.5). We have three such interferometers, positioned 
at the optical ports at z = 57, 91, 142 cm. 

A very attractive alternative to the interferometer is the laser deflectometer, which also 
accomplishes electron density measurements, but it does so by detecting microscopic deflections of 
a laser beam passing through the density gradients that exist within our finite sized plasma. It uses 
a low cost diode laser, and a four quadrant laser spot detector in combination with a custom-made 
amplifier circuit to measure the laser deflection as a function of time. Special care needs to be 
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given when designing the frequency response of the circuit to make in insensitive to low frequency 
mechanical vibrations, as well as very high frequency electrical noise. In addition to low cost, it has 
the advantage that it is relatively easy to align, and is well suited to detecting small-scale structures 
that have steep density gradients, such as would occur in shock waves, and possibly in turbulence. 
The dcflcctomctcr is currently being developed by Samuel Brockington as part of his Ph.D. project. 

The plasma actively emits visible light, and there is a lot that we can do with this fact. 
One route is to measure the time dependence of the emission within some range of wavelengths, 
and localized along a narrow chord in the plasma using optically filtered photo-multiplier tubes 
(PMT). The plasma light is collimated using a pair of small adjustable apertures mounted at the 
ends of a long tube. Before the light can be detected by the PMT, it must first pass through an 
optical bandpass filter, allowing only a small range of wavelengths through. The PMT is needed 
because its high sensitivity and dynamic range allows detection of the very small fraction of plasma 
light that makes it through the apertures and narrow filters. This has been very useful in quantifying 
individual atomic processes within the plasma. The use of the filtered PMT assemblies to measure 
the excitation of singly ionized helium will be discussed in chapter 6. 

We can also take digital photographs of the plasma using a variety of high speed imaging 
devices, including a high-sensitivity CCD, a micro-channel plate intensified CCD, and a fast framing 
Imacon tube. These can also be used to view the plasma light through optical filters in order to 
distinguish different atomic processes. Fast imaging is primarily useful in reveling the existence of 
global structures within the plasma that would have gone unnoticed using only point measurements. 

Lastly, we can use spectrometers to quantify how plasma light is distributed in wave- 
length space. Typically the visible plasma spectrum is dominated by line emission from atomic and 
molecular transitions. We have used two different spectrometers on CTIX, that differed in spectral 
resolution and sensitivity. The details of these will be discussed in chapter 6. 
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1.6.3 System monitors 

A number of current and voltage monitors are used to diagnose what is going on within the 
CTIX accelerator circuit. These included the current applied to operate the formation gas valve, 
and formation solenoid. We routinely measure the voltages across the formation and accelerator 
gaps, as well as the currents flowing in the formation and acceleration electrodes. Figure 1.4 shows 
the waveforms generated by some of the primary system diagnostics. 

The vacuum pressure within the vessel is continuously monitored using low pressure ion 
gauges that can operate below 10 -4 Torr, as well as Convectron gauges the are sensitive above 1 
milliTorr. The higher pressure gauges are part of a vacuum interlock system that will close gate 
valves to protect the turbo-pump if the pressure somehow climbed too high. Typical base pressure 
is about 5 x 1CT 7 Torr. 

1.6.4 Data acquisition 

The primary architecture for handling time-domain diagnostic data consists of 10 CAMAC 
crates controling fast analog-to-digital converters (digitizers) located inside the experimental area, 
which are connected by a fiber optic serial highway to the main control room computer via a GPIB 
card. The data from up to about 120 digitizer channels is stored to disk immediately after each 
shot, with a total read time of about 5 seconds or less, depending on the number of channels being 
stored. A variety of different digitizer modules are used, and typically we run at sampling rate of 
10 to 100 MHz (1 digital sampling every 100 ns or 10 ns). 

The data acquisition system software, with integrated accelerator system controls and data 
viewing was implemented in the form of a custom made user interface designed in Lab View. The 
binary files recorded by Lab View can be read and processed using a separate data analysis program, 
also implemented in Lab View, or they can be converted into delimited text files suitable for reading 
into a spreadsheet application, or other analysis programs such as Igor, Mathematica, Matlab, etc. 

Image and spectral data are handled in a more decentralized manner than the time data. 
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Each camera has its own dedicated acquisition computer, with data and control commands being 
exchanged between them using a proprietary data line (either fiber optic, copper bus, or BNC) and 
a corresponding data acquisition card installed on computer. 



* Exodus As briefly mentioned already, the CTIX-DDT facility has been taken apart, carried by 
truck, forklift and tired backs, and reassembled about 3/4 of a mile away. Somehow we lost DDT 
along the way. Moving to the new experimental facility located in the UCD Hertz Hall building, and 
departing from our historic venue in B435 of LLNL, began in November of 2004, with experimental 
operations ceasing on Sept. 1st 2005, complete vacancy of the premises by Sept. 30th 2005, and 
expected first plasma in February 2006. 
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Chapter 2 

Compact toroid 
magnetic field geometry, 
experiment and theory 

2.1 Force- free magnetic field model 

Let us begin by examining the basic constraints imposed on a self-confined magnetized 
plasma ring existing within a cylindrically symmetric conducting vessel. The coaxial geometry of 
CTIX will require that we apply the following boundary conditions at both the inner and outer 
walls. It is natural to work in a cylindrical coordinate system (r,#, z). 

The time varying part of the magnetic field B of the compact toroid must satisfy the 
boundary conditions 

n-B = 0 nxB = /i 0 K (2.1) 

where K is the surface current density. To satisfy the first condition, B r must vanish at the wall. 
From the second condition we see that B z is nonzero if there is a surface current Kg in the azimuthal 
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direction, and Bg is nonzero if there is a current K z in the axial or z direction. 

Thus measurements of B z and Bg at the wall as a function of axial position and time (z, t) 
will let us know the state of the magnetic field and the current density at the boundary of the 
plasma. 




Figure 2.1: Relationship between surface current and magnetic field components 

Much work has been done on the theory of force-free equilibrium of plasmas in cylindrical, 
toroidal, and coaxial geometries, [ref] The relaxation of a magnetized plasma to a state of minimum 
magnetic energy, while global hclicity is preserved, results in the Woltjer- Taylor force- free equilibrium 
(j x B = V P = 0) given by 

VxB = AB (2.2) 

where the constant A = Hoj/B. Within a conducting vessel of coaxial geometry, and with the added 
assumption of zero axial surface current, K z = 0, an axisymmetric solution exists that is separable 
and is expressible in terms of Bessel functions according to 

B r = BbyCos[A: a (2 ! -^)][Ji(fc r r) + /yi(Arr)] (2.3) 
B z = -B 0 ^sm[k z (z-z 0 ))[J 0 (k r r) + fY 0 (k r r)] (2.4) 
B e = -B 0 sin[fc z (2-2o)][Ji(£vr) + /Fi(M] (2-5) 

where k% + k z = A 2 relates the radial and axial wavenumbers to the magnetic eigenvalue A. 
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Radial coordinate r [cm] 

Figure 2.2: Plot showing the radial dependence of Bessel function solution of the magnetic field. 

To match the first boundary condition, values of k r and / must be found such that 
B r (r,z) — 0 at the inner and outer walls, r = r\ and r = r 2 respectively. The radial term is 
zero when 

/ = -Ji{k r ri)/Y x {k r ri) = -Ji{k r r 2 )/Y 1 (k r r 2 ) 

and so k r must satisfy 

Ji{k r r 1 )Y 1 (k r r 2 ) ~ Ji{k r r 2 )Y 1 {k r r 1 ) = 0 

For the CTIX vessel, n = 3.34 cm and r 2 — 7.62 cm in the acceleration section and so 
k r = 75.18 to -1 . In order for B z = 0 at z — z 0 and z = z 0 + Az, k z = mir/ Az. The axial extent 
of the compact toroid within CTIX tends to be about Az — 1 meter and should be dominated by 
m —1 mode, which implies that k z — 3 m _1 . Using these values the Taylor state eigenvalue should 
be A = 75.24 m _1 . From this result we would also conclude that the ratio of the current density to 
magnetic field is 

j/B = X/fio = 5.98 x 10 7 H 

where H = 1 Henry = C 2 /kg m 2 . 
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This means that for a modest CT surface field of 0.1 Tesla we would have a current density 
of about 6 x 10 6 Amps/ m 2 or 60 kA/cm 2 . This is the rough agreement with of other estimates 
based on experimental data. 

One paper on this general topic is a 1993 paper discussing the Marauder plasma accelerator 
and their experiments [ref]. Their paper tackles the issue of MHD equilibrium for compact toroids, 
and presents a combination of analytical results, simulation and experimental results. Overall it is 
a very useful review of many relevant topics, and is the only published account I have been able to 
find that gives an analytic solution of Vx B — AB in a coaxial geometry. 

Despite their thoroughness, they made a few problematic mistakes in their presentation of 
the theory of force-free equilibria that are worth notice. The least of the errors was a simple typo, 
they gave the magnetic helicity as A x B, when it should be A • B 1 . 

In the same section they made a more confounding mistake by writing down a believable- 
looking yet incorrect solution for the magnetic field of a force- free state Vx B = AB in a coaxial 
geometry. Lastly, they listed a key reference that seemingly does not exist, even after much effort 
investigating a very wide margin of possibility for typographic or factual error. 

As it relates to my present discussion, relying on the correctness of their published formulas 
caused some difficulty when I first attempted to make quantitative comparisons between theory and 
our experiment. The effect is that their formula gives a much larger B z than Bg inside the CT, 
whereas experiment showed that the two components are roughly equal (within a factor of 2). Their 
formula for B z had a coefficient (/c 2 /fc 2 ), but when the curl is taken and field components are 
compared, I found that B z should just depend on (k r /k z ), not squared. The corrected formulas 
listed above (2.3 - 2.5) actually satisfy Vx B = AB, while their counterparts in the 1993 paper 
did not. When this is corrected, the ratio of Bg/B z approaches unity within the CT, basically in 
agreement with experiments. 



See Appendix A for a discussion of Helicity of a force-free state 



22 



A further point of refinement can be made, which will be examined more closely in 
section 2.2, if the bending of the field lines near the probe port hole is taken into account. This 
bending results in an effective depth of the probe measurements that is approximately 0.5 cm into 
the interior of the plasma. The Taylor state model predicts that at this depth, the magnitude of 
Bg should be very close to 1/2 of the value of B z . It also makes a much broader prediction that 
throughout volume of the CT, the ratio Bg/B z should be only a function of radial position simply 
because the z dependence divides out. 

When compared to our experiments, this first prediction that our probes should measure 
Be = \B Z turns out to be reasonable accurate when considering the peak value of B z at the center 
axial position of the CT. On the majority of shots we observe a noticeable dip in the Bg signal that 
is coincident with the CT center, whose amplitude is very close to 1/2 the peak value of B z . My 
interpretation is that this central dip in Bg is due to Taylor state fields. 

However, the behavior of Bg elsewhere in the CT is in complete disagreement with the 
second prediction of the Taylor state model: that the ratio Bg/B z is not a function of z. Our 
measurements show a significant ramping up of Bg, going from zero at the leading edge of the CT 
to about 1/2 to 3/4 of the value of B z (t pea k), by the time the back edge of the CT has passed the 
probe. Then when the probe crosses over into the pushing field, Bg increases until it is near the 
value of the peak of B z . 

This linear ramp-up of Bg inside the CT is not accounted for by the Taylor state, and some 
kind of modification to the simple force-free state is required. 

Many interesting questions remain that involve exactly how the CT magnetic geometry 
varies away from a true force-free Taylor equilibrium when it is being strongly accelerated. We will 
briefly examine one possible modification to this theory, after asking some hard questions about the 
physics that underlies our magnetic probes, in order to get them up to the task of testing these 
various theories. 
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2.2 Magnetic field measurements 

The measurements of the axial and toroidal magnetic field components are done with special 
magnetic probes, each component is measured with its own 1 cm x 1 cm square loop coil of several 
turns. The two loops are wound at right angles to each other through small holes at the tip of a 
long plastic rod (see figure 2.3). 

The magnetic probe operates like any standard Faraday inductive pick-up coil. As the net 
magnetic flux through the loop changes when the magnetized plasma flows around it, a corresponding 
current is induced in the wire that is proportional to the time derivative of the instantaneous magnetic 
field at the loop. This current is put through a passive RF integrator that effectively takes the time 
integral of the signal, and the output signal is digitized. 

This recorded signal has some further numerical corrections applied to it in order to com- 
pensate for an inherent droop effect of the integrator circuit. The final product is a well calibrated 
measure of the magnetic field strength as a function of time at the position of the loop. 

In order to make a quantitative comparison between theory and our experimental measure- 
ments of the magnetic field, we need to split a few technical hairs. In particular, the exact shape of 
the vacuum vessel in the vicinity of the probes will play a critical role in formulating a correct in- 
terpretation of the experimental signals. The results of this analysis will be essential in determining 
the systematic errors involved with our magnetic measurements, needed in later chapters for placing 
constraints on estimated plasma resistivity and fluid velocity fields. 

Each magnetic probe is inserted into a special port that consists of a 0.874 inch hole in 
the stainless steel outer conductor, with a stainless tube (0.874 inch inner diameter) welded to the 
hole and extending radially outward (about 1 inch) to a flange that makes an o-ring vacuum seal 
with a 1/2 inch diameter ceramic tube. The ceramic tube is 6 inches long with one closed end that 
is inserted into the vessel so that the tip of the ceramic tube is approximately flush with the inner 
surface of the outer conductor. The other end of the ceramic tube is open to the air, and a plastic 
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rod with the magnetic probe coils at its end is then inserted fully into the ceramic tube so that the 
coils are as close to the plasma as possible (see figure B.6 in appendix B for more details). 

The ceramic tube acts in three ways: 1) as high voltage insulator to protect the probe 
coils from the ~ 10 kV potential of the outer conductor, 2) as a vacuum seal and vessel wall that 
separates the atmospheric pressure from the 10 -6 Torr vacuum inside CTIX, and 3) as a dielectric 
boundary that lets the magnetic field at the edge of the plasma enter the probe coil and be detected. 

Thus the probe coils are actually inside a finite-sized conducting cylindrical well, approx- 
imately 0.874 inches in diameter and 0.765 inches deep. The boundary conditions are significantly 
different than if the port was not there. As the magnetized plasma flows by a port, it dips into 
the annular gap between the outside of the ceramic tube and the inside of the stainless steel. Even 
if plasma flow effects are ignored, the edge magnetic field will expand to fill the conducting well 
of the port and there will be a gradient in the field as components change to match the boundary 
conditions. So the result of the measurement depends strongly on the precise depth of the probe 
position inside the port well. 

The radial position of the probes have been carefully measured. Here we list the difference 
between the tip of the probe coil and the inner surface of the outer conductor, measured in inches, 

Ar (t 'outer ^ probe tip) 

Ar 57 = -0.151", Ar 91 = +0.0125", Ar li2 = +0.0075". 

Here, we took into account the thickness of the wire used for the coil (0.0625"), so that the "probe 
tip" really corresponds to the top position of inner area of the coil, not the tip of the wire surface. 
Deciding how to split these hairs is not easy, but the small scale of the system we are considering 
requires that we make precise definitions, sometimes without any good intuitive guidance from simple 
idealizations. 

One key characterization that is needed to properly interpret the magnetic measurements 
is the calculation of the effective depth of the experimental magnetic measurements. The idea is 
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that because of the well geometry of the port, the probes are effectively measuring the fields at some 
point in the interior of the plasma, some small but not exactly zero distance away from the outer 
wall. 

This would help reconcile the fact that the analytic Bessel function solution to the Woltjer- 
Taylor force- free equation has zero Bg at the walls, while our probe measurements show it to be 
non-zero. If we accept that the probe is not actually measuring exactly "at the wall" but rather at 
the effective depth into the interior, then the question of consistency is relaxed. The comparison 
between model and experiment would then become a quantitative comparison of the ratio of field 
components at the effective depth. 

A reasonable approximation can made by considering the magnetostatic problem with an 
ideal conductor in the shape of the magnetic probe port. This ignores any plasma dynamics that 
might effect the measurement, but in the case of a low beta plasma, like our CT, the magnetic field 
is boss. And although we consider a magnetostatic problem, this isn't such a bad approximation 
given that the relevant MHD transit timescales are many orders of magnitude slower than the speed- 
of-light transit timescale for the probe port system. Still, we hope to greatly improve upon these 




Coil 



Figure 2.3: Sketch of magnetic probe port and the effect it would have on a uniform field. See 
appendix B for a more detailed drawing of the probe port. 
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preliminary estimates in the future, and in part, I am documenting the details of my work here in 
order to provide a starting point for future work, including some 3-D MHD simulations proposed 
for the near future. 

We want to find the geometric effect on field lines near the port, given that the field is 
uniform at infinity. There will exist some field line that intersects the tip of the probe coil. If we 
follow this field line away from the port sufficiently far, it will level off uniformly. The effective 
depth of the probe would then be defined as the distance from the wall of this field line when it is far 
from the port. Additional calibration factors due to field expansion must also be taken into account. 
The effect is that the measured field is very strongly dependent on the exact radial position of the 
probes. The interior fields of the CT are several times larger than what we would expect without 
these additional calibration factors. After certain brief theoretical considerations have been made, 
we will compare predicted and measured calibration factors. 

In order to tackle these issues I considered a variety of approaches. The mathematical 
details of these methods are presented in Appendix B. However, there are several key results from 
these methods that are relevant to the current topic and should be summarized. 

Because the diameter of the probe coil is relatively small compared with the diameter of 
the port hole, an accurate 2-D approximation can be made for the magnetic fields near the central 
plane of the probe port well. I have found that field expansion effects should be more faithfully 
recreated by a 2-D solution of the right cross-section geometry than a 3-D solution of the wrong 
cross-section geometry. For example, there exists a standard, fully analytic result for the magnetic 
field near a plane with a circular hole in it. This is attractive because it has some features that 
are qualitatively similar to our probe port; however, the key difference is that in this simple case 
the field is allowed to expand out infinitely into to void beneath the plane, instead of just into the 
finite sized well that we actually have. This means that the calibration factor / B f °° — due to field 

\ / probe 

diminishment caused by expansion, would be somewhat less in the hole geometry than for a finite 
well. Also, each field line that enters the hole will go deeper than it will in the case of a finite well. 



27 

Thus for a given probe position, the resulting effective depth will be noticeably larger for the simple 
hole than it would be for the finite well. 

There is also one significant practical simplification that we should use to our advantage. 
The magnetic field vectors themselves are not really important to find the effective depth and 
expansion factor. Really we just need to know the shape of the flux surfaces on which the magnetic 
field lines reside. The simplifications arises from the fact that the magnetostatic potential 4>m is the 
harmonic dual of the electrostatic potential cpE- To go from a potential function to its harmonic dual, 
we just have to interchange the role of equipotentials and flux surfaces, and the resulting potential 
will still satisfy the Laplace equation, but with the dual boundary condition, [ref] 

The magnetostatic potential satisfies the Neumann boundary condition at the surface of the 
conductor, while the electrostatic potential for the same conductor geometry satisfies the Dirichlet 
boundary condition. Therefore, we can find the flux surfaces of the magnetic field by actually 
computing the equipotential surfaces of the dual electrostatic potential problem. This is an advantage 
because Dirchelet boundary conditions are intrinsically simpler in that you don't have to evaluate 
derivatives to apply them, and they lend themselves to numerical solution by iterative relaxation 
methods. 

Once the flux surfaces have been evaluated using your method of choice, the effective depth 
can be found by tracking field lines away from the well and then finding how high above the wall 
they are after going several hole-radii away from the well. Field expansion ratios can be calculated 
by comparing the distance between two neighboring field lines, when they are far from the well, vs 
when they are at the center of the well. 

Here we are using the fact that the magnitude of the magnetic field is proportional to 
the density of field lines, which gives us the ratio of the far field to its corresponding value of the 
distorted near field. To measure the average flux through the probe coil, we just choose one field 
line that passes through the upper tip of the coil, and a second that passes through the lower tip of 
the coil. If the distance between these two field lines far from the probe well is Sz and the height of 
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the probe coil is Az co u then the calibration factor due to field expansion is C exp = z £° il . 

Following the procedure outlined above, I have implemented a simple 2-D Gauss-Seidel 
method to find the flux surfaces of the magnetic field in the region near the probe port, and then 1 
evaluated the calibration factor and effective depth as a function of probe position and hole diameter, 
[ref jackson 1.13 pg 47-50] (see appendix B for details) 
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Figure 2.4: 2-D Relaxation simulation result of magnetic field near a conducting port well. Image 
on left uses a cyclic gray-scale to visualize flux surfaces. Image on right represents the magnetic 
field vectors using a color-wheel type coloring of the complex plane, which corresponds to the (y, z) 
directions in physical space. Red indicates a vector pointing directly in the -y direction, while other 
colors represent vectors at other directions as indicated by the color legend in top left. Vector 
magnitude is represented by a color's brightness, with black for the zero vector. 



The result of the magnetic field computation is shown in fig 2.4, using two different visu- 
alizations. On the left, the potential function has been displayed in a cyclic gray-scale to visually 
identify level sets of the potential, which correspond to magnetic flux surfaces (and also field lines). 
On the right, the magnetic vector field has been displayed in a novel color-based vector representa- 
tion, which is described in appendix E. Also, a set of field lines have been overlayed to guide the 
eye. 

This is a nice result, but we should try test its predictions against some real probe data. 
This can be done using an independent consideration regarding actual magnetic measurements that 
leads to a relative calibration between the primary set of magnetic probes. Our goal is to test 
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Figure 2.5: Probe position dependent trends in the expansion factor (left) and the effective depth of 
measurement (right) for magnetic field expansion into full volume of well. Position zero corresponds 
to the probe coil tip being flush with the inner edge of the port hole. 



Case 1. 


Calibration factors 


C57 


3.67827 


C91 


2.52267 


Cl42 


2.32853 


Effective depth of probes (inches) 


P57 


0.11045" 


P91 


0.17906" 


Pl42 


0.18732" 



Table 2.1: Calibration factor (C) due to field expansion and effective depth (p) for each magnetic 
probe, using the results after 2000 iterations in the case of the field expanding into the full volume 
of the probe well. 
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if the relative difference in probe tip position is agreement with the resulting difference in probe 
measurements. 

The basis of this comparison relies on the observation that for late times (t > 14/xs) during 
the acceleration process, the Bg pushing field behind the CT relaxes toward an equilibrium state 
that is nearly uniform as a function of z. 

This effect is clearly seen in the time signals for the toroidal magnetic field component 
Bg, and it very closely agrees with the simple kinematic model of a uniform magnetic slab that is 
decreasing in field strength due only to its uniform expansion. 




n 1 i 1 i 1 i ' i ' i ' i ' i ' i ' i ' i ' i ' i ' r~ 

2 4 6 8 10 12 14 16 18 20 22 24 26 28|jsec 

time (microseconds) 

Figure 2.6: B z and Bg signals at z = 57 cm as a function of time. The fit to the expansion model 
is shown as a black curve. 



The formula used to fit the expansion of Bg(t) is 

f(t) = foAz ° (2 6) 

JyJ Az 0 + Av(t-t a ) + h 2 Aa(t~t 0 y y ' J 

We hold constant to, which is the starting time for the region of the curve that is being fitted, 
and also hold constant the initial signal level fo = f(to), and the initial axial extent of the slab 
Azo = z pro b e . The free parameters for the fit are the difference in velocity and acceleration between 
the front and back of the slab, Av and Aa. 
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In order to compare total calibration factors we just need to look at the ratio of the post-CT 
Bg signals as measured by different probes at the same time. If these ratios equal 1, then there is 
no problem with the relative calibration between probes. Any difference between the probes can be 
compared to predictions of the field expansion analysis given above. 

When we look at the actual Bg data, there is a small but systematic discrepancy between 
the measurements at the three locations. This could mean that there exists a highly reproducible 
non-uniformity of Bg(z,t) that depends on z in an exactly fixed way, or that previous calibration 
factors for the probes were slightly off by differing amounts. 

The average value (over 60 shots) of ratios of the three possible pairs of probes are: 

Be [571 B e [571 B e \91] , . 

4^ = 1-02529 ±0.10296, = 1.21117 ± 0.17172, ^ = 1.19293. (2.7) 

This shows that Bq [57] and Be [91] are nearly balanced, while £?#[142] is about 20% di- 
minished relative to the other two probes. Comparison with the data in table 2.1, shows a clear 
discrepency. If the only unaccounted relative difference between the probes was due to the depth- 
dependent field expansion effect, then we would expect the probes at 91 and 142 to be almost exactly 
balanced (because they are inserted at nearly the same depth), while the probe at 57 cm would have 
a reduced signal compared to the other two (because is is more recessed into the well). Using the 
computed values for the field expansion calibration factors we would expect that the signal ratios 
would be 

|4^4 = 0.68583, 0.63305, *U 0.91903. (2.8) 

B e [91] B 0 [142] B e [U2] v ; 

The disagreement between (2.7) and (2.8) could be due to a variety of reasons. There could be errors 

in my 2-D approximations, or in the accuracy of the probe calibrations the we have been using. Or, 

the existing probe calibration may be actually better than anticipated and it has already taken these 

relative scalings into account, implicity by the nature of the measurements. This would bring the 

ratios closer to unity, but there is still some of room for uncertainty about the matter. 

An alternative way of reconciling this issue is to consider that during the CTIX discharge 
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there might be a stationary plasma surrounding the probe, acting as conductive filling that excludes 
the magnetic field of the passing CT. The field of the compact torus would only expand into the 
interior of the ceramic tube itself, which is the same for all the probes. In this case, the details 
of precisely how far the probe is recessed into, or protruding out of a larger port hole, would be 
unimportant. Effectively, all the probes would have the same exact geometry for the conducting well; 
one made of the plasma surrounding the ceramic tube, not the stainless steel vessel. To calculate the 
effective depth and field expansion calibration factors, we would use a well with the same diameter as 
the ceramic (0.5") and measure the conditions at a position that is recessed into the well a distance 
equal to the thickness of the ceramic wall at the probe tip. Figure 2.7 shows the results of the 
analysis in this case with stagnant plasma in the annular gap around the probe (call this Case 2). 

Case 2. Magnetic field calculations 
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Figure 2.7: Probe position dependent trends in expansion factor (left) and effective depth of mea- 
surement (right) for magnetic field expansion into only partial volume of well due to field exclusion 
by stagnant plasma. Position zero corresponds to the outer surface of the ceramic probe sleeve. 



Over the course of these considerations we have found that value of the effective depth of 
any of the probes should be between 0.079" and 0.187" (2 mm to 4.7 mm), while the additional 
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Case 2. 



Calibration factor 



(^expansion 4.20797 



Effective depth of probes (inches) 



Peff 



0.07919" 



Table 2.2: Calibration factor (C) due to field expansion and effective depth (p), which are the same 
for all the magnetic probe in the case of the field expanding only into the smaller volume of the 
ceramic tube. We use the results after 2000 iterations. 



calibration factor due to field expansion should be between 2.3 and 4.2. However, the precise 
values within these ranges are still too indeterminate given the existing data from experiment and 
simulations. Direct experimental measurements of these constants should be done soon, and this 
could be complemented by some better 3-D simulation, possibly including plasma effects. For now, 
we can let this range of values serve as error bars on our subsequent estimates. 

As it relates to our primary question of how the compact toroid magnetic field can be 
modeled, the effective depth is the critical variable for comparison with experiment. Our estimates 
for the effective depth can now be used to quantitatively compare the ratio -g 2 - between theory and 
experimental measurements. Using the Bessel function solution of Vx B = AB in a coaxial geometry 
with n = 3.34 cm and r 2 = 7.62 cm we can evaluate the magnetic field components at any specified 
depth. We find that 



(p e ff = 2mm) 



El 

B, 



= 0.1379, {peg = 4.7mm) 



= 0.3818 



(2.9) 



Conclusions Based on our preliminary observation that most CTs have a central dip in Bg equal 
to about 1/2 the peak of B z , we can now say with some certainty that the Bessel function model 
is in moderate agreement with our measurements, provided that the true effective depth is at, or 
above, the upper end or the believable range. However, if the physically real p e g is as small as 2 
mm, then this is not the case; the measured Bg would be about 4 times too large to be accounted 
for by the Bessel function model. On the other hand, if we could stretch our errors enough to let 



34 

p e ff — 6 mm be considered reasonable, then we would have a perfect match between theory and 
experiment. 

Further progress can be made toward quantifying this comparison. In addition to the 
improvements that are possible with understanding the magnetic expansion effects in the vicinity of 
the well, we can also better quantify this "dip-to-peak" ratio in actual data and gather the statistics 
of this quantity over a large data set. Clearly, this is a task that would be best postponed until there 
ample time to do it properly (when some future student feels compelled to take up the cause). In 
the meanwhile, we can trust that there is truth in idea of the CT existing in some sort of minimum 
energy state, and rely upon this idea as an approximation if it helps to make difficult problem more 
tractable. 

Next, we will introduce a very useful "quasi-tomographic" method, which I call Lagrangian 
interpolation, whose purpose is to squeeze as much information out of sequence of three probe 
measurement as is humanly possible. We will use this method on several occasions throughout this 
study, for both visualization and deduction. 

2.3 Lagrangian interpolation of probe signals 

In order to bridge the gap between theory and experiment there is much progress that can 
be made simply by improving our understanding of existing experimental data. Here we present 
a method of reconstructing the approximate spatial dependence of the magnetic field within the 
plasma, using the data collected from three magnetic field probes as the plasma rushes past them 
during the acceleration process. The result is an accurate estimation of the magnetic components 
B z and Bg as a function of axial position at each instant in time. These lend themselves to being 
turned into nice little movies of the CT going down the accelerator, making it easy to visualize the 
spatial structure of the magnetic fields. 2 

More importantly, once we know how these magnetic field components depend on (z, t) we 

2 Magnetic field movies can be downloaded from our website http://ctix.das.ucdavis.edu 
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can take derivatives with respect to space and time, allowing the full repertoire of electrodynamics 
to be applied. We can then develop certain quantities that have been derived from our experimental 
data, and compare them to theoretical predictions. 

I call this method Lagrangian interpolation because it follows the trajectories of fluid ele- 
ments through the system, and performs an interpolation of the field quantity using the measured 
values at the points where a trajectory crosses the probe locations, doing so along the path of each 
trajectory. This interpolation relies on first making an estimation of the velocity field throughout 
the system, usually based on the apparent time-of-flight kinematics of the magnetic pulses. A more 
detailed account of the method is presented in appendix C. Here we will examine those specific 
results relevant to resolving the issue of the magnetic geometry of the compact toroid. 

One way to estimate the flow field on CTIX is to track one very distinct feature as it travels 
down the pipe, and then apply the same relative kinematics to the rest of the fluid elements in the 
system. We have found that the most stable such feature on the CTIX waveforms is the back edge of 
CT. This is the junction point between the CT and the pushing field, and we can reliably define the 
arrival of this point as being the time when the B z signal of a given probe crosses the half-maximum 
level on the trailing edge. The high accuracy of this method is due to the very reproducible and 
cleanly steep slope that occurs at the CT back edge on virtually every shot. 

In this way we can find the times that the CT back edge passes by the three accelerator 
probes, £57, tgi, £142- These immediately allow us to find the average velocities due to the transit 
times between the probes 



91 - 57 



142-91 



V57-91 = 



tai — t 57 



Vgi-142 = 



^142 — *91 



Then from the change in velocity we can find the average acceleration of the back edge 



a = 



2(V91-142 — V57-9l) 
(tl42 — £57) 



We can subsequently infer an initial velocity and initial position, 



Vo — V57_9i — a^ 



(t 57 + tgi) 

2 
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1 _ 2 

z a = 57 - v 0 t57 - o a i 5 . 



57 



This is based on the assumption that the CT back moves with constant acceleration 

zcr{t) = zo - v 0 t - - a t 2 . 

For positive values of acceleration there will always be a minimum z- position, z m in, that the CT 
never goes below, according to this simple description. 



This is also the point on the trajectory where v(t) = 0. The problem is that for early times, before 
the fluid element has reached z m i n according to this description taken literally, it would have a 
negative velocity, which is unphysical for our system. Instead we will only apply these uniformly 
accelerated models to the region above z m i n , and during times of positive velocity. 

If we work with a simple model in which all the fluid elements have the same acceleration 
throughout time and space, they will all have the same z m i ni but they would pass through it at 
different times. Based on this, we arrive at a velocity field of the form 



Notice that this depends on space but not on time. We can think of this like a Taylor series 
approximation which is most accurate near the central point, and decreases in accuracy as you move 
away. In this same way, this velocity field will be most accurate in the neighborhood of the CT back 
edge. But long after the CT has passes by, the acceleration dynamics will certainly change, and this 
approximation will no longer be a good one. 

Some modification to this velocity field must be done to make it fit certain basic physical 
constraints. For the case where z m i n > 0 we clearly need a different way to handle the velocity in the 
region z < z min . Since this will always occur in the region before the first probe at z = 57 cm, where 
there is a scarcity of real velocity information, the most pragmatic approach is to set the velocity 
equal to a constant value that matches the accelerated velocity curve at some point Zfu > z min . For 




V(z,t) = ^2 a (Z - Z min ) 



(2.10) 
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instance, we have found good results by defining the position of the fit Zfu = min(|,2 mi „, 57cm), 
then defining the velocity in a piece-wise fashion. 



y/2 a (z - z min ) if z > Zf it 
v(--.l)={ (2.11) 

y/2 a (zfu - z min ) if z < Zf it 
An additional modification is needed because we begin recording the signal about 12 \is 
before the shot begins, and then there is a period of about 4 [is before the accelerator back fires, 
during which time the velocity field is driven only by the formation bank. This leads to a natural 
division of the time domain into three parts. 1) Before the formation back fires, the velocity should 
be zero everywhere. 2) Once the formation current begins to flow, the velocity field should be 
determined by the kinematics of the leading edge of the CT, where we detect the first early rise 
above the zero level of B z . 3) Lastly, once the accelerator current turns on we can use a velocity 
field like (2.11), based on the motion of the CT back edge. 

This is not the only possible velocity field that can be extrapolated from the kinematics 
of the back edge of the CT. However, when it is applied to actual data it seems to give the best 
agreement with our expectations of how the magnetic fields should behave globally. For instance, 
velocity models that go to zero as z — > 0 will produce a quite unphysical "precognition" effect near 
the origin, where there is a build up of fluid trajectories that results in an axially compressed version 
of the future wave pattern that seems to be "waiting in the wings" until it hears its cue and makes 
its debut passing by the first probe and out onto the theater of our observation. Instead of this, 
some moderately large initial velocity is needed to be in agreement with the fact that there is a high 
rate of magnetic flux input from the external circuit, as well as an ongoing formation of new plasma 
by ionization of the steady flux of neutrals from the slowly closing gas valve. In this way, velocity 
fields like (2.11) are good because they are qualitatively consistent with these facts, and therefore 
prevent the fluid elements from "hanging out backstage" too much. 

Once we have a reasonable estimation of the fluid velocity field, we can interpolate the 
real probe signals across the trajectories, filling the whole (z,t) domain with values of the field 
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Figure 2.8: The estimated velocity field on the (z,t) domain. The trajectory of the CT back is shown 
as a black curve. 

components. We can then sample this 2-D array of interpolated values along slices of constant time, 
and see the shape of the fields as a function of z, at any instant of time. Typically we performed 
this interpolation using 200 points in the z direction, one point per cm, and 1024 points in the time 
direction at a rate of 0.1 \xs per point, taken from the scaling of the raw time signals. 

Shown in figures 2.9 to 2.14 on the following few pages, is the resulting interpolation 
of the z and 0 components of the magnetic field at several instants in time. Before computing the 
interpolation, we applied the calibration factors discussed in the previous section to the probe signals. 
The safest approach is to use the relative calibration determined from comparison of the experimental 
ratios of Bg in the expanding pushing field, and combine that with an absolute expansion factor 
determined from the magnetostatic simulation. These factors had values of: 

4 2 4 2 4 2 

C 57 = — = 3.4682, C 9 i = — ^ — = 3.5207, C U2 = — = 4.2 
57 1.211 ' 1.1929 ' 1.0 

As can be clearly seen in these figures, the axial dependence of the magnetic fields is 
clearly more complicated than the simplistic sin(fc z (z — zq)) and cos(fc z (z — zq)) of the force-free 
Bessel function model. However, there is some evidence to support that the radial dependence of 
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Figure 2.9: Jig and B z as a function of axial position at time t = 7.5 (is for shot 59842. 
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Figure 2.10: Tig and H z as a function of axial position at time t = 8.5 (is for shot 59842. 



the fields is reasonably well described by the Bessel function model. We may combine the two if we 
assume that the field is separable in terms of radial and axial dependence. Doing so will at least 
allow us to make a comparison between experiment and a minimally modified version of the Taylor 
state equilibrium. 

In the following discussion we will use this Bessel function combination so frequently that 
it is notationally useful to define 

Fi{x) = [Ji(x) + fYi(x)}, ie {0,1,2,...}. 
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Figure 2.11: B# and B z as a function of axial position at time t = 9.5 [is for shot 59842. 




Figure 2.12: Be and B 2 as a function of axial position at time t = 10.5 fis for shot 59842. 



This merger requires having a good estimation of the effective depth of the probe mea- 
surements. The z-component of the field would then vary throughout the full domain of the system 
according to 

B z (r,z,t)=B z (z,t) /^ krr \ (2.12) 
F 0 (k r r eff ) 

where B z (z,t) is the interpolated probe signal as a function of axial position and time, and r e // is 
the effective radial location of the probe measurement. All other terms are the same as have been 
used previously. The form of this is necessary to match the probe signals B z (r e ff,z,t) — B z (z,t). 
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Figure 2.13: Bg and B z as a function of axial position at time t = 11.5 fis for shot 59842. 




Figure 2.14: Be and B z as a function of axial position at time t = 12.5 fis for shot 59842. 



Likewise, the ^-component has the form 

Bg(r,z,t)=B e (z,t) f^ krr \ (2.13) 
l<i(K r r e ff) 

We can use this radial dependence and the non-existence of magnetic monopoles (V- B = 0) 
to evaluate the r-component of B. 

dz F 0 (k r r ef f) 



dB z (z,t) 1 F^Kr) 
dz k r F 0 (k r r e ff) 



(2.15) 
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Figure 2.15: plot ol B r , B Zl Btheta at a single time. 



We can evaluate Ampere's equation Vx B = /ioj in cylindrical coordinates with assumed 



axisymmetry (J^ = 0). The resulting components of j are: 



1 dB a 



Mo £> z 

1 / dB r dB z 
38 = — 



3z 



fiQ \ dz dr 
1 d rBg 



(2.16) 
(2.17) 
(2.18) 



Hor dr 

Notice that both j r and j z depend only on Bg, while jg depends on both B r and B z . 

We find that at the effective depth of the probe, the following generalization appears to 
hold. The radial current density has a nearly constant proportionality with the radial magnetic field 
component, j r oc B r , in agreement with j x B — 0. Yet the other two components appear to be 
cross- wise related, j z oc Bg, and jg oc — B z . This leads to a significant deviation away from a force 
free state. 

The proportionality between j r and B r , is interesting and significant because they are 
derived from independent experimental measurements. The radial current density j r comes from 
Bg(z,t), while B r comes from B z [z, t) using the assumption of the Bessel function radial dependence 
only to fix the magnitude of the signal, its z-dependence is determined purely from B z (z, t). So this 
is a true symmetry inherent in the actual magnetic field configuration, and not some artifact of the 
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Figure 2.16: Trio of overlay plots, j r and B ri j z and B Zl jtheta and B t heta at a single time. 



reconstruction. 

On the other hand, the relations between the z and 9 components arises from an actual 
lack of symmetry between the z and 9 magnetic fields. We must conclude that j x B ^0 within the 
CT, and that this departure has a noticeable effect on the magnetic equilibrium. This effect is really 
nothing more than what we have known all along, that the CT is being dramatically accelerated by 
a large rail-gun Lorentz force density equal to j x B. The uncertainty up to this point has always 
been that we were unsure of exactly where in the system j x B has nonzero values, and to what 
extent we can apply the standard model of force-free equilibrium. 
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2.4 The Hall term, j x B 



Based on the Lagrangian interpolation of the probe signals we can calculate the j x B vector 



throughout the system. Shown in figure 2.17 is the magnitude of j x B, as well as the resulting angle 



between j and B, both calculated from the data taken on shot 59842. 




Figure 2.17: The magnitude of j x B and the angle between j and B, shown as a function of axial 
position at time t = 7.8 (is for shot 59842. 



The primary qualitative observation is that j x B is large throughout the whole CT, and 
not limited to just the interface with the pushing field. The plasma may be in some relaxed minimum 
energy state, but is not anything close to force-free. We see that within the center of the CT, the 
current is approximately at right angles to the magnetic field, while at the edges the angle drops 
down to about 45°. 

We can break down the j x B vector into (r, 9, z) components, and consider their individual 

effects. 

We see that most of j x B is due to a component in the negative radial direction. This 
should be in agreement with the standard MHD equilibrium condition j x B = VP, where P is 
the hydrodynamic pressure of the plasma. Due to the large radial gradients in density, the pressure 
gradient at the outer wall should also be directed inward, in the -r direction. Exact quantitative 
comparison between this value of (j x B) r and reasonable estimations of (V P) r , still remain to be 
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Figure 2.18: The radial, toroidal, and axial components of the j x B vector, shown as a function of 
axial position at time t = 7.8 /j,s for shot 59842. 



done. It is hopeful that their equivalence may serve as an independent constraint on the calibration 
factors used in the magnetic measurements. 

We may also be able to relate the volume integral of (j x B) z with the net force on the CT 
that results in its acceleration in the axial direction. Ball-park estimations of these quantities show 
promising agreement. It is interesting to note that (j x B) z has a strong positive peak at the back 
edge of the CT, while it has a lesser negative peak at the front edge. This is what is responsible for 
the CT maintaining its shape as it travels. The net axial force on the CT should correspond to the 
pressure differential across the CT. 

The most surprising result is the effect of a nonzero (j x B)#, which should result in toroidal 
spin-up of plasma. The apparent front-to-back sign reversal of (j x B)g would imply that the front 
edge is counter-rotating relative to the back edge, with a strong axial shear in the toroidal velocity 
field. The fact that (j x B)g is comparable in magnitude to (j x B) z would seem to imply that the 
toroidal component of the velocity may be as large as the axial velocity of the CT (5 to 30 cm//xs). 

Anticipating this possibility, we did make an attempt to measure this toroidal effect via 
Doppler measurements of the He II ion line at 4686 A, using the high resolution spectrometer 
described in chapter 6. We did this by acquiring the light from the z = 57 cm optical ports at 60° 
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and 300°, using multiple fibers inputs imaged on the same intensified CCD at the dispersion plane 
of the spectrometer. One fiber collected light in the + 9 direction, while a second fiber collected 
light in the - 9 direction. Any difference in the center wavelengths of the line profiles should be due 
to a net toroidal rotation of the ions. The preliminary measurements of toroidal rotation were null, 
vg < lcm/fis, but this may have been due to the relatively long exposure times (2 to 5 /us) needed to 
overcome problems with low light levels in this spectrometer configuration. Long exposures would 
result in a blending of the positive Doppler shift at the front edge of the CT with the negative 
shift due to the back edge. On average, over a few fis, there would be zero net toroidal rotation 
observed. In order to accurately measure a differential Doppler shift like that predicted by this result 
for (j x B)g, the measurements would need to have exposure times of less than 0.5 /is, preferably as 
short as 0.1 fis in order to resolve a change in velocity across the length of the CT. 

Now going back to our original question, "How do we describe the magnetic equilibrium 
of a compact toroid?", we find that all of these considerations lend their support to the validity 
of a more recent theory that is a generalization of Taylor's relaxation idea, but includes the effect 
of a non-zero Hall term (j x B) in the formula for resistivity of the plasma (Ohm's Law). This is 
ultimately part of an entire field of study often called Hall Magnetohydrodynamics (HMHD), which 
can yield quite different descriptions of the same phenomena treated by regular MHD (which ignores 
the Hall term), simply due to this change in the character of the resistivity of the plasma. 

The idea first proposed by Turner [ref], is that magnetic relaxation can still occur, but when 
the Hall term is included in Ohm's law, a different quantity that he calls the hybrid helicity becomes 
an invariant of the system. The conservation of the both the magnetic helicity and the hybrid 
helicity will constrain the magnetic field so that it can never fall all the way down to the Taylor 
state, but instead reach a different equilibrium which is not force-free. There are actually several 
distinct generalizations to the Taylor relaxation theory, but Turner's idea is overall the simplest, 
and most applicable the CTIX plasma. 

The derivation of an equation for this HMHD minimum energy state will require us to 
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consider the velocity field on equal footing with the magnetic field. In fact we consider the vorticity 
of the flow, and the cyclotron-normalized magnetic field. 

eB 

uo = Vx v, uj c = — 
mc 

We can then define a hybrid vorticity according to 

fi = u + UJ C 

and this evolves in time via 

on 

— = Vx(vxtt) 
If we note that Jl is expressible in terms of a vector potential V 

fi = VxV 

and consider a solenoidal form of V, (i.e., V- V = 0), we can see that 

eA 

V = v+ — + VX 

mc 

where A is a particular instance of the magnetic vector potential, and \ is an arbitrary scalar 
function that satisfies Laplace's equation V 2 \ = 0. The hybrid helicity is defined as the volume 
integral of the dot product between fi and V 

Hhybrid = I V n d 3 x 



It can be shown that this is an invariant of an ideal MHD fluid with a non-zero Hall term. 

A two-parameter Lagrange multiplier analysis can be performed to solve the variation 
problem in which total energy is minimized while the system is constrained by both the magnetic 
and hybrid helicities being held constant. The result is that each of the Lagrange multipliers serve 
as eignevalues of a pair of (coupled) minimum energy equations, 

Vxa) c = (Ai + A 2 H + ~\ 2 oj (2.19) 

Vx(uj + uj c ) = ^- (2.20) 
A2 
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where Ai,A 2 are the multipliers after having been normalized, A, = ;j-Aj. The solution to this 
system of equations for w and uj c is what we call the Hall minimum energy state. 

In the future, we would like to find an exact solution to these equations for the coaxial 
geometry and compare them to our experiment. Additionally, the main features of this theory 
provide good motivation to attempt improved measurements of the velocity and magnetic field 
profiles by experimental means. 
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Chapter 3 

Determination of plasma resistivity 

Of all the properties of the compact toroid plasma, it is the electrical resistivity that 
provides the most fundamental characterization of its nature as a magnetofiuid. Resistivity measures 
a material's intrinsic tendency to block the flow of electric current, by way of some energy dissipation 
mechanism. For a conductive fluid, the resistivity determines how electric and magnetic will couple 
to the velocity field of the fluid. The resistivity also plays a central role in determining the value 
of the magnetic Reynolds number used to characterize the nature of turbulence in a magnetofiuid, 
as discussed in chapter 7. We need to consider a few questions. Yet despite this important role, 
the resistivity and related magnetic diffusion coefficients remain poorly understood. For most of our 
estimates we will rely on standard formulas, despite their limited applicability. We will also explore 
some methods to improve our approximations. 

3.1 Collisional Spitzer resistivity 

One of the longest standing results in in plasma physics is Spitzer's formula for the resis- 
tivity of a magnetized plasma. It is derived by analyzing the dissipative energy loss (drag) due to 
long range collisions mediated by the electrostatic Coulomb force between charged particles. The 
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resulting formula is dominantly dependent on electron temperature T e , 

t] ± = 1.03 x 10~ 4 77 3/2 Z In A [O • m] (3.1) 

Where Z is the effective charge state of the plasma, which is very nearly 1 for our hydrogen CT's, 
and the term In A is the called the Coulomb logarithm, which is defined as the natural log of the 
ratio of the maximum interaction distance (r max ) over the minimum interaction distance (r m i n , also 
called the impact parameter), averaged over a Maxwellian distribution, 

InA = In(^^J ~ ln(127m e A?,) = 24 - In {n\ /2 T" 1 ) for T e < 10 eV (3.2) 

V ( r min) / 

The maximum and minimum interaction distances are related to the Debye length, \d of 
the plasma, which is the length over which electric fields are shielded out by free charges in the 
plasma. Xd = \J kT / 47m e e 2 

The presence of a magnetic field acts to impede the flow of charged particles (and hence elec- 
tric current) in the direction perpendicular to the field. This is manifested as directional dependence 
on the resistivity, with the plasma being more resistive in the cross field direction. Under normal 
conditions there is a simple proportionality between the transverse and parallel values, i]± = 1.96 rj\\ . 

The underlying assumptions that make this derivation possible are: 1) The plasma is in an 
steady state (all quantities are constant in time). 2) The magnetic field is uniform in space. 3) The 
electrons are well magnetized, i.e., the gyroradius is less than the mean free path p e <C A m f p . 4) The 
electron motion is dominated by collisions, A m f p <C ^system- 5) The electric field is small enough 
that only small departures away from Maxwellian electron distributions will occur, i?A m f p < T e . 6) 
Collisions between electrons and neutrals can be ignored. 

We can consider the validity of each of these assumptions. The steady state assumption 
(1) is rather hard to apply in any global sense to our experiment given the rapid acceleration that is 
occurring. However, on the scale of a microsecond or two, and from the perspective of a frame that 
is comoving with the plasma, there is some local steadiness to the plasma's behavior, and this will 
be good enough for us. With regards to (2), clearly the magnetic field inside the CT is not uniform, 
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Figure 3.1: The transverse Spitzer resistivity and its dependence on electron temperature. 
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Figure 3.2: The weak but non- vanishing dependence of transverse Spitzer resistivity on plasma 
density. We see that rj± only changes by less than a factor of 2, when n e varies over 4 orders of 
magnitude, due to the weak dependence of the Coulomb logarithm on density. 



and under many circumstances the finite gradients in B will create radical departures away from 
the conditions needed for the Spitzer model to work. Such as the case of magnetic reconnection 
(discussed in chapter 4), we will need to consider more comprehensive models of the resistivity. As 
mentioned in the introduction (table 1.1) the electron gyroradius is at most a few millimeters, while 
the mean free path at typical parameters can be much longer, 3 cm < A m f p < 1 m. This ensures 
that assumption (3) will hold. 

For the CT plasma, the assumption (4) that the electron drag is dominated by collisions 
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Figure 3.3: The electron mean free path as a function of electron temperature and density. 



is just on the edge of validity. Essentially, we can regard a plasma as being dominated by collisions 
when the mean free path is significantly smaller that the size of the system. By inspecting figure 
3.3, we can see that for compact toroids that fall toward the colder, denser end of the spectrum, the 
collisionality assumption will hold well. But for compact toroids that happen to be relatively hot 
and diffuse, the mean free path can easily exceed 1 meter. 

Assumption (5) depends on the electric field within the plasma, and it uncertain exactly 
how large this really is on CTIX. This may hold over most of the system, yet be violated in localized 
pockets where E is large, such as near a reconnective current sheet. The validity of assumption (6) 
is also conditional. When interacting with a dense neutral gas puff, electron-neutral collisions might 
become important. Otherwise it is a fairly safe assumption to make. 
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Overall, the questionable validity of the Spitzer model on CTIX implies that its resistivity 
formula is best use as a lower bound for the resistivity of the plasma. When the full menagerie of 
phenomena are taken into account, the effective resistivity can only be increased about this standard 
value. 

3.2 Theory of magnetic analysis 

The results established in chapter 2 have given us the ability to reconstruct the approximate 
spatial and temporal dependence of the magnetic field within the CTIX accelerator. This can now 
be put to use as a tool to explore some fundamental physics. 

How does the magnetic field evolve in space and time? How does that evolution constrain 
the effective electrical resistivity of the plasma? How does that measured value of resistivity compare 
with standard theoretical predictions. If the magnetic field and plasma velocity field are known 
functions of space and time then, in principle, the effective resistivity of the plasma can be deduced 
using the MHD evolution equation for B 

<9 t B — Vx(v x B) = — V 2 B (3.3) 
Mo 

Where r] is the electrical resistivity of the fluid. Alternatively, we could more directly find the 
resistivity using a generalized Ohm's law. 1 

E + v x B = 77J (3.4) 

It is clear that these are equivalent to each other since the curl of (3.4) equals (3.3). 

To find the resistivity we would need to evaluate the given combinations of the known 
quantities, being careful not to mix-up vector components, and solve for rj. 

But first we must take notice of the vector nature of the above equations. In truth 77 is not 
really just a number. In the presence of a strong magnetic field its true properties are best described 

1 For the simplicity of presentation we are temporarily neglecting the Hall term e ^ (j X B), but wc can add it 
back in at any point. 
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as a tensor that takes into account the difference of plasma conductivity in the direction parallel to 
the magnetic field, as compared to its value in the direction perpendicular to the field. This tensor 
can be expressed in terms of just two scalars, the parallel resistivity r)\\ , and the transverse resistivity 

To describe this, let us choose a coordinate system {ei,e 2 ,e 3 } that is locally aligned with 
the magnetic field vector, 2 in such a way that B = £>e 3 = (0, 0, B). The subsequent choice of the 
orientation of {e l7 e 2 } is completely arbitrary, and can be made however is most convenient. 

In this coordinate system Ohm's law becomes 

E3 = V\\h, E 2 - v\B = t]± j 2 , E 1 +v 2 B = i]± j 1 (3.5) 

If we define an effective electric field E = E + vxBin order to simplify our notation, then 
we have the standard form for the resistivity tensor. 
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One practical difficulty is that we really need to write our equations with respect to the 
laboratory- framed cylindrical coordinates (r, 8, z), in terms of which all our measurements have 
been made. The most straight-forward method is to use Eulcr angles (a, (3, 7) (with standard x- 
convention), [ref] where j3 is the angle in space between the local direction of the magnetic field (e 3 ), 
and the z-axis of CTIX. The angles (a,j) are determined by your previous choice of {ei,e2}. 

Thus, the orientation of the local magnetic field relative to the lab-frame determines a 
transformation (r, 6, z) — > (ei, e 2 , e 3 ), which we formulate in terms of the Eulcr angles. 

For example, if we are able to measure in the lab the (r, 8, z) components of the effective 

2 In general, there may not exist any one global definition that provides a single continuous coordinate system that 
covers the entire volume without being multivalued somewhere. To deal with this problem in a fully rigorous way, a 
set of many coordinate systems might be needed that cover the space in overlapping patches, along with the functions 
that smoothly relate them to each other (i.e., the calculus of manifolds). This level of rigor is clearly beyond the 
scope of this discussion, and we will simply keep in mind the local nature of the above transformations as we try to 
get some testable results. 
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electric field, we can apply a transform to get the magnetically aligned components of E, 
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cos 7 cos a — cos (3 sin a sin 7 cos 7 sin a + cos (3 cos a sin 7 sin 7 sin /3 
- sin 7 cos a — cos (3 sin a cos 7 — sin 7 sin a + cos (3 cos a cos 7 cos 7 sin (3 
sin /? sin a — sin (3 cos a cos (3 
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(3.7) 



and then use (3.6) to find the resitivity. 

In the regions near the wall where B r is negligible, especially at the peak in B z and in the 
trailing plasma, a practical simplification can be applied. Here, the magnetic field vector resides 
only in the two dimensional (z,9) plane, therefore a, 7 = 0, and the transformation depends only 
on one angle, (3. 
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(3.8) 



Applying this transformation to both j and E we have the following relations 



Ei = E r ji = j r 



E 2 = cos (3 Eg + sin (3 E z j 2 = cos (3 jg + sin (3 j z 



E 3 = - sin )3 Eg + cos (3 E z j 3 = - sin (3 jg + cos (3 j z 



(3.9) 



(3.10) 



(3.11) 



Overall, these relations lead to two distinct formulas for r]j_, and one formula for r/|| . Let us 
examine the simplest of these first, equation (3.9). Recalling (3.6) immediately gives us an expression 
for the transverse resistivity. 



E r E r + v g B z - V z Bg 
V± = — = : 

Jr Jr 



(3.12) 



In this expression, the different terms are known with different levels of precision. The 
axial velocity of the plasma v z is very precisely measured using a robust time-of-flight analysis. The 
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values of the magnetic components B z , and Bg are also well calibrated, and our analysis has shown 
that we can make firm upper bounds on the error of these diagnostics. 

Next in line is the radial component of the current density j r , which can be estimated using 
Ampere's Law, and our Lagrangian reconstruction of Bg , according to 

Jr = ^— 3.13) 

oz 

The reconstructed axial variation of Bg(z 7 1) is most accurate in the neighborhood of the probe, and 
so this is where we will evaluate the entire expression for r]±. 

As a check of our math if nothing else, we can also establish some general upper and lower 
bounds on the magnitude of the radial current density by dividing the measured total current flowing 
down the center conductor by an estimated cross section area equal to the surface area of the con- 
ductor where the circuit is completed by the radial conduction through the plasma, j r ~ I/A est i m ate- 
In order to be in qualitative agreement with the basic nature of the CTIX discharge, the overall 
axial extent of this region of return current should be in the range 5 cm < Az < 50 cm, and so 
^estimate = 27IT2 Az should be between about 478 cm 2 and 2394 cm 2 . Using a typical value of peak 
accelerator current of 120 kAmps, we find that we can conservatively limit the radial current density 
to be in the range: 50 Amps /cm 2 < j r < 500 Amps /cm 2 . If our reconstructed value of j r fell 
outside of these limits, then we would know there was a serious error. 

Alternatively, if we use a more educated estimate of the area that supports radial current, 
based on the axial extent of the region where Bg(z,t) has a slope in the z-direction, then we could 
make an independent evaluation of the average j r using the measured accelerator current. The 
results of these two methods can be compared, and their difference would be a good indication of 
the uncertainties inherent in the measurement. 

We have quite a different situation when it comes to evaluating the radial electric field E r . 
It is less constrained by existing measurements. We can make a simple estimate based on the voltage 
difference between the inner and outer electrodes (at radii n and r 2 ), in a vacuum the radial electric 
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field is 




rln(r 2 /ri) 



The first problem we encounter is that the voltage V acc is a non-trivial function of position and time 
during the CTIX discharge. In the future it would be best to use a good circuit model to estimate 
V acc (z, t), but for now we are limited to making a back of the envelope estimation, and then evaluate 
the vacuum expression for E r at the outer radius. 

Voltage can be estimated from Vacc(t), but it may be better to assume that E goes to zero 
after going through the plasma sheath at wall. Sheath I ~ Ad which is smaller than effective depth 
of probe 

The toroidal velocity vg is also largely unmeasured. Because our preliminary attempts to 
measure the toroidal ion velocity using the high resolution Doppler spectrometer yielded null results 
(see chapter 6), we strongly believe that vg be small compared with v z , yet it may not be exactly 
zero. A reasonable estimate is that 0 < vg < lcm/ [is. 

One final constraint is that any physical value of resistivity must be a positive number. 
This matters because each of the terms in formula 3.12, could in principle be positive or negative. 
So if j r is negative at some point, then (E r + vgB z — v z Bg) must also be negative at that point, so 
that the ratio is positive. 

We can now put all of these factors together, each with an upper and lower value that has 
been derived from actual CTIX data. 



dificulty with differnce of two similar large numbers, with errors larger than the difference 
can use the Spitzer model to 



As mentioned earlier, there is a second way to find rj± using equation (3.10). Solving for 
r]± yields the result 



where the components of the effective electric field are Eg = Eg + v z B r and E z = E z — vgB r . 




(3.14) 
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However, we are only examining the regions where B r ~ 0 and so these simplify Eg = Eg, E z = E z . 



Also according to the definition of (3, tan/3 = -jf- 



. It follows that 



*l± = 



E e (l 



(3.15) 




From here we can use Faraday's law to find Eg in terms of the time derivative of B z . 
This derivation here is motivated by the measurements of the transverse resistivity within collisional 
plasma undergoing magnetic reconnection on the MRX device [refer to MRX paper] , which found a 
very good agreement with their measurements and the classic Spitzer value of transverse resistivity. 
We can equate the closed line integral around a toroidal loop of the electric field with minus the 
rate of change of magnetic flux through that loop. 



If we rely on the assumption of approximate axisymmetry then the line integral is exactly propor- 
tional to Eg, and Faraday's law reduces to 



where n is the radius of the inner conductor and r the radial location where Eg is being 



measured. Use Bessel function force-free radial dependence (or leave it arbitrary and see if it is 
guaranteed to divide out) 

[figure, showing geometry of toroidal path integral in space] 

discuss why equivalent method for E z would be difficult (instead consider 4.1?) 

Other model on CT back edge where Bz = 0 write down matrix, solve, Any better? 

The positive result on the MRX device may be due to their simpler magnetic geometry. 




(3.16) 




(3.17) 
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Chapter 4 

Reconnection and compression 
experiments 

Reconnection within magnetized plasma occurs when two separate regions of magnetic field 
come in contact in such a way that the magnetic fields of the two regions are oppositely directed in 
the vicinity of the interface surface (separatrix). A change of the topology of the magnetic field then 
occurs in which the separate field lines merge at the reconnection zone, and become a single larger 
field line that now connects the two regions. 1 Reconnection of magnetic fields in a conducting fluid is 
only possible if the field is able to diffuse relative to the fluid. However, due to the high conductivity 
of laboratory and astrophysical plasmas, classical resistive diffusion rates are far too slow to account 
for the incredibly rapid rates of field line merging observed in actual reconnection processes such as 
the tearing mode instability within tokamaks. The challenge is to explain these anomalously fast 
diffusion rates occurring within a wide class of laboratory and astrophysical plasmas. 

Reconnection has been identified as an important mechanism in a variety of processes 
throughout the full extent of the CTIX accelerator system. It plays a critical role in the dynamo 

1 However, the notion that field lines arc indistinguishable from one another leads to a subtle fact, that the time- 
evolution of a single field line is not uniquely defined. See appendix D. 
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formation process of the CT, it limits the acceleration timescale, and is important when considering 
injection/refueling into a tokamak type reactor. 

During the formation of a compact toroid, magnetic reconnection occurring on a small scale 
(possibly via microturbulence) allows internal flux cancellation and enables the plasma to relax into 
a minimum energy configuration. 

Next, it is a larger scale reconnection process in the radial direction that allows the CT to 
pinch off at its back edge, forming a closed magnetic topology, just like blowing a soap bubble. 

Then, during the acceleration the system consisting of the pushing field, the CT, and the 
radial current sheet at the junction between them can be modeled as a region where right-angle 
reconnection occurs. Across this separatrix, the field difference is actively maintained by the radial 
current driven between the electrodes by the accelerator circuit. However, current sheet does broaden 
during the shot as it resistively diffuses away from it central plane. 

Reconnection across this interface may also play a role in the blow-by acceleration instabil- 
ity, [ref Kevin's] When the pushing magnetic field significantly exceeds the strength of the internal 
magnetic field of the CT, the current sheet junction between the two fields can destabilize near 
the inner conductor, where the the pushing field is the largest and the CT field is weakest. When 
this happens, the pushing field can "blow-by" the CT, and the result is poor accelerator efficiency. 
Onset of this instability may be due to reconnective flux cancellation of the larger Bg eating up the 
smaller Bqt- This can occur at a rate on the order of the Alfven transit time. Further analysis 
is needed to determine if this mechanism could maintain a positive growth rate for long enough to 
significantly alter the CT magnetic geometry. The blow-by effect can also be explained in terms of 
a Rayleigh- Taylor instability between the relatively heavy CT, and the lighter pushing field. Re- 
connection could add to this effect, and complicate the magnetic structure that would occur during 
this kind of blow- by event. 

Lastly, reconnection provides a very efficient method for the deposition of a fuel plasma 
into the center of a reactor field. This is the primary engineering application of CT1X, and there 
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is still much to be understood about this process. As discussed in chapter 5, we have found that 
during this kind of reconnection the plasma will thermalize, and the electron energy distribution 
function will become modified away from a simple Maxwellian equilibrium. 




Figure 4.1: 3D Rendering of axial field coils on CTIX . 



This chapter presents the results of a series of experiments that studied the interaction 
dynamics that occur when an accelerated compact toroid is injected into a solenoidal target magnetic 
field, which can be aligned either parallel or anti-parallel to the internal field of CT at its leading 
edge. Depending on the alignment direction, basic resistive MHD predicts that the interaction 
should result, respectively, in either compression or reconnection of the two fields. The results of 
the three independent diagnostics provide consistent and complementary new information about the 
MHD interaction physics within the CT plasma. 
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4.1 Experimental overview 



Four external magnetic field coils have been installed on the CTIX accelerator for the 
purpose of applying magnetic perturbations to CT within the accelerator region. The coils are 
solenoids that are concentric with the center axis of the accelerator and produce an axisymmetric 
B field that is directed along the axis of CTIX near the middle of the coils, but also has a radial 
component near the ends of the coils. Because the magnetic field produced by the coils is primarily 
directed axially, we refer to these coils as "axial field coils" . 
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Axial Field 



Figure 4.2: Parallel axial field results in a compression effect. 




The direction of the external magnetic field produced by these coils can be aligned either 
parallel or anti-parallel to the internal field of CT at its leading edge. Depending on the alignment 
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direction, basic resistive MHD predicts that the interaction should result in reconnection of the two 
fields for the case of anti-parallel alignment, or a compression effect in the case of parallel alignment. 

These effects have been measured on CTIX using the three magnetic field probes located 
along the accelerator, a laser interferometer chord to measure plasma density at the 142 cm position 
on the accelerator, and also by using a fast framing camera that records a sequence of 12 images of the 
plasma along an axial line of sight. Using these images, the spatial and temporal dependence of the 
visible light emission of the CT during the interaction is compared to its emission intensity without 
the external field. The results from the framing camera data are consistent with the independent 
measurements of plasma density and magnetic field. Overall, we observe that magnetic field strength, 
plasma density, and acceleration kinematics of the CT are modified by when we apply an external 
magnetic field, by an amount that depends on the strength and polarity of the external field. Figures 
4.2 and 4.3 show side views of the external axial coils and fields in the two geometries. 

4.2 Description of external coils 

The four axial field coils are located at z = 65 cm, z = 100 cm, z = 120 cm, and z = 148 
cm, as measured from the start of the formation section to the start of each coil. The coils are 
composed of 25 turns of 2 cm diameter welding cable. Each coil is 15.24 cm long and has an inner 
diameter of 20.32 cm, and produces 1.85 Gauss per Amp at coil center. 

Insulating acrylic support spools encircle the 16 cm outer diameter accelerator section of 
CTIX. The spools are suspended by an independent support structure that carries the weight of the 
coils and also maintains a 2 cm gap between the inner diameter of the acrylic spool and the OD of 
the stainless steel accelerator. This gap is needed because the accelerator section is charged to high 
voltage (up to 15 kV) and can heat up to as high as 140° F as a result of repeatedly carrying high 
pulsed current loads (120 kA for ~ 20 /is) during shot runs that can last many hours. 

The welding cables have been wound around the spools to form the coils, and connected 
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in series to a variable DC power supply. Each coil can be individually connected or disconnected 
from the circuit to produce a variety of field geometries. Only 3 of the coils were used for the 
reconnection/compression experiments, those being the coils at z = 100 cm, 120 cm and 148 cm. 
The coil at z = 65 cm was disconnected because early tests indicated that it disrupted the CT 
formation process when run at high current. 

Typically we cycled the power applied to the axial coils so that they would not not overheat. 
They would be ramped up to full current over the course of one second, reaching a flat-top about 2 
seconds before the shot, then ramped down starting about one second after the shot. In total there 
was a 3 second window for the shot to occur in, which was always long enough. For 20-second shot 
cycles, the on-time was about 1/5 th of the cycle period. 

4.3 Magnetic fields during reconnection and compression 

The external fields generated by the axial field coils are relatively small compared to the 
peak field of the CT, yet they have a noticeable effect on the plasma. The peak field of the CT is 
a high as 12 kGauss, while the strongest external fields we applied were only 400 Gauss near the 
center, a ratio of about 30:1. 

First, let us consider the effect that the external field has on the strength of the internal 
field of the compact toroid plasma as the two fields are driven into each other. The results gathered 
from the data set consisting of (x?) shots are displayed in figure 4.4. 

The applied axial magnetic field (Baf) was scanned from - 400 Gauss to + 400 Gauss, 
and we plotted the resulting average value of the peak B z field at each setting along with error 
bars representing the variance of the data at each setting. Reconnection occurs for negative values 
of Baf, while magnetic compression occurs for positive values. The data shown in figure 4.4 does 
not include any of the proposed calibration factors discusses in chapter 2, and so is probably an 
underestimate of the true values. 
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Figure 4.4: Effect of external field on peak B z of CT 




Figure 4.5: Effect of external field on B z at z = 142 cm 
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When the axial field is aligned anti-parallel to the internal field of the CT (when Baf 
is negative) we observe a significant decrease in the peak value of the B z field as it propagates 
into the axial field. This is consistent with magnetic flux cancellation that should occur during the 
process of magnetic reconnection. The amount of reconnective flux cancellation that is observed in 
our experiment is directly proportional to the strength of the applied external field, given the same 
initial CT field strength. 

Conversely, when the axial field is aligned parallel to the internal field of the CT, the peak 
value of B z is observed to increase as the CT enters axial field. This is consistent with magnetic 
flux amplification due to a compression of the magnetic field in the axial direction. 
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Figure 4.6: Effect of external field on CT velocity 



The magnetic probe signals can also be used to measure the time-of-flight kinematics of 
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the CT. In figure 4.6, we plot the average velocity of the CT as a function of applied axial magnetic 
field. Here the average velocity "at 91 cm" can be more precisely defined as the average velocity of 
the CT from 57 cm to 91 cm, v 91 = (91 - 57)/(i 9 i - i 57 ), while v 142 = (142 - 9l)/(t 142 - t 91 ). 

In the case of the reconnection field geometry (anti-parallel Baf < 0) we observe an 
increase in the acceleration of the CT that is directly related to the applied external field. This 
means that as reconnection occurs the magnetic energy of the CT is converted into kinetic energy 
of the plasma's forward motion. 

When compression occurs, it does so at the expense of the forward acceleration. Dynam- 
ically, an opposing external force is required to compress the plasma and embedded fields. The 
time-of-flight analysis shows dramatic decrease in rate of acceleration, down to zero acceleration 
when B AF = +400G. 

Some insight regarding the interaction dynamics can be gained by examining the axial de- 
pendence of the magnetic field under these different circumstances (using Lagrangian interpolation). 

The closed magnetic structure of the CT results in a localized pulse in the B z signals that 
propagates down the accelerator. Under normal operation (Baf — 0), probe signals show B z fields 
that remain approximately constant in size and shape as the CT is accelerated. When the axial field 
is applied, noticeable changes in CT shape occur. 

The toroidal component Bg is affected in an unusual way by the reconnection process. 
During reconnection, the BT component of the CT advances into the external field, ahead of the 
region where Bz is localized, which is normally considered to be the front edge of the CT. This effect 
can be understood if the 3D nature of the reconnection is taken into account. Since it is the poloidal 
components of the CT and external field that are undergoing reconnection, the toroidal component 
of the CT field (and the associated poloidal current density) is unconstrained in the reconnection 
zone and can freely expand into the region ahead of the CT. bringing some diffuse plasma with it. 

When thinking about the interaction dynamics it is important to remember that the ex- 
ternal field lines pass through the walls of the vessel, and they are securely pinned to the wall as the 
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Figure 4.7: Effect of external field on CT at t = 3.4 (is 




Figure 4.8: Effect of external field on CT at t — 7.2 (is 




Figure 4.9: Effect of external field on CT at t = 8.4 (is 
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CT passes because of the wall's high conductivity. As mentioned in chapter 2, the time varying part 
of B must not have any component that is normal to the wall. However there can be a DC magnetic 
field (constant in time) that intersects the wall at any angle, such as the external field created by 
the axial field coils in this experiment. As the internal field of the CT interacts with the external 
field, this DC component which intersects the wall must remain, while new time-varying fields are 
created in superposition with the original. 

4.4 Interferometric measurement of density increase 

For this experiment we used the He Ne laser interferometer to measure the time evolution 
of line averaged electron density at the z= f 42 cm port position on the accelerator. The processed 
interferometer data is shown in figure 4.10. The signals from each data set were averaged together, 
and the result was smoothed with a moving average filter to eliminate high frequency noise. 

Effect of Axial Magnetic Field 
on Plasma Density at z = 142 cm 
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Figure 4.10: Effect of external field on electron density of CT, as measured by the interferometer 
chord at z = 142 cm 

There are several effect that are evident in this data. Most clearly, we have a direct obser- 
vation of the magnetic compression effect in which the line averaged electron density as measured 
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at z= 142 cm position is greatly increased (by a factor of 1.6 for Baf — 400 Gauss), and the arrival 
time of the peak density is delayed. These are both in agreement with the idea that there is a force 
acting in the -z direction, due to the applied external field. We also see that this density build up 
scales in a nearly linear fashion with the strength of the applied external field, when running in the 
polarity that causes compression. 

We also see that there is almost no effect on the value of peak density when running with 
the external field polarity that causes reconnection. It remains the same as with no axial field 
applied, however, the signal length does seem to be stretched out for the case of Baf = -400 Gauss. 
In this case, it seems that the reconnection effect was strong enough to keep some of the plasma 
back in the accelerator. While at the lower setting, the density profile was mostly unaffected. 

One last observation can be made. It seems that with all the external field configurations, 
the leading edge of plasma density comes a few microseconds later with axial field on, than with no 
field at all. It is unclear exactly why this would be so. This leading edge has lower density, and 
weaker internal fields, and would be more strongly effected by the external fields. It seems that in 
both the reconnection and compression geometries, the interaction at the leading edge causes either 
a delay or a diminishmcnt of plasma density at the 142 cm position. 

It is possible that these are actually opposite dynamical effects that result in the same 
final measurement. In the reconnection case, a very high rate of acceleration might occur for the 
low mass leading edge, resulting in an early diminishment of plasma density because it has greatly 
expanded while being swept far ahead of the CT center. On the other hand, the compressional case 
would result in simply a delay of the arrival of the density, and when it does come it has been built 
up to a higher peak value. 
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4.5 Optical measurement of reconnection and compression 



The fast framing camera diagnostic is particularly useful because it can give us information 
about the radial distribution of plasma and how it changes with time. The Imacon 790 camera we 
used can take a sequence of 12 images of the CT during the interaction. It has an adjustable frame- 
rate that can varied in fixed increments from one frame every 2 /is to frame every 0.25/is. Using 
these images, the light emission intensity of the CT during the interaction can be compared to its 
emission intensity without the external field. For the reconnection experiments the framing camera 
is particularly useful because it simultaneously provides a time history of the axial position of the 
CT, as well as the radial and angular variation in light emission. 



Perspective drawing of main visible features 
of CTIX vessel within image frame 
(Drawn at an exaggerated angle) 




t = 5.8 us 



I I I I I I I I I 

10 20 30 40 50 60 70 80 90 



pixel number 



Figure 4.11: Diagram of CTIX vessel structure seen in framing camera image, and grey-scale image 
of framing camera output at t = 5.8 /zs. 



For this study we used an exposure of 400 ns, with an inter-frame time of 2 [is. Visible 
light emission is primarily from impurities in the trailing plasma, and so the framing camera is able 
to image the plasma sheath in the transition region between the CT field and the purely toroidal 
field behind the CT that is caused by the railgun acceleration current. 

Neutral atoms cannot be co-moving with the CT because they are not directly acted on by 



72 

the accelerating fields. However, they may be scraped off the wall as the CT goes by. Once ionized, 
the story is different, and impurity ions can follow behind the CT for a period of time. 

The framing camera is well suited for measuring the level of angular variations in plasma 
density and current density in the CT, as an indicator of the quality of the equilibrium. The line of 
sight bounces off of an internal front surface mirror, then through a glass vacuum window, and then 
via a sequence of external alignment mirrors, and into the f 4 zoom lens mounted onto the framing 
unit. The output phosphor plate of the Imacon is viewed by a 384 x 576 pixel CCD camera that 
makes a digital image of the shot. 

Some care needs to be given to properly consider the line of sight used by the framing 
camera, and what parts of the vessel are within the frame of the image. The view is looking down 
the barrel of the accelerator, but is slightly off-axis (fig. 4.11). We can clearly see the circular 
end of the center conductor, situtated at the center of the image and taking up about 1 /3 rd of 
the total diameter. Because the line of sight is slightly tilted, we can also see some of the side of 
center conductor shaft. The larger diameter of the formation section center electrode is also visible, 
especially at early times when the formation region is brightest. Some of the outer wall can be seen 
on one side. Most of the structure we see in the image is due to the plasma hugging the wall, and 
being affected by the edges of the vessel. The plasma near the wall has a higher impurity content 
and so it emits more visible light than the cleaner, fully ionized core. 

One feature (which has advantages and disadvantages) is that all the light is integrated 
axially along each line of sight. On one hand, this is good because it means that we can see all the 
way back to inside of the formation section. On the other hand, it is bad because it can be hard 
(or sometimes impossible) to determine what axial position the light is coming from, whether it is 
axially localized and moving toward the camera, or spread out uniformly along the entire line of 
sight. 

The radial variation in vessel geometry does provide some help with this issue. The clearest 
example of this is the diagnostic opportunity provided by the right angle cylindrical edge of the center 
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conductor end. At early times in the discharge the surface at the electrode's end appears dark, simply 
because no light-emitting plasma has made it down the accelerator that far. When plasma begins 
to come off the end the accelerator, it will normally fill in the region in front of the end of the center 
conductor, and begin to light it up. Thus we can be certain that when we see light from in front of 
the end-disk of the center conductor, that plasma has passed that point, about z = 165 cm. 

A similar right-angle radial step occurs at the end of the formation section, which also 
allows us to see when the plasma is entirely in the formation region, and then watch as it flows out 
to fill in the region in front of the step. 

The overall result of inspecting such features is that the images taken during the reconnec- 
tion experiment demonstrate a reduced amount of light from in front of center conductor. Whereas 
during the compression experiment, we saw an increase in light emission from back in the accelerator 
region. 
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Figure 4.12: Effect of reconnection with axial field on light emission from CT 
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Figure 4.13: Compression effect of axial field on light emission from CT 

We have compared the average of 20 images taken with no axial magnetic field to the 
average of 20 images taken when the CT is injected into an axial field of - 400 Gauss (in the anti- 
parallel direction) . We have also compared the average image taken with no axial magnetic field to 
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the average image taken when the CT is injected into an axial field of + 400 Gauss (in the parallel 
direction). In the time-sequence of images below, the red pixels show where the plasma was brighter 
with the axial magnetic field turned off, and the blue pixels show where the plasma was brighter with 
the axial magnetic field turned on. For black pixels there was no difference in plasma brightness. 

When the CT is fired into an anti-parallel magnetic field, we observe a decrease in light 
emission from the region in front of the center electrode, implying plasma is prevented from filling 
in this central region. After reconnection, the expected mechanism for plasma fill-in of the region in 
front of center conductor is cross-field diffusion. Less plasma is seen in this region for Baf — — 400G 
than for Baf = 0. These framing camera images have demonstrated the efficient trapping of CT 
plasma on the external field lines, and have provided information about the radial dependence of 
reconnection effects. 

When the CT is fired into a parallel magnetic field, we first observe an increase in light 
emission from the region between the center electrode and the outer electrode, implying that the 
plasma density and corresponding luminosity greatly increase before the CT is able to leave the 
accelerator section. For Baf = +400G it appears that most of CT plasma does not leave the 
end of the center conductor. During the compressional interaction its brightness increases due to 
the increase in density. And consequently, there is a delay while the internal field increase due to 
axial compression and the trailing plasma is prevented from getting past the compression interaction 
region. Thus it also make sense at later times during the discharge that there is less center electrode 
fill-in when compression is occurring. 

It is very satisfying that the conclusions based on the framing camera data are consistent 
with the data taken on the same set of shots from surface magnetic probes that and the interfero- 
metric measurements of plasma density. 
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4.6 Model of two interacting magnetic fields 

When the Compact Toroid enters an axially directed magnetic field, there is either an 
attractive force or a repulsive force of interaction between the two fields. In certain simple ways, 
this is just like two magnets being brought together along the axis of their poles. Whether the 
interaction is attractive or repulsive depends on the polarity of the external axial field relative to 
the internal field of the CT plasma. 

For these experiments, the magnetic polarity of the CT is fixed, as it is determined by the 
electrical polarity of the CT formation circuit and formation seed field, while the externally imposed 
magnetic field can be easily varied in polarity and strength by adjusting the DC power supply, with 
no consequence on the quality of the compact toroid formation. 

To understand this qualitative effect of either attraction or repulsion, it is simplest to 
consider the poloidal component of the internal field (the component which is created by toroidal 
plasma currents), since its geometry is analogous to the field of a simple permanent magnet. In 
standard operation the CT is accelerated with the North end of its poloidal field at the leading edge, 
so that the poloidal magnetic field vector (near the center axis) points in the direction of motion, 
which we call the positive z direction. 

Then if we were to apply, via a solenoid, an external axial field with its North pole pointing 
backward toward the oncoming CT, we would expect the two fields to repel each other (North repels 
North). With the addition of this repulsion force, the acceleration of the CT would be decreased as 
it enters the external solenoid. 

On the other hand, if we had oriented the external field so that its South pole was facing 
the CT's North pole as it approached, then the CT will be attracted toward the external solenoid 
(South attracts North). In this configuration the linear acceleration of the CT will be enhanced 
above the normal value caused by the rail gun circuit alone. 

While this simple intuitive model of magnets attracting or repelling one another is correct 
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on a basic level, and does predict the approximate outcome of actual experiments, it leaves out many 
important and complicating subtleties. 

The first detail that must be considered is that a compact toroid is a plasma, not a rigid 
body, ft is a conducting fluid with a nearly frozen-in magnetic field, and is capable of compression 
and expansion, as well as vortex-like and complicated shear flow patterns. 

When the CT is accelerated into a repulsive target field we would expect that the CT 
will be compressed in the z direction since it is being pushed on from both ends, from behind by 
the pushing field of the accelerator and in front by the repulsive axial field. The CT density and 
internal field will increase until some sort of pressure balance is reached. Although this compression 
is relatively straightforward, there are some interesting effects that are possible. Radial variation of 
the external field can provide a weak spot for the CT to squeeze through. If the external field is 
high enough reflection of the CT may occur. 

A very different process occurs when the CT is propelled into an attractive field. In that 
case, the full complexity of the plasma shows itself during the process of magnetic reconnection. This 
is a topological change in the structure of the magnetic field in which the effective rate of diffusion 
of magnetic field lines relative to fluid elements typically exceed what would be expected based on 
the bulk resistivity of the plasma. 

The two separate magnetic fields merge as the outermost field lines of the compact toroid 
reconnect with adjacent field lines of the external field; each pair of closed lines becomes a single 
longer field line that encircles the still distinct inner field lines that have yet to reconnect. This is 
accompanied by (or mediated by) a layer of strong electric current at the surface of reconnection. 
The key condition for reconnection is that the magnetic field changes sign in at least one component 
as the surface of interaction is crossed. 

Overall, the most interesting complications come about when considering the geometry of 
the internal and external magnetic fields and the boundary conditions that the conducting walls of 
the vessel apply to the plasma and fields. 
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Figure 4.14: Sketch of flux surfaces during reconnection at three times A. t = 3 /is, B. t = 6 /is, 
C. t = 9 /is. Red field lines are members of the external field from the axial coils, before the have 
reconnected with the CT field, which is depicted in blue. The green curve in each figure shows the 
location of the reconnective current sheet at the separatrix between the opposing fields. The current 
in this sheet would be flowing in the toroidal direction, coming out of the page. Behind the CT the 
toroidal pushing field is shown as a set of x marks indicating a magnetic vector going into the page. 
The thick orange field line is the last field line that intersects the both the inner and outer wall. 
The smaller field lines that are nested inside the orange curve will only intersect the outer wall. The 
large turquoise ellipse represents the approximate location of the original plasma of the CT as time 
passes. 



Shown in figure 4.14 is a sketch that illustrates the qualitative behavior and geometry of 
reconnection of the CT with the axial field. Certain features are disproportionate, such as the ratio 
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of internal CT magnetic flux to external flux, and the drawing is not to scale. It approximately 
maintains a divergence free B field while satisfying the nontrivial boundary condition that force the 
external field lines to be pinned to the wall on the experimental timescales. This leads to several 
effects, the most obvious is that the field lines that do intersect the wall will at certain times become 
compressed against the wall, greatly increasing the magnetic flux density (B) until it matches the 
strength of the larger internal fields. 

Given an internal field that is larger than the external field, there will be a final flux 
surface of the CT that will undergo reconnection with the external field. Reconnection will procede 
until... After the last flux surface has reconnected, the core of the CT will essentially "slip past" 
any remaining opposing external flux until it comes out the other end of the coils, where is will meet 
with a parallel field that interacts compressively. This is depicted in frame C. of figure 4.14. 

To preserve the divergence- free nature of B (equivalently, by Gauss's Law) the contact 
angle between the field lines and the wall on the post-reconnection side of the CT, (the back side) 
must obey 

.Baf (wall) , „ „. 

tanem ^ cr (wall) (41) 

which will lead to very small values of Contact due to the small ratio between Baf and Bqt- This 

is constrained by the fact that Baf (wall) must maintain the same value as it had before the plasma 

discharge began. On the other hand, -Bct (wall) will be evolving in time as the CT passes by. 

This simple geometric model correctly predicts the interesting feature of negative B z mea- 
sured by the probes after the CT has passed by (see fig. 4.5), which is the result of the radial 
compression of the remnants of the anti-parallel external axial field, as shown in fig. 4.14 C. in the 
region between z = 57 cm and z = 91 cm. 

The best conceptual picture for the reconnective increase in the acceleration of the CT is 
that ultimately, the presence of an anti-parallel field is simply a catalyst to unlock the large amount 
of stored internal magnetic energy of the CT. The gain in kinetic energy of the CT comes mostly 
from the magnetic energy stored in the outer layers of its volume that are shed as it reconnects. This 



79 

fraction of its total volume can be estimated very roughly using a few approximations. In the future 
this could be turned into a more precise calculation. The idea is that the magnetic energy lost due 
to reconnection should be equal to the gain in kinetic energy of the CT. If wc use the approximation 
that without reconnection, the total magnetic energy of the CT is roughly equal to the total kinetic 
energy, then we can express the energy balance in terms of a propotional change in energy. Stated 
in this way, we see that the mass of the CT divides out, 

AU m (reconnect) AKE (reconnect) vfj — Vq 

U m (total) = ATS (total) = v^ ^ ' ' 

We can simplify the magnetic energy ratio if we assume a nearly uniform distribution of energy 
density, 

A U m (reconnect) = L ihc f dV = AV shc ii . . 

U m (total) ~ J B 2 dV Vtotai ( ' ' 

V total 

where AV s h e ii is the volume of the outer region of the CT that attaches to the axial field during 
reconnection and V to tai is the total volume of the CT. Then the fraction of the CT volume that 
gets sloughed off can be calculated from the measured final velocity in the case of reconnection 
(vfl <~ 19 cm/fis), and that without (vo ~ 18 cm/ us). 

^-"^ 

In other words, approximately the outer 10% of the CT will be separated from the core during 
reconnection with an external magnetic field of 400 Gauss. Alternatively, we can say that 90% of 
the CT makes it all the way through the external field intact. 

Given that this magnetic energy is released in a finite period of time, we can estimate the 
average rate of flux annihilation. The total process is at most about 5 fis in duration and the total 
magnetic energy of the CT is about 500 J, which means that the outer 50 J of magnetic energy was 
liberated from the CT, resulting in about 10 J / [is or 10 megaWatts of reconnection power. 
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Chapter 5 

Temperature measurements 

5.1 Particle velocity and energy 

In this chapter we will discuss measurements of the electron energy distribution function of 
the plasma. By taking the first moment of this energy distribution we arrive at the average thermal 
energy of the electron population, which is proportional to the electron temperature of the plasma 
for a Maxwellian equilibrium, (Ethermai) = \nKT e . Here n is the number of degrees of freedom 
that the plasma electrons can move in, which is 3 for our system, and K = 1.38 x 10 _23 J/°K is 
Boltzmann's constant, and so KT has units of energy. Because the electron and ion populations are 
distinct interpenetrating fluids, it is possible for them to not be in thermal equilibrium with each 
other and exist with different energy distribution functions, and different temperatures. 

When thinking about the kinetic energy of the plasma particles it is important to distinguish 
between their drift and thermal velocity components. The drift velocity of the plasma is defined 
to be the average velocity of all the particles with a macroscopic fluid element, relative to the rest 
frame of the laboratory. In other words, the drift velocity is the fluid flow velocity of the plasma. 
It may vary with position, since the fluid elements may be flowing along the streamlines of some 
non-uniform velocity field. 
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In the CTIX accelerator, the plasma drift velocity is dominantly in the direction of the axial 
acceleration, although there may be some smaller toroidal or poloidal component of the velocity of 
individual fluid elements. On the largest scale, however, if we average over all the particles in a thick 
annular slice at a given axial position, any poloidal rotation of the fluid will average to zero over 
the cross section, while toroidal rotations may exist but are expected to be much smaller than the 
dominant axial velocity. This is what we often call the CT velocity, and it is what we can measure 
with time-of-flight analysis of sequential magnetic probe signals, or with Doppler spectroscopy under 
certain conditions. 

On the other hand, the thermal velocity component of any particle is simply the velocity 
vector that remains when we subtract the drift velocity from the particle's velocity relative to the lab 
frame. In this way, a particle's velocity vector is the vector sum of its thermal velocity component, 
and the drift velocity of the fluid. 

Vparticle — ^ drift "T ^thermal 

This thermal velocity will vary randomly in direction and magnitude as we go from particle to 
particle. We can characterize this random variation by the distribution function of the thermal 
velocity. We will also separately take account of the drift energy, and the thermal energy of each 
particle, 

Ednft = -^nv 2 dnft , E therma i = -mvt hermal 
where m is the mass of the particle. 

5.2 Thermalization of electrons 

A thermalization effect occurs when a compact toroid is injected into the vacuum magnetic 
field of the Davis Diverted Torus. Experiments were conducted in which the CT was injected 
transversely (_L) into magnetic fields of 105, 150, and 180 Gauss, each defining a distinct data set. 
Due simply to the net motion of the CT the ions have a drift kinetic energy in the range of 100 eV 
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to 200 eV. The electrons, in contrast, gain little kinetic energy from the net motion of the plasma, 
owing to their much smaller mass (1/1856 of the mass of a proton); the net drift contributes no 
more than 0.1 eV to the electrons, which is insignificant compared to the electron thermal energy of 
10 eV or greater. 

During the process of collision and reconnection with the magnetic field of DDT the initial 
velocity field cannot stay directed axially. As the center of mass of the CT is brought to rest, the 
original directed kinetic energy is converted into thermal energy as the anisotropy of the particle 
velocity field is lost to a more isotropic velocity field in which entropy is maximized. While this 
basic concept provides a general description of thermalization, it leaves many questions unanswered. 
A more complete and predictive model would need to delineate the mechanisms that enable the 
conversion of particle energy to field energy and vice versa, as well as coupling between particles, 
such as the energy transfer from the ions to the electrons. These mechanisms could include kinetic 
effects, possibly via collisions, wave growth instabilities, and kinetic damping. Fluid effects could 
also be important, such as a turbulent cascade of fluid energy down to microscopic dissipative scales, 
shock heating, and magnetic reconnection. The theory behind these issues is far too extensive to 
tackle within this dissertation. However, the basics can be outlined, and several new experimental 
findings can be reported, which may help future work toward confirming existing theories, or guiding 
the invention of better theories. 

The fundamental concept needed to describe the process of thermalization is the distribu- 
tion function of the thermal energy of the particles; we will notate the energy distribution by Je(E). 
The energy distribution function Je{E) expresses the electron population density as a function of 
energy. If it is normalized to unity (the convention that I usually prefer), then Je(E) also represents 
the probability that any single electron will be found to have a thermal energy in the neighborhood 
of the energy E. 

And of course there are other relevant particle distribution functions. For example /e{E) 
is related to the speed distribution by Je(E) = Jf v (v(E)) where v is the speed, or magnitude of 
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the velocity of a particle, and J is the Jacobian of the transformation v — > E(v) = \mv 2 . Then 
J = \dv(E)/dE\ = |(2niE)~ 1/2 |- 

Maxwellian equilibrium For reference purposes I would like to include some relevant formulas 
here. The velocity distribution function for a Maxwellian equilibrium of temperature T, with n 
spatial degrees of freedom, normalized to unity is 

where the n components of v range over the domain from — oo to oo. Although the velocity distri- 
bution depends only on the speed v = |v|, it is critical to remember that the distribution function 
for the speed v of a particle (i.e., the speed distribution) is not the same as the velocity distribution. 
The Maxwellian speed distribution function in the case of n = 3 is 

where the domain of f v (v) is 0 < v < oo. From the speed distribution we can write down the energy 
distribution function for a Maxwellian system (n = 3), 

E 

f E (E) = ^-^(KTyim 2 VE e "kt (5.3) 

where the energy is, of course, always positive; 0 < E < oo. 

On the CTIX device, we have successfully measured electron energy using a pair of electro- 
static energy analyzers. We found that immediately before the collision/reconnection the electron 
thermal energy was approximately 30 eV to 40 eV and afterward the thermal energy increases to 
greater than 70 eV. Ion thermal energy is also expected to increase during the collision, however 
experimental verification of this effect remains incomplete. An attempt was made to repeat the 
energy measurement on the ion population using a larger 4-grid Ion Analyzer Probe. This was 
ultimately not successful. The measurement of ion thermal energy is more challenging due to the 
larger diameter ion Larmor orbits and certain plasma sheath effects. 
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A note on distribution functions A quick comment should be included here regarding a differ- 
ence of terminology. The function that physicists call a distribution function (as I have been using 
it above), is the same thing that mathematicians somewhat more precisely refer to as a probability 
density function (or pdf). The math nomenclature does a better job of conveying the meaning of 
the function because the integral of a pdf over a given interval results in the probability that the 
continuous random variable will take on a value within the limits of integration. In this way the 
value of the pdf at any given point represents the density of probability of random occurrence near 
that point. It reminds us that we need to integrate the pdf to get a value that is actually a prob- 
ability. To add to this confusion, mathematicians have given the name 'distribution function' yet 
another definition, F(X) which is the integral from — oo to X of the pdf f(x). The mathematician's 
distribution function measures the probability of the event that the random variable will take on 
a value less than X. This will never be what I mean when I refer to the distribution function of 
some variable. I will always use 'distribution function' as the physicists do, equating it with the 
mathematician's pdf. And since my audience is primarily composed of physicists, I will follow their 
convention and mainly speak in terms of distribution functions. 

5.3 Experimental overview 

In our experiments we measured the electron energy distribution before and after it interacts 
with the DDT target magnetic field using electrostatic energy analyzer probes (EAP) in a variety of 
configurations. For all of these experiments we measured the initial velocity of the CT as it entered 
the injection region, as well as the final velocity of the reconnected plasma as it propagated along 
the field lines of DDT. 

We wanted to know how the initial and final states scale with injection energy, and so we 
repeated the same measurements for three different accelerator settings, CTIX accelerator voltage 
at 9.5 kV, 12 kV, and 15 kV. For this sequence of data sets, we adjusted the strength of the tokamak 
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field to follow the increased accelerator voltage, (-Bddt = 105, 150, and 180 G). We did this with 
the intention of maintaining the CT stopping condition. As we will see, because of some accelerator 
inefficiency, the CT injection velocity did not depend on V acc according to a simple, monotonically 
increasing function, as we naively expected. Instead, the velocity peaked near 12 kV, and then fell 
off as the applied accelerator voltage was further increased. The result is that we explored a much 
broader variety of stopping conditions than we originally planned, [tabulate stopping parameters 

for the three settings, to show that 9.5 kV data set should not have been stopped by the field] 
The details of the design and operation of gridded electrostatic energy analyzer probes is presented 

below, which includes an analysis of simulated data to arrive at an estimate of the experimental 
error, as well as some statistical results regarding the systematic error due to the analysis algorithm. 

5.4 EAP construction 

The probe tip is a hollow cylindrical assembly, 2 cm in diameter, made from anodized 
aluminum with a moderately large 1 cm diameter aperture to let plasma inside. Within the probe 
tip are the following components: two circular grids, (the ion repeller grid and control grid), and a 
circular collector plate. They are stacked together in a sandwich configuration with mica spacers to 
insulate them from each other and from the case of the probe tip. The probe tip is supported by a 
28" long stainless steel shaft which slides through an o-ring vacuum seal to allow the probe to be 
inserted into the vacuum vessel over a range of depths. 

The two grids are tungsten mesh of high transparency, and the collector is tantalum. The 
two grids and the collector each have a thin nickel wire spot-welded to them that runs back through 
the shaft and connects to the electronics that operate the probe. By way of these wires the grids can 
be charged to different potentials and a current can be drawn from the collector plate, and recorded 
by a digitizer. We will ultimately derive the electron energy distribution function from measurements 
of how the electron current drawn by the collector plate depends on the applied repeller grid voltage. 
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Figure 5.1: Exploded view of EAP showing grids and collector plate 



For our experiments the energy analyzer probe tip was positioned at the center of the 
vacuum vessel with the aperture oriented toward the incoming plasma to have the best opportunity 
to intercept a good signal. For the injection experiments involving collision/reconncction on DDT 
target field, the EAP was mounted on the top port 90° from the injection point on DDT; then for 
the pre-injection measurements, the probe was positioned in the drift section. 



5.5 Grid biasing 

The functionality of the energy analyzer probe is achieved primarily by charging the control 
grid to a negative potential with respect to ground. This potential on a probe is often referred to 
as the bias voltage or simply the bias. Negatively charged electrons that initially enter the aperture 
of the probe will feel a repulsive force from this grid. The electric field of the control grid can be 
thought of as an electrostatic hill standing between the electrons in plasma and the collector plate. 
This hill will act to separate those electrons which are fast enough to make it over the hill from 
those which are too slow to get over the hill. All of the electrons that do make it over this hill are 
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then promptly absorbed into the metal collector plate, causing an AC-coupled negative current to 
flow across a bias capacitor and then back to the digitizer. 
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Figure 5.2: Schematic of EAP biasing circuit 

The magnitude of this current is directly proportional the total number of electrons in the 
vicinity of the probe which had enough energy to make it over the hill. In Fig. 5.3 we have shown 
the collector current as a function of time over the course of a typical shot. 

1 OOO Shot Averaged EAP Signal 




Figure 5.3: Time dependence EAP collector current, parameterized by control bias voltage. As the 
negative bias is increased, a larger electron fraction is repelled. The first strong negative peak in the 
signal corresponds to the arrival time of the compact toroid. 



Now, as we vary the bias voltage applied to the control grid, we will repel a variable amount 
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of electrons in a manner that is energy dependent. At a larger negative bias, more electrons are 
repelled by the control grid, letting only the very fastest electrons get to the collector plate, hence a 
small current is drawn. When we decrease the magnitude of the bias, very few electrons are repelled 
hence a very large current. At zero bias, all the electrons can get by and so the current has leveled 
off to a maximum value. 

The other grid in the probe, the ion repeller grid, is situated to be the first to encounter 
incoming plasma and is charged positive in order to prevent ions from entering the interior of 
the probe. In the simplest possible biasing circuit, the positive bias on the ion repeller grid is 
accomplished using a high impedance battery that provides a DC voltage yet blocks the fast signals 
caused by interaction with the transient plasma. 



Everything that I have already described about the operational mechanism of the EAP, and 
how it can be used to measure the electron temperature of a plasma can be succinctly summarized in 
the following integral relationship between I(V) and Je{E). Because the electron energy distribution 
/e(E) measures the probability density of electrons in each infinitesimal slice of energy space, we 
can express the collector current I(V) as an integral of the electron energy distribution. 



Here, q is the charge of an electron, it works as the conversion factor between V in Volts and E 
in electron- Volts. This integral counts all of the electrons which have an energy above (qV), which 
are collected as electron current by the EAP when the bias voltage is V. Randomness enters our 
measurements primarily from the electron flux rate T = n e v e A co u ector which varies from shot to shot 
with an approximately Gauss-normal distribution. Here n e is the electron number density, v e is the 
average electron speed (dominated by thermal motion), and A co u ector is the area of the collector 
plate. The physical source of this variation is mostly due to irreproducibility in the total gas output 



5.6 Relation between I(V) and Je(E) 




qV 



(5.4) 
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of the formation gas valve, and also partly due to spatial variations in plasma density and how that 
density intercepts the probe. Now, it is the current that is experimentally measured, initially as the 
set of data points making up the I-V scatter plot, one (I, V) point per CTIX shot. We first need to 
extract an approximation of the function I{V) from the raw data by some averaging or curve fitting 
procedure. Once we have the measured /(V), we can solve equation 5.4 for the unknown Je(E) by 
differentiating both sides. 

ME) = -L"£> (5.5) 



qT dV 

q 

We also see that V — E/q is the bias voltage which is able to repel an electron of energy E. 

Lastly, the average thermal energy is computed by taking the first moment of the energy 
distribution function. 

/•OO 

(Ethermal) = / Ef E (E)dE 

Jo 

So the experimental task is to find a feasible way to measure the I-V curve with the EAP so that 
we can differentiate it to get the local electron energy distribution function. In our setup we do not 
try to measure the energy distribution of the plasma for each individual shot because the desired 
plasma conditions in the vicinity of EAP are so transient. The plasma wave passes over the probe 
for only a few microseconds and then is gone. While it is physically possible to sweep the bias over a 
few hundred volts at a MHz rate, we believe that the cost and inherent difficulties in implementing 
and operating such a system would overshadow any benefits for measuring the fast plasmas on our 
experiment. 

It is more experimentally feasible to build up the I-V curve point by point, one point for 
each shot. To make this pointalistic composite of the I-V curve we collected 20 shots per bias 
voltage setting, with a total of 20 settings from 0 volts to -190 volts in increments of -10 volts. This 
gave a total of 400 shots per data set in the original post-reconnection measurement of electron 
temperature. 

We were careful not to take all 20 shots for a single bias voltage at once. We broke up the 
bias settings into even and odd multiples of -10 V. We would first go sequentially through the even 
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Figure 5.4: Analysis method using moving- window smoothing 

set, with 10 shots at each setting, then do a round on the odd set with 10 shots per bias setting. 
Then we repeat this to get the second 200 out of 400. 

The reason for dividing the scan into a several repetitive cycles was to separate out any 
systematic trends in system performance from the real electron temperature effect. For example, it 
might be a problem if all the low bias voltage shots were taken at the beginning of the day, and all the 
high voltage shots at the end of the day. There would no way of knowing that as the device warmed 
up over the course of the day, variation in the gain of the electronics or in the performance of the 
accelerator was the real source of the observed trends in the data. Based on the comparison of these 
divided data sets, we found that there was no measurable systematic variation in the performance 
of the EAP diagnostic over the course of several weeks of runtime. 
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5.7 Analysis methods 

Method 1: Two Parameter Nonlinear Fit to a Maxwellian form I(V) Curve. This method 
assumes that plasma is in thermal equilibrium with an effective number of degrees of freedom that 
must be chosen by the user. The fitting function returns a temperature and a normalization factor. 
The idea is that if the electrons had a Maxwellian distribution, then the I(V) curve for the EAP 
should obey some analytically expressible formula. If we evaluate equation 5.4using the Mawellian 
energy distribution given in 5.3, then we find 



Method 2: Non-analytical Free Fitting Method. Using no analytical assumptions about 
the form of the Energy Distribution Function, this method is able to reconstruct deviations from 
thermal equilibrium that may exist within the plasma. The Raw Data is averaged, smoothed and 
then differentiated numerically to yield a discretized plot of the Energy Distribution Function. Both 
Methods have been implemented and compared. 

5.8 Experimental results 

[figure, graphs (f(E), and velocity histograms, for each DS)] with a few comments, and a 

table 

Each of the 10 colored line-plots in the 1000 shot Average is the point-by-point time- 
dependent average of the 100 shots that were all taken at the same Bias Voltage. On the right 
are individual EAP signals for two shots that contributed to the 1000-shot average graph. The 
important features here are first, the strong negative peak that occurs at the time of CT arrival, 
and second, the weaker negative pulse that is due to the trailing plasma eventually getting to the 




(5.6) 
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Figure 5.5: Raw data {•} with average- value curve {violet}, iteratively smoothed curve {black}, and 
Maxwellian fitted I-V curve {blue}. The x-axis measures the absolute value of the EAP bias voltage, 
and the y-axis measures the resulting EAP collector current (arbitrary scale). 



Table 5.1: Summary of results from 9.5 kV data set 



V acc = 9.5 kV, B DDT = 105 G 


Pre-injection 


Post-reconnection 


Fluid velocity 
Net Ion Kinetic Energy 
Electron Thermal Energy 


13.65 cm/ (is 
97.2 eV 
unknown 


5.95 cm/ (is 
18.5 eV 
70.1 eV 



Table 5.2: Summary of results from 12 kV data set 



V acc = 12 kV, B DDT = 150 G 


Pre-injection 


Post-reconnection 


Fluid velocity 
Net Ion Kinetic Energy 
Electron Thermal Energy 


18.4 cm/ ' (is 
176.7 eV 
unknown 


14.95 cm/ (is 
116.6 eV 
86.4 eV 



probe. The top signal was taken with electron repeller bias voltage set to zero and you can see that 
both the CT and the trailing plasma cause electron current to be collected by the probe. When we 
increase bias voltage to the -90 V in this lower shot, it is clear that almost all electrons have been 
repelled. 

Pre-injection electron measurements Method 1: Maxwellian Fit Average Thermal Energy 
= 33.47 eV Errors are unknown 
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Table 5.3: Summary of results from 15 kV data set 



V acc = 15 kV, B DDT = 180 G 


Pre-injection 


Post-reconnection 


Fluid velocity 
Net Ion Kinetic Energy 
Electron Thermal Energy 


12.14 cm/ \is 
76.9 eV 
unknown 


10.26 cm/ lis 
54.9 eV 
80.3 eV 



Method 2: Free Fit Average Thermal Energy — 43 eV The Error due to Data Set Size and 
Quality is cr\ hermal Energy = ±2.02788 eV The Error due to Smoothing is s2Smoothing = ± 2 eV 

5.9 Error analysis of simulated data 

We shall closely examine two cases, one simulation done with a Maxwellian input dis- 
tribution and one done with non-Maxwellian input distribution. Here is the Maxwellian energy 
distribution, [figure ?] It has an average thermal energy of 1.4 (arbitrary units). We integrate this 
analytically to create the ideal I-V curve, from which we can generate a simulated data set which 
constitutes the "measured" I-V curve. This is shown below. 

This data has 400 points in it, just like a real data set. To make the points scatter, for 
each bias setting, we generate 20 normal-random numbers from 0 to 2 which have the same standard 
deviation as the real data set does and we use these numbers as the vertical distance of the data 
points away from the ideal I-V curve. To do this we just multiplied the 20 normal-random numbers 
by the height of the ideal I-V curve to get the scatter of points in right relative frame above and 
below the ideal curve. 

If the input and output distributions were wildly different in the simulation, then we would 
expect that the experimentally measured distribution would also be wildly different from the ideal- 
ized, physically real distribution. If the input and output distributions were very similar, and this 
similarity held for a wide range of input curves then we would be guaranteed an equivalent level 
of similarity between our measured energy distribution and the actual energy distribution of the 
plasma. We then accept the simulated data as being "real" and do our analysis on it. 
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Below is the energy distribution reconstructed from the data. It has close to the same 
average thermal energy of 1.44 which differs from the average value of the original by about 3%. 
The distribution's shape has spread out quite a bit due to the effects of randomization in creating 
the simulated data. This is not an error in the average-smooth-differentiate algorithm because when 
simulations were done with 2000 data points in the set, we reconstructed the proper shape of the 
distribution quite well. 

As the number of data points shrinks, the reconstructed distribution becomes correspond- 
ingly more round in shape. This spreading is due primarily to our data size of 400 points. In 
addition, different randomizations produced slightly different data sets, and hence slightly different 
reconstructed energy distributions. This kind of fluctuation was examined informally, and it was 
observed to be relatively small in the Maxwellian case with 400 points. No simulated sets produced 
humps where they were none to begin with. 

In this way all of the Maxwellian reconstructions were easily identified as being Maxwellian. 
For comparison, we tried using a non-Maxwellian distribution as the input for the simulation. Here 
we have the original distribution, its average thermal energy is 3.3. Using the same simulation algo- 
rithm as before, we create the "measured" I-V data set as shown below. We then try to reconstruct 
from this data set and yield the reconstructed energy distribution. While its shape is quite different 
from the original, its average thermal energy is still very close, 3.46 which is within 4.8% of the 
original. Many other functional forms were tried as the input for the simulation. Similar results 
were obtained for all simulations using the same statistics as our actual data sets. In all but a few 
of these runs, the average value of thermal energy was correctly reconstructed to within 5% of the 
original. 
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5.10 Statistical error analysis 

The main result from the simulation work was that the average value of thermal energy 
had a relatively small total error for the size of data sets we are working with. In this section we 
take this result a step further by deriving an exact relationship between the data set dimensions and 
the resulting error in the average thermal energy of the plasma. 

We consider a situation where the probe current I(V) has been measured ny times at the 
k bias voltage settings, where in general, ny could take different values at different bias voltages, 
hence the V index is used to denote this possible dependence. A non-constant ny could easily occur 
in a real life experimental data set, where system failures of one sort or another could require that 
data from certain shots be excluded from the dataset. The result would be a slightly patchy scan of 
the parameter space,where more good data was gathered for some voltage settings than for others. 



To get the Energy Distribution the data is averaged, smoothed until it is a well-behaved 
monotonic function, and finally differentiated to yield the Energy Distribution. 



For now, we will ignore the smoothing step, which only acts to decrease the error by 
bringing the data closer to the real distribution. Without the smoothing step we will be looking at 
upper bounds for the errors involved with this sampling process. The variable T was observed to be 
normally distributed about the mean density. 

Because each I(V) measurement depends on the random variable T and a monotonic func- 
tion of V, the variance of the I data at V obeys 



JqV 

Note that the variance at V = 0 should be equivalent to the variance that describes how 
the electron flux rate T is randomly distributed (i.e., ct 2 (0) = Op). 



[figure I-V raw data, caption: ny samples of probe current k bias voltage settings] 



[refer to figure 5.4 for analysis process] 




(5.7) 
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As mentioned earlier, the random variation in I(V) is due almost entirely to fluctuations 
in T, while the stable dependence on voltage is due entirely to the integral of Je(E), and so the 
variance will decrease with increasing \V\. 

This condition is useful as an indicator of whether or not the probe is really working. If 
<J 2 (V) doesn't decrease then there is a problem with the probe. Error bars for the values of the 
/(V) curve at each voltage can then be estimated from the measured <J 2 {V) and the number of data 
points n v at that energy. For each voltage setting, the variance of the average value of the (n v ) 
probe current measurements is related to the variance of the measurements themselves, according 
to 

«4.(V) - ^ (5-8) 
fly 

Of more importance to this study is an upper bound on the error of the average thermal 
energy, quantified by the variance first moment of Je (E) . We used moment generating functions to 
make the transformation from the variances of the individual measurements to the variance of the 
error in the average thermal energy. The resulting total variance is just the sum of the variances of 
the average value at each voltage V\ , V2 , ■ . . , Vk- 

_ a^V,) a\V 2 ) a\Vk) 

thermal error _ _ _ W'^/ 

n Vl n V2 n Vk 

The key feature in this expression is that the error will tend to grow linearly with k (where k is the 
number of bias voltage settings in the scan, also the number of terms in this sum) . It must be noted 
that it does not work to try to minimize the error by taking k =1, since your choice of that single 
voltage setting will dominate the final answer. Enough voltage settings must be included in order to 
have a minimum amount of energy resolution to consider the final result a reasonable approximation 
to the distribution function. With that in mind, k should be at least 10, and then n must be chosen 
to be large enough to bring down the error, (k =10, n=200 are good values to plan an experiment 
around). 
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5.11 Summary of results 

The simulation has shown the limitation of accuracy of our experimental method when 
only 400 data points are taken. The exact shape of the distribution cannot be trusted, except that 
if it has some lumps, it is very unlikely to be from a Maxwellian plasma. If it has no lumps, the 
physically real energy distribution still has a chance of being non-Maxwellian. On the other hand 
we can put more confidence in the average thermal energy of the distribution. 

[table of EAP results] 

The simulation has shown that our analysis on 400 data points is fairly robust with respect 
to this average value, and we can expect that with our real data we have correctly reconstructed the 
average thermal energy of the plasma to within 5%. 

When we compare the results of the EAP study to the results of Time of Flight study we 
find that the final electron thermal energy is comparable to final ion kinetic energy. This implies 
that electrons and ions have gone through some process of equilibration during the collision. 

We notice that of the three data sets gathered in this study, the 15 kV and 12 kV sets 
have similar, non-Maxwellian energy distributions as well as post-collision loss of net flow velocity. 
The 9.5 kV set had a much greater loss of flow velocity and in addition appeared to have a more 
Maxwellian shaped distribution. This supports the model that in the 9.5 situation the interaction 
of the CT with the DDT field results in the plasma coming to a dead stop at the point of injection 
and then slowly diffusing out into DDT. 

Whereas with the 12 kV and 15 kV sets, a definite post collision plasma wave propagates 
through the vacuum field of DDT carrying significant mass and kinetic energy. The primary variable 
which must be used to account for this difference are gun velocity and DDT field strength. For 12 
kV and 15 kV the field strength was 150 Gauss and 180 Gauss, respectively. For the 9.5 kV set 
the field was only 105 Gauss, yet the gun velocity was higher than for the 15 data set. This could 
indicate that in the 9.5 kV case the field was not sufficiently strong to capture the CT before it 
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collided with the back wall of the DDT vacuum vessel. If CT-wall collision occurred then the post 
collision plasma would likely have completely different thermal properties, 
comments about apparent non-Maxwellian nature of results. 



99 

Chapter 6 

Spectroscopic measurements of 
ion velocity 

6.1 Previous low resolution spectroscopy 

One of our standard instruments has been general purpose visible survey spectrometer. It 
can measure the entire visible spectrum 300 nm to 670 nm on a single exposure of a gated/intensified 
linear detector array (1024 pixels). This allows diagnosis of the light emitted by impurity species 
within the plasma. Two modes of operation are possible. We can integrate over the entire shot, or 
gate the spectrometer to measure light from a small time interval within the shot. Time evolution 
of spectrum during the acceleration process can be built up over a sequence of shots in the gated 
mode by stepping the delay time of the start of the spectrometer gate pulse relative to the time 
of CT plasma formation. The spectrometer can view the plasma axially via an internal mirror, 
or transversely through glass window ports. Focusing optics are required for time-gated operation 
in order to collect enough light into the fiber optic bundle to make detectable light signal at the 
diffraction plane. 
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Figure 6.1: Time sequence of low resolution spectra. 



6.2 High resolution visible spectrometer 

A transmission grating spectrometer appropriate for visible light (3500 A to 6800 A) was 
built by Alex Graf 1 in collaboration with EBIT group and was used to study the plasma created by 
the CTIX machine. The fast time resolution and high spectral resolution are possible due to special 
LLNL fabricated transmission gratings, which the design is based around. 

The transmission grating is circular with a diameter 150 mm, groove spacing of 3944.5 
+/- 0.5 A. Due to its large diameter this grating allows faster integration times, it has an f/3 
which is much greater than the f/10 in a conventional spectrometer. It has a demonstrated spectral 
resolution, A/AA, of up to 14500 for a wavelength of 5790.7 A while using a 50 [ira slit (Fig. 1 and 
2). The optics within the spectrometer allow for off axis (± 1 cm) input light, [ref RSI paper] 

My part in the spectrometer project was primarily the development of practical implc- 

^^Graduate student in UCD Physics Dept. 
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Optical parameters 

Slit width = 50 una 
Resolution AX = 0.3A to 0.8A 
Range of A. : 3 650 A. to 6750A 
Disp ersion = 11.4504 A/mm 
or 0. 1 4648 A per pixel 

Transmission Grating: 
2535 lines per mm 
150 mm diameter 

Lenses 
Focal length = 50 cm 
f = 3 . 8, 1 30 mm diameter 

CCD 

Pixel size 12.8 x 12.8 um 
Gating time: 0.5 us to 5 us 



Input fiber optics 



Center Wavelength 
K = d(sin a + sin (3 ) 
d = 3945 A 




Transmission 
Grating 



Figure 6.2: Operational schematic of high resolution spectrometer 




1 . Transmission grating 

2. Achromat lenses 
■3. Detector 
■4. Input fiber optic 



Figure 6.3: Rendered 3-D cut-away view of high resolution spectrometer 
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mentations of fast detector arrangements suitable for viewing the CTIX plasma. We constructed a 
sequence of custom-made adjustable supports for mounting any available fast CCD camera that we 
had at our disposal into the interior of the spectrometer so that it could record the spectrum image 
at the dispersion plane. All together, we took data using a sequence of 5 distinct detectors in as 
many months of run time. 

The first experimental task that we completed was a time integrated impurity survey of the 
light from a helium CT plasma. The spectral resolution is far superior to that of the former device 
allowing a proper separation and identification of the many features present. This high- resolution 
survey spectrum spanned most of the visible wavelengths from 3500 A to 6800 A. This spectrum was 
compiled from about 50 individual spectra, and it mapped out the location and relative intensities 
of hundreds of impurity lines, as well as the primary helium and hydrogen lines. Next we selected 
one of the brighter spectral lines for further study. The best candidate was the line transition at 
4686 A of singly charged helium (He II). 

The largest shift expected for a 200 km/s CT is around 4 A well within the spectral 
bandwidth of typically 150 A, while the smallest shift we can detect is at least ~ 0.4 A (less if the 
light is sufficiently bright). 

We measured Doppler shifts of this line using an axial line of sight, where the plasma 
is rapidly moving toward the point of observation. We compared this to spectra that were taken 
along three different transverse chords that are orthogonal to the direction of plasma motion and so 
experience no Doppler shift. The spectra along all four of these lines of sight were recorded using 
a time gating method that mapped out the time evolution of the spectral lines and their Doppler 
shifts in 1 microsecond increments from the moment of plasma breakdown until 30 microseconds 
after breakdown. During this period the CT plasma is accelerated from 0 to 200 km/s. 

The resulting time evolution of Doppler measured plasma velocity was a bit of a surprise 
at first. Using time gating of the CCD we are able to watch the evolution of ion velocity using the 
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Figure 6.4: Survey spectrum of helium plasma near He II line 



measured Doppler shift at each time step (see fig 6.6). He II ion velocity peaks very early at 1.5 lis, 
with a maximum of ~50 km/s and then dwindles to only 10 km/s during most of the acceleration 
of the compact toroid. A second peak of ^72 km/s occurs at 14 /is, which is coincident with the 
moment when the CT leave the end of the accelerator, at this time a re-strike occurs at z = (0 to 
20 cm) due to the ringing of the railgun circuit. 

We never saw the plasma reach anything close to the peak velocity of 200 km/s as measured 
by the time of flight between magnetic probes, and we would like to know why. 

After completing a subsequent experimental and analytical investigation, we found a satis- 
factory explanation for this effect, which relies on new evidence of a dramatic variation of ion density 
as a function of position and time that can be accurately modeled within a simple framework. 
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Figure 6.5: Doppier shift of helium II line at A = 4686 A 



6.3 Measurement of average Doppier velocity 



We are looking at a comparison between the transverse spectra and the axial spectra of 
the line emission of He II (Z=l) at A = 4686 A. 

The transverse spectra have no detectable shift, while the axial spectra show a Doppier 
blue-shift because the plasma is moving rapidly toward the input fiber. The apparent wavelength 
of the blue-shifted light is: 

A = Ao(l-«/c) 

where v/c is the ratio of the ion velocity to the speed of light, which is always less than 10 -4 for the 
CTIX plasma. Experimentally however, we cannot measure the Doppier shift of individual ions. In 
our setup, the light from nearly all the ions in the entire volume of plasma is collected by the axial 
line-of-sight fiber optic. The measured shift in center wavelength is actually a weighted average of 
the Doppier shift throughout the vessel. And so we can only infer an average fluid velocity based 
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on this average Dopplcr shifted wavelength. 



<■„„(/) = < [1-^) , J ((i.l) 



A 0 

where X ave (t) is the experimental average of Doppler shifted center wavelength as measured 
by the axial line of sight and A 0 is the unshifted center wavelength of the line emission as determined 
by the transverse measurements. X ave (t) is found by taking the first moment of the measured line 
profile. 

Ao+A 
Ao-A 

where the normalization of a single line of the spectrum at each instant in time is given by 

Ao+A 

iVj(i) = J I(X,t)dX 

Ao-A 

Here, I(X,t) is the output signal of the spectrometer when viewing the axial line of sight, 
calibrated for wavelength but with intensity in arbitrary units. The range of the integration, ±A 
away from the center wavelength, is just large enough to include all of the line profile, without 
including any neighboring lines. If we assume a Gaussian line shape for emission from an infinitesimal 
fluid element we can approximate the spectrum as the integral over the axial coordinate z over the 
length of the accelerator L 

L 

I(X,t) = -^= I I(z,t)e-^- x ^ 1 '^^ t ^ c '>'> 2 ^ 2 dz (6.2) 

vv n J 

o 

where I(z,t) is the intensity of emission at the center wavelength from the fluid element 
at position z and time t, and a is the line width. Here v(z, t) denotes the actual fluid velocity of 
the plasma, which varies in time and position throughout the accelerator. The velocity function of 
the plasma is not known for all (z,t) but it is strongly constrained by measurements of the average 
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velocity determined by the time of flight between signals of magnetic probes at three axial positions 
along the acceleration section. 

Also, the results of 2-D MHD simulation of the accelerator dynamics of CTIX provide 
accurate approximations to v(z, t). However, one of the goals of this spectroscopic work is to find an 
independent measure of v(z, t) that can be compared to other results as a check of accuracy. Using 
the above expression for I(X,t) we can evaluate v ave (t). After some simplification we find 



The emission intensity I(z, t) is of key importance. We have endeavored to measure it as 
directly as possible with the use of narrow bandpass filtered photomultiplier tube measurements 
made simultaneously at several positions along the length of the accelerator. These measurements 
can be analyzed using the Lagrangian interpolation technique, in synthesis with the magnetic field 
data. Additionally, our understanding of the basic ionization and excitation mechanisms can be 
used to make a simple predictive model of the accelerated helium ions in our system. Ultimately 
we would like to compare the average velocity evolution measured via Doppler spectroscopy to the 
results from the two other independent methods, the PMT-Magnetic composite of I(z, t) the yields 
v a ve(t)i and a purely theoretic approach. 

In the light of this possible cross-comparison of methods, the geometric complication of 
the axial integrated data is ultimately a good thing since it means that the Doppler measurements 
provide some constraints on a model of ion density field rii(z,t), which up to this point had been 
uncharted on the CTIX device. 

6.4 Ion velocity evolution 

Using time gating of the CCD we are able to watch the evolution of the average ion velocity 
using the measured Doppler shift at each time 1 [is step. 



L 




(6.3) 
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Figure 6.6: Average Doppler velocity as a function of time. 



This analysis shows that He II ion velocity peaks very early at 1.5 (is, with a maximum 
of ~ 50 km/s and then dwindles to only 10 km/s during most of the acceleration of the compact 
toroid. This unusual result requires some further explanation. Careful examination of the transverse 
data show that the He II state becomes extinct as it propagates down the accelerator. Very little 
is observed at the z = 91 cm, and absolutely none at all is seen at z = 142 cm. A second peak 
of ~ 72 km/s occurs at 14 (is, which is coincident with the moment that the CT becomes fully 
detached from the end of the accelerator, at this time a re-strike occurs in the formation section 
(0 < z < 20cm) due to the ringing of the rail-gun circuit. Error bars in figure 6.6 represent the 
standard deviation of the measured average velocity at each time step. There were 10[?] data points 
per time step. 
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It is clear that the maximum velocity measured by the Doppler spectrometer is much 
less than the CT velocity of 200 km/s as measured by the magnetic probes. The reason for this 
discrepancy is all in equation 6.1. We have come to understand that because the electron temperature 
is sufficiently hot, the population of He II does not last more than a few microseconds during the 
discharge before it becomes fully ionized. Also, He II is limited in space to the region close to the 
formation section, because by the time the He II ions get to the end of the accelerator they have 
had enough time to turn entirely into a particles (see ionization diagram 6.12 on page 114). 

And so the Doppler spectrometer using the axial line of sight to view the plasma will never 
see He II ions moving at the maximum plasma velocity It can only observe them early on in the 
acceleration process, when velocities are small. The helium ions do in fact get up to the high speeds 
of 200 km/s or higher, we just can't see them once they loose their last electron. This is a very 
reasonable interpretation of the data, however it is important to carefully compare the timescales 
involved in the ionization and acceleration processes and make sure the numbers work out. In the 
next section we investigate several models and compare them to experimental results. 

6.5 Fast visible emission measurements 

When a charge state is being depleted during the acceleration process, there will be some 
diminishment of its light emission with increasing axial position. Fast silicon photodiodes with low- 
noise RF amplifiers were used to make the first preliminary measurements of this effect, by looking 
at the visible light emissions from the CTIX plasma. The light output of the three accelerator ports, 
as well as the axial line of sight were monitored in this way. 

Optical band-pass filters were used to isolate the light from distinct impurity species and 
find its evolution as a function of axial position and time. The most general observation we can 
take away from the photodiode measurements is that across the spectrum, all visible emission was 
decreasing exponentially with increasing axial position. We believe this is primarily due to thermal 



109 



electron impact ionization. 

I was involved in the task of using the results from the high resolution survey spectrum 
to select a set of optical bandpass filters that allowed light from a limited region (10 nm wide) of 
the plasma spectrum to be measured by amplified silicon photodiodes. A summer student, Michael 
Rogers 2 constructed three amplified photodiode assemblies with moderate width (10 nm) optical 
bandpass filters and measured the time evolution of light at the three diagnostic windows along the 
accelerator section. The visible emission was measured at wavelengths of 465 nm (includes He II, 
Oil, Nil), 435 nm (includes Oil) and 655nm (includes Ha line). For comparison, the total visible 
light emission was measured without the use of filters. 

Using time-of-flight analysis on these measurements, we could find a difference in the rate 
of acceleration between the trajectory of peak magnetic field of the CT, and trajectory of the peak 
light signal, for each spectral filter. We could also compare the kinematic results between the 
measurements using different filters, in order to determine if some impurities were lagging behind 
the others. This measurement served as a simple test of the hypothesized centrifuge effect where 
ions of different mass become separated in space due to a difference in acceleration. And most 
importantly, the measurement of the decay of ion luminosity as a function of axial position gave us 
the ability to track the time evolution of atomic processes such as ionization, and accumulation of 
neutrals due to scrape-off from the wall. 

As the CT moves down the accelerator it emits less and less light. Fast photodiode mea- 
surements show an exponential decrease of: Total visible light Emission near 465 nm (includes He 
II, Oil, Nil) Emission near 435 nm(includes Oxygen II) Emission near 655nm (includes H a line) 

These measurements gave some circumstantial evidence that the ionization model I pro- 
posed was correct. However, they had the problem that the filters took in light from a variety of 
different line transitions, and so it was not possible to isolate exactly what the He II line at 468.6 
nm was doing. 

2 Undergraduate in Physics at UC Santa Cruz, visiting us as part of the National Undergraduate Fellowship for 
Plasma Research, Princeton Plasma Physics Lab. 
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Figure 6.7: Log of peak radiated power of filtered light (blue) and unfiltered light (red) as a function 
of axial position, for filters centered at 435 nm, 460 nm, 465 nm, 655 nm. 

To improve upon these encouraging results, we purchased a set of custom made optical 
band-pass filters with very narrow transmission range. These made it possible to isolate the light 
coming from a single line transition. And with the right detectors we were able to record the 
magnitude of an individual line emission continuously in time. 

Once I began working with the narrow filters, I ran into a number of technical complications 
that we had not encountered with the wide filters. Initially I had simply plugged the narrow filters 
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Figure 6.8: Fractional distribution of power between wavelength ranges as a function of axial position, 
for filters centered at 435 nm, 460 nm, 465 nm, 655 nm. 
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Figure 6.9: Transmission curve of narrow bandpass He II filter. 



into the original photodiode housing, but that suffered from the serious problem that light plasma 
could enter the filter at all angles from 0 degrees to as much as 30 degrees off normal incidence. 

The critical effect is that the center wavelength of the transmission curve becomes shifted 
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with increasing angle. So non-helium light from neighboring spectral lines could still get through 
the filter if they enter at some significant angle. To solve this, I first quantified this effect and found 
that there is a shift of the center wavelength by approximately 0.1 nm per degree away from normal. 




-9degrees8 -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7 8 9 

A 9 nltt[ ( degrees ) 



Figure 6.10: Effect of tilt-induced wavelength shift 
Include formula for center wavelength as a function of angle. 

Using the measured rate of wavelength shift, the filter must have an effective index of 
refraction of 6 [check?] . 

Based on the survey spectra there were neighboring oxygen lines about 0.5 nm on either 
side of the He II line so we would need to limit the angular acceptance to well below 5 degrees. This 
has been accomplished with the construction of 9 cm diameter collimator tubes that have a pair of 
adjustable apertures separated by 55 cm, and internal baffles that block internal reflections off of 
the tube wall. 

The narrow bandpass filters are mounted in micrometer adjustable universal-joint swivel 
mounts that allow precise alignment of the plane of the filter, in order to orient them exactly 
perpendicular to the axis of the aperture pair. The apertures can be varied in size to optimize the 
trade-off between angular acceptance and light signal level. 
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Micrometer adjustable Photomultiplier Collimator Assembly 

universal swivel mount 

for f Iter ^ var j a y e Fixed Fixed i r i s variable 

aperture aperture aperture aperture 




Incoming light 



black fabric enclosure 



Figure 6.11: Drawing of collimated PMT assembly 

Under typical plasma conditions the angular acceptance can be limited to less than half a 
degree, while still allowing good light signal. However, many orders of magnitude of the plasma light 
are thrown away when using the combination of the collimator and the narrow bandpass filter. The 
resulting light signal was far too dim to be detected by the original silicon photodiode and amplifier 
circuit. A much more appropriate choice of detector is a photomultiplier tube (PMT). They have 
nanosecond response times and can be sensitive enough to count single photons. We happened to 
have an old RCA model 8852 PMT and I was able to borrow two more of the same model from the 
SSPX group (Sustained Spheromak Physics experiment). 

[figure of time signals of PMTs at the various positions with explanations in captions] 

6.6 Ionization model 

At the instant of plasma formation there is already a large cloud of moderately high density 
neutral helium filling the formation section and diffusing lightly out into the accelerator section. Once 
a small population of energetic electrons becomes liberated within the static electric field, ionization 
of the neutral helium proceeds rapidly and the electron and ion densities rise primarily through the 
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process of single electron impact ionization. 




Figure 6.12: Helium ionization and excitation process 



Figure 6.12 diagrams the excitation and ionization process that occurs in the CTIX helium 
plasma. Note that the excitation and ionization energies are comparable. This leads to the general 
effect that if the plasma electron population is hot enough to cause visible line emission from a 
species, then it also has enough energy to ionize that species into the next higher charge state. 

As discussed in chapter 5, electron probe measurements show an electron temperature of 
T e = 54 eV within the accelerator. This is high enough to rapidly burn through neutrals and singly 
charged ions leaving only an a particle plasma that does not emit in the visible. 

In this section we will model the evolution of the charge state populations once ionization 
has begun, given known cross sections. First we will consider the case where the ionization is 
occurring within a dense cloud of neutral helium, such as in the formation section. Secondly, we 
examine the case when the neutral (He I) density is zero, but there is a population of He II, 
which subsequently evolves into a population of fully ionized helium without replenishment from 
the neutrals. This case should describe what happens at the leading edge of the CT plasma where 
a population of singly charged ions has been accelerated out of the neutral cloud, and will evolve on 
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its own from then on. 

In both of these models we begin with the following assumptions: 

• Constant electron temperature — > Constant ionization cross sections 

• Closed system — ► No particles enter or leave the system. 

• Ignore vessel geometry effects — > Zero-dimensional + time 

• Ignore recombination, double impact ionization, and photo ionization. 
Recombination rates are small (a ~ 10~ 14 cm 3 s _1 ) and so will be ignored on the fast 

timescales in the formation. In fact, ions will not recombine before the experiment is over (25/zs) 
because the lifetime of ions (before they recombine) is l/(o n e ) <~ 1 second. Because of this, the 
ionization fraction of the plasma is a monotonically increasing function of time. 

The rate of double ionization due to electron impact is more than an order of magnitude 
smaller than for single ionization, and does not have a significant effect on the lifetime of the He 
II in the CTIX plasma. First, He II can only lose one more electron before it is fully ionized, so 
double ionization does not affect He II at all once it has formed. Second, the small effect that does 
exist acts primarily to slightly reduce the maximum value of He II density that is obtained during 
the evolution by allowing a small fraction of the Helium neutrals to bypass the He II state and 
go straight to the fully ionized state. This does not significantly affect the lifetime of the He II 
population, and so double electron impact ionization will be ignored. 

[figure of cross section as a function of electron energy] 

Comparison of cross sections 

Give numerical value for cross section used in these calculations, [ref], describe method of 
integrating 

Formation Region Model In the formation region, where the neutral density is high, the evolu- 
tion of no(t) and n»(t) are coupled through two first order nonlinear ODE's and an algebraic equation 
for n e (t). The density functions can be formulated in terms of the constant N = n 0 (— oo) which is 
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the number density of neutrals before any ionization takes place. 

At time t = 0 the initial conditions are n 0 (0) — N(l — e), rij(O) = n e (0) = Ne where e is 
a number much less than one. This model is closed in that no particles enter or leave the volume. 
However, if we assume that the ionization rates So and Si (for neutrals and singly charged ions 
respectively) to be constant in time then some large amount of energy must enter the system to 
keep the electron temperature constant. 

If it is desired, the doubly charged ion density (a -particle density) can be easily accounted 
for using conservation of nucleons once no(t) and rij(i) are found. The evolution equations are: 



= -S 0 n e n 0 = -S l n e n i + S 0 n e n 0 n e (t) — 2N - 2n 0 (t) - n^t) 



We can rephrase this system in terms of dimensionless variables 

T = tS Q N = t/t Q K = Si/So a(T)=n 0 {t)/N (3{t) = m{t) / N 7 (r) = n e (t)/N 

The evolution equations for a, (3, 7 are 

da d/3 „ 

— = — a-f — = — A/57 + a l 7 = 2 — 2a — p 

dr dr 

These have initial conditions a(0) = 1 — e, (3(0) = 7(0) = f«l 

If you solve for 7 and then combine the two ODE's you can find the path in (a, (3) phase 
space by solving 

da a 

Provided that K ^ 1 then the curve /3(a) is given by 



f3 =T ±^(a K -a) 

The evolution of the system can then be fully described in terms of a single variable ip = a~ x 
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dip 1 da 7 
dr a 2 dr a 

_„__!_) (6 .4) 

We do not know if a general closed form solution exists for this equation, and so we have 
mainly worked with numerical solutions instead. However there is a special case, when K = 1/2, in 
which equation 6.4 can be transformed into a linear equation and solved exactly. 



For K= 1/2, iP(t) = (Ae T + l) 2 where A = 




This corresponds to a density evolution of 

no(t) = N{Ae s ° Nt + 1)- 2 (6.5) 
m(t) = 2NAe s « m (Ae SoNt + l)- 2 (6.6) 

This happens to be a reasonable approximation to the evolution of the helium plasma 
created in CTIX with its value of K = 0.57 based on the calculated cross sections. 
[ figure showing evolution curve in (no, n^) phase space] 

[ figure showing no, ni,na, n e as function of time for typical CTIX plasma parameters] 
discuss any key features 

Accelerator Region Model Because the neutrals are not effected by the electromagnetic fields 
that accelerate the ions, the slow moving neutrals don't have enough time to make it into the 
accelerator section during the acceleration of the CT plasma. The absence of neutral helium simplifies 
the evolution equation for the singly charged ions while they are being accelerated. 

We will now consider the evolution of the ion density with respect to a new time variable 
t' = t — t a where t a is the time the fluid element has just entered the acceleration section. The 
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initial condition is rn(t' = 0) = TVj which is the final ion density rii(t a ) using the formation region 
ionization model. 

Because some doubly charged ions were created in the formation region, the accelerator 
initial condition for electron density is n e (t' = 0) = N e = Ni + 2Nu, again using the formation 
model to find the doubly charged ion density at the time of entering the accelerator Nu = nu(t a ). 
The system evolves according to: 

dft ■ 

-^- = -Sin eni n e {t') = N e + N - m(t') 
which can be stated more compactly as 

fin ■ 

^ = -S i n i (2(N i + N ii )-n i (t')) 
Restated in terms of a few dimensionless variables, 

r' = t'SM R = Nu/Ni /3(t') = n^/N 

The evolution equation becomes a logistic equation: 

^ = (3 2 - 2(1 + R)(3 0(0) = 1 

This can be solved to give an exact solution of the ion density evolution in the accelerator 

region: 



m(t') - 2JVi(l + R)[(2R + l) e 2^(iVi+iv«)*' + i]-i 
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6.7 Reconstruction of I(z,t) from PMT data 

[figures, typical graphs vs time of raw signals of: Axial view, 57, 91, 142, ALREADY HAD 
THIS IN 6.4] 

velocity field deduced from magnetic probe signals. 
[Figure, I(z,t) shown as a function of z for various times] 

6.8 Comparison of methods 

Compare theoretical evaluation average velocity with PMT reconstruction 
[figure, overlay plot of 3 different methods] 
briefly discuss size of errors, any problems if any. 

6.9 Conclusions 

This has been a cross-confirmation of several different physical effects, and as a result we 
have established several facts about the formation and acceleration process, as well as developed 
some sophisticated analysis tools that we can use again in future work. 

As diagnostic, a Doppler spectrometer has some unique advantages, and certain limitations 
that are worthwhile to take into account when planning an experiment in this field. One key practical 
result is that the fast rate of ionization imposes a limitation to the size of the region where Doppler 
measurements can be made. 

It happens that the A = 4686 A line of He II is one of the brighter lines in the plasma. 
Its brightness is good for visible spectroscopy because it allows short integration times, however the 
high magnitude of brightness also means that He II very rapidly becomes fully ionized long before 
the experiment is over, making it impossible to observe the maximum velocity of the CT via Doppler 
shifted line emission. 
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Unfortunately this effect also puts some constraints on the usefulness of higher Z impurities 
for visible spectroscopy with our plasma. For a plasma with a moderate electron temperature, there 
will be a highest charge state that can be obtained from electron impact ionization. For example, 
consider Ar V (Z=4) which is the most strongly charged argon ion that can be created in a T e = 100 
eV plasma. If recombination is negligible, then the Ar V will persist until the plasma discharge is 
complete. This terminal state has bound electrons, unlike fully ionized helium, and so can be made 
to emit visible light. 

The difficulty is that because the excitation and ionization energies are comparable, Ar V 
only emits very dimly in the visible. If there is not enough energy to ionize, then there is not enough 
energy to excite. Attempts were made to measure visible line emission of argon ions on the CTIX 
plasma with no real success so far. The best results were with helium, and in spite of its limitations, 
the data was very revealing about the basic plasma physics. 

6.10 Future work 

Plans are underway to improve these Doppler measurements with a new spectrometer, and 
allowing better optical access using periscopic mirrors placed within the ports on the accelerator 
section. The goal is to observe the red shifted light only from the plasma that has passed the port 
and is accelerating away from it. 
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Chapter 7 

Determination of 
plasma Reynolds number, 
hydrodynamic and magnetic 

We would like to characterize each plasma discharge on CTIX according to certain di- 
mensionless parameters that describe the flow. In particular we will use previous results regarding 
resistivity to estimate the magnetic Reynolds number of the plasma. 

7.1 Definitions 

Hydrodynamic Reynolds's number 

Re = LV/v (7.1) 
= Inertial force/Viscous force (7.2) 

Where v = fluid viscosity which will be estimated from a formula because we currently have no way 
of measuring it on our system. 
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why does high Re correspond to turbulence, why not low Re, why does is matter at all. 



Magnetic Reynolds's number 



(7.3) 



(7.4) 



= Flow velocity/magnetic diffusion velocity 



(7.5) 



where r\ = resistivity, (a — conductivity), L = length scale of system, V = flow velocity of fluid 
(plasma in our case). 

R m is a dimensionless parameter that determines the degree to which the magnetic field is 
effectively frozen into the fluid j only partially true 

Use of Rm in determining reconnection rates. 
Use in Shock physics 

Different MHD systems with the same R m will have geometrically similar flow behavior 
that can be scaled to match each other under a suitable transformation of length and time units. 

How does this matter? Are we going to scale it up to other devices? or scale it down? No 
were not really. 

We can measure r\ with some confidence on our system, as described in chapter 3. Here we 
will use these measured values to estimate the magnetic Reynolds numbers that are attainable on 



CTIX. 



7.2 Range of accessible values 



brief discussion 



tabulate answers 
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7.3 Onset of turbulence 

Flows with high values of the Reynolds numbers typically exhibit turbulent behavior, while 
low values are indicative of simple laminar flows. 
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Chapter 8 

Gas puff imaging technique 

Gas puff imaging is a relatively new plasma diagnostic technique that makes it possible to 
photograph fully ionized plasmas using very short exposure times (1 microsecond or less). The idea 
is to make the plasma pass through a cloud of neutral gas, thereby causing the plasma electrons 
to collisionally excite the neutral atoms, which then emit enough light to be imaged by high-speed 
high-sensitivity cameras. 

As we discussed in chapter 6, by the end of the accelerator, the CT plasma is composed of 
fully ionized hydrogen (or helium) and emits visible light very dimly due to trace amounts of higher 
Z, partially ionized impurities. We have implemented neutral helium gas puff imaging systems, 
located on the transverse viewing chamber, as well as the z = 91 cm port of CTIX. The neutral gas 
is injected about a millisecond before the formation of the hydrogen CT plasma. 

In combination with existing diagnostics, wave phenomena and turbulent flow states are 
being observed via the density fluctuations they create. 

A target chamber with large windows (0.5 m by 0.2 m) was constructed that allows trans- 
verse imaging of a large volume of plasma. The gas puffing in the target region greatly enhances 
plasma brightness, and spatial variation of image brightness can be used to infer the spatial fluctua- 
tions of the plasma electron flux density. These fluctuations have been observed with fast intensified 
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and non-intensified cameras, with and without the use of narrow bandpass optical filters. 



8.1 Practical considerations 



timing sequence, valve and pulser design, Late puffing causes helium to be localized near 
the outlet. 

valve output vs applied voltage, density estimates 




700 710 720 730 740 750 760 770 Volts 

^bias 

Figure 8.1: Peak gas puff neutral density in viewing chamber as a function of gas puff pulser voltage 
Vbias for a sequence of test shots. 



[figure, sketch of gas valve mechanism, with circuit] 

second window on back side with black curtain as optical dump to lessen reflections 



8.2 Excitation dynamics 

discuss conclusions drawn from ionization model in chapter 6. We see that light emission 
as a function of time correlates well with electron density as measured with Langmuir probes (LMP) 
at floating potential, drawing electron current. The top graph shows the light emission measured 
along two different chords in the transverse viewing chamber. Photomultiplier tubes (PMT's) were 
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Figure 8.2: Circuit schematic for Silicon controlled rectifier(SCR) switched gas valve pulscr. 



used with a pair of adjustable apertures in order to limit the input light to a narrow chord. The first 
large peak in both graphs corresponds to the passage of the CT plasma through the drift section, 
while the second broader peak is the trailing plasma, [ref] 

quote estimate of lifetime of excited state. 

relate to PMT LMP measurements in 8.6 



8.3 Implementations on CTIX drift section 



brief discusion 



(1) Top gas puff 
Shots 54997 to 56986 



History of gas puff 
locations 



CTIX 




(2) Side gas puff 
Shots 56987 to 58349 



(3) Bottom gas puff 
Shots 58350 to 58513 



(4) Insulated 
center line gas puff 
Shots 585 14 to 60498 





Fast opening 
gas valve 



Figure 8.3: History of gas puff locations, summarized. 

Top mounted gas puff [show example, not kink] 
[Show looking up at puff] 

Bottom mounted gas puff [show example, toroidal field interaction?] 

Insulated gas puff glass tube [show example, with clear localized puff] 

Accelerator gas puffing comments acceleration of ions, transit velocity 
[figure, line drawing indicating position at 91 cm] 
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8.4 Cameras 

highspeed sensicam 
PI Max? 

Imacon (why it wasn't used) 

resolution 

ion diffusion blur 

motion blur 

8.5 Spectral measurements 

basic low res showing presence of neutral helium line, compared with and without gas puff, 
[figure, side by side spectra, no GP, with GP] 

8.6 PMT and LMP measurements 

describe as basic existence indicators, and use for timing of cameras 

[figure, overlay of PMT transverse chords and LMP (electron/ion current) to show the 
correspondence between density an luminosity] 

discuss future work of making improved absolute calibrations of LMP density, and plasma 
Luminosity to quantify their ratio, and dependence on any other parameters. 
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Chapter 9 

Wire target perturbations 

9.1 Target design 

The final set of experiments conducted on CTIX focused on an analogy with grid turbulence 
of a classical fluid in a wind tunnel. The idea is that we create a localized obstruction to the flow, 
that has a known k-value, and we watch the plasma's reaction to this perturbation as it flows 




Figure 9.1: Photograph of wire target 
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downstream. The large-window transverse viewing chamber was designed with this experiment in 
mind. 

We constructed a rotatable wire array composed of 1 mm diameter stainless steel wire 
mounted to an insulating Pyrex support tube. The wire array was designed to be rotated relative 
to the plasma flow, this was implemented without the need for a complete vessel opening by (very 
carefully) releasing most of the clamping pressure on the o-ring vacuum seal, and then rotating the 
end of the support tube with one hand while steadily increasing the o-ring clamping pressure again 
with the other hand. The support tube is hollow and allows a magnetic probe to be inserted along 
its entire depth. 

An second Pyrex tube was positioned at the tips of the wires, in order to provide a sym- 
metric configuration and also to allow magnetic measurements at a second location in the target 
region. 



Rotatable 

= (-56 to +230 degrees) 



Pyrex tubes (hollow) 

(30.5 cm X 1.27 cm diameter) 



Wires hold onto 
the support tube 
using a gentle 
spring force 



Stainless steel wires 

(10 cm X 1mm diameter) 



Magnetic field 
probes can be 
inserted over a 
range of depths 

- Probe coils at tip 



Figure 9.2: Drawing of wire target detail 



rotatablity 

magnetic probe access 
description of effect of wire array 
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Figure 9.3: Top view of drift section arrangement 



description of effect of glass support tubes, explain implantation of hydrogen, 

spectroscopic confirmation of light source [figure, (2 graphs x 2 graphs Combo EPS) images 
of before and after really bright flare up, and graphs of corresponding spectra ] 

What are the properties of the fluctuations caused by interaction with the target? How do 
these fluctuations grow or decay in the flow after the target interaction? 

pros and cons of design 
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Figure 9.4: Angular positions of wire target 




Figure 9.5: Images showing shadow effect of wire target on shots 59928 and 60139 
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Figure 9.6: Images showing residual gas on target 
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Figure 9.7: Spectra showing composition of residual gas on target 
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9.2 Observation of coherent density waves 

This graph shows the location of observed plane waves and their approximate wavelength 
A and direction. For some shots distinct waves were observed in multiple locations and so are 
listed separately. There was a noticeable amount of reproducibility of waves occurring in the same 
location, implying that the waves do not form at random, but instead are caused by some identifiable 
mechanism. 



point 


# Shot 


0 


59841 


1 


60091 


2 


60235 


3 


60276 


4 


60295 


5 


60299 
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60306 
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60438 


10 


60439 


11 


60448 
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60486 
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20 


60497 




Figure 9.8: Plane wave locations with labels according to shot number. 



Careful examination of these images gives some indication that the regions of highest 
luminosity are sometimes very close to the glass window. For example, on many of the images 
exhibiting the wire shadow effect, when the wire target was rotated to the —30° position so that the 
tips were touching the glass, the profile of the shadows closely 
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Figure 9.9: Coherent waves in plasma during shot 60438 




200 400 600 800 1000 1200 



Figure 9.10: Coherent waves in plasma during shot 60295 

Displayed in figures 9.9 - 9.12 are some examples of very localized Fourier modes, which 
we visually identify as plane waves. 

Exposure smear length is the distance the plasma should move during the camera exposure 
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Figure 9.11: Coherent waves in plasma during shot 59841 
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Figure 9.12: Coherent waves in plasma during shot 60235 



time (of 1 /is), based on time of flight measurements of plasma velocity (L exp = vqt i~ exp ). Here we 
compare this length to the wavelength of the plane waves. 

The curious fact is that if these were plane waves that were co-moving with the plasma at 
a velocity between 3 to 15 cm/ [is, then for an exposure of 1 full jis, we would not be able to resolve 
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Figure 9.13: Localized Fourier modes of plane waves, in k x , k y Fourier space (color scale represents 
magnitude of FT) 
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Figure 9.14: K-space distribution of Plane waves 

the clear wave structures on the order of 1 cm that we do see. If these were co-moving plane waves, 
the images would be mostly a blur due to the relatively long exposure time, and the fast motion of 
the plasma. We also see that there is no apparent trend that relates wavelength and plasma flow 
velocity, such as might be expected from a simple traveling wave with A/ = v^. This again implies 
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Shot number in plane wave dataset 



Figure 9.15: Ratio of expected exposure smear length to wavelength is displayed here vs shotnumber 
of data set {L exp = v C t r exp ). 



that the wave speed is not the same as the plasma flow speed. 
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Figure 9.16: Wavelength vs plasma flow speed 



The only viable conclusion is that these waves are essentially standing waves in the labora- 
tory frame of reference, with wave velocities of less than 1 cm/ /is. The mechanism that is responsible 
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for generating these density waves is not fully understood. The problem here is an over abundance 
of candidates due to the many important interaction effects that are conceivable within our system. 

One possibility is that the wire targets in the plasma did just what they were intended to 
do, and the result of the perturbation is that coherent, single-mode density waves are created if the 
plasma conditions arc right. 

Include PMT data here, examine shots where waves were seen in vicinity of chords, [figure, 
PMT time signals in drift section during shot with waves] 

These waves may be generated by the two stream instability between the moving CT 
plasma and a stationary singly ionized helium plasma. Other possible mechanisms include shock 
wave driven ion acoustic modes or that the electron density variation is due to the magnetic islands 
that form as the plasma passes through the wire array. 

Further analysis on the existing data will be required to develop a better model of what 
causes these, and other wave phenomena in the CTIX plasma. 
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9.3 Observation of turbulent waves 

Observation of turbulence during accelerator region gas puff experiments. 

One mode of operation was found to be especially useful for visualizing turbulent flow in 
the plasma. For example, the image taken on CTIX Shot 60062 (Fig 9.17) displays an apparent 
transition from laminar to turbulent flow as plasma passes wire target and support tubes, traveling 
from left to right. 




200 400 600 800 1000 1200 



Figure 9.17: Transition to turbulent flow on shot 60062 

The camera exposure for this shot was 1 fis, with a delay of 30 [is from the formation of the 
CT plasma. This image is of the relatively uniform Bg pushing field plasma behind the CT, after the 
CT has already moved through the viewing chamber. The helium gas puff was injected into the z = 
91 cm accelerator port at approximately 14 milliseconds before the shot, which primarily served to 
increase the plasma electron density by about a factor of 3 for this shot. The secondary effect of the 
gas puffing is that helium ions (which serve as the dominant source of light for fast imaging) have 
a chance to be captured and accelerated by the magnetic fields. By the time they enter the drift 
section they are co-moving with the plasma and therefore the light emission from these captured 
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ions is indicative of the structure of the flow, and evidence of density fluctuations. When gas puff 
conditions are just right, this mode of operation is ideal for flow visualization, roughly analogous to 
using smoke to trace streamlines in wind-tunnel experiments. Injection velocity of SCT for this shot 
was 3.35 cm//xs. The dominant wavelength of post-target fluctuations for this shot was 2.5 cm. 
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Figure 9.18: Fourier transforms of the image taken for shot 60062, (top) in laminar region before 
plasma has interacted with target, (bottom) turbulence after interaction. 
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Top graph in Figure 9.18 is the 2-D Fourier transform of the plasma luminosity in the 
region before it has interacted with the target. Bottom graph is the FT of the plasma after it 
has interacted with the target. The color scale represents the magnitude of the FT for any given 
wavevector (k x ,k y ), and the units are arbitrary, but the same for both top and bottom graphs. A 
significant increase in fluctuation amplitude is observed, as well as wider distribution of waves in 
the k y direction. Also the k-spectrum is noticeably more choppy and irregular for the turbulent 
flow than for the incoming laminar flow. The dominant wave modes in the turbulent spectrum are 
localized near k = (2.5, 0) [1/cm], where k = 2n/X. A custom-made discrete Fourier transform 
algorithm was used instead of a prepackaged FFT algorithm in order to provide more precise control 
of the resulting resolution in k space. 

Let us examine several other examples of apparent laminar-turbulent transition. The time 
at which the images were taken are listed in parentheses. 



Shot 60053 (t = 26 /is) 



Shot 60055 (t = 30 /is) 
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Figure 9.19: Radial profiles of luminous flow, before and after target for shot 60055. 



Shot 60061 (t= 27 fis) v\A2 = 2.6947 cm/fis Xdominant = 3.5 cm Exposure smear ratio for 
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Figure 9.20: Image of shot 60055, showing vertical lines where sampling was done. 



this image is 0.76, which allows for direct imaging of waves at Xdominant = 3.5 cm that are at most 
76% blurred due to fluid motion. If the velocity decreases in the drift section to below 2.6947 cm/ [is 
then there will be even less bluring. For this shot it is possible that the apparent wavefronts in 
luminosity are oriented transverse to the stream lines of the flow. 
Shot 60091 (27 (is) 

Discuss measurements of initial conditions of plasma, 

show possible stream lines overlayed on image, (drawn by hand) 



9.4 Comparison of experimental results with 
plasma turbulence theory 

It is also worth noting that more plane waves of the type discussed above were observed 
almost exclusively during target region gas puffing. In these shots the CT plasma collided with a 
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Figure 9.21: Wavelength vs plasma flow speed 




Figure 9.22: Possible long wavelength streamlines for shot 60061 



stationary neutral population and the result was the appearance of coherent spatial variations in the 
total plasma luminosity. On the other hand, during shots with accelerator gas puffing, the neutrals 
had time to become ionized and accelerated by the fields by the time luminous regions of plasma 
crossed through the transverse viewing chamber. Images of the plasma during these shots had a 
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Figure 9.23: Wavelength vs. plasma flow speed for plane wave and turbulent cases 



fundamentally different character than those with the target region gas puffing. 
Discuss distinctions that apply to this data 

Concept of onset condition for turbulence (critical Reynolds number) 

concept of transient behavior during transition to turbulence 

concept of anisotropy 

Concept of energy cascades 

Distinction between 2-D and 3-D systems 

Issue of non-uniquenss of field line flow 

Instabilities? growth rates 

examine any trends in velocity scan (use data set that ranges over fullest extent (lo to hi)) 
plot probability of waves vs. renolds number? 
plot dominate wavelength vs. reynolds number? 

plot reduced 1-D Energy spectrum, derived from 2-D fourier transform? 

Discuss limitation on highest resolvable k, due to ion diffusion and motional blur. 

examine other diagnostics, LMP, B-probes, PMT's 
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Compare to hydro turbulence around a cylinder, (include image from NSF proposal) 
concluding remarks about what this data says about the critical re 
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Appendix A 



Magnetic helicity of a compact 



toroid 



The key idea in magnetic relaxation is that short timescale effects such as microturbulence 
and reconnection cause the magnetic fields to decay and become reconfigured in space in such a way 
as to find a state of minimum total energy. This minimum is non-zero because the total magnetic 
helicity of the plasma is a conserved quantity on the MHD timescales, and this provides a constraint 
on the system. The magnetic field can only become reconfigured in ways that preserve the value 
of the magnetic helicity A • B integrated over the volume of the plasma. In order to measure this 
constant 



we need to determine the vector potential A. The magnetic field is curl of the vector potential, and 
so A is determined only up to the gradient of a scalar function. 

For the case of a force-free field it is sufficiently general to formulate A in terms of the 
magnetic field and two scalar functions 7 = 7(x) and tp = ^>(x) according to 




(A.l) 



A = 7B + V ip 



(A.2) 
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The helicity is gauge invariant for a bounded plasma in a conducting vessel, and the helicity integral 
simplifies to 

f A-BdV = [ -fB 2 dV 
Jv Jv 

The V i> term was dealt with by applying a vector identity in combination with V - B = 0, followed 
by the use of the divergence theorem and the fact that B • n = 0 at the walls of the conducting 
vessel. We see that the grad tp term vanishes. 



\/<P-BdV= / V-(i/>B)dV= / ipB-dA = 0 

V JV JdV 

For a force free magnetic field configuration satisfying V x B = AB the multiplicative factor 
7(x) that determines A can be found by taking the curl of A 

B = VxA = 7VxB + (V7)xB 

so 

AB = I B -^xB 

7 7 

Then taking the dot product with B eliminates the cross product term and we see that 

XB 2 = -}-B 2 
7 (x) 

so 7(x) = 1/A which is a constant. So the vector potential is uniquely determined (up to a scalar 
gauge) by the magnetic field and the force free eigenvalue A. The vector potential for a force free 
magnetic field is simply 

A = \b (A.3) 
A 

The magnetic helicity is then directly proportional the total magnetic energy 

H = I / B 2 dV (A.4) 
A Jv 

This result implies that once a plasma has reached a force-free state, the total magnetic energy 
cannot be reduced any further because it is constrained by the conservation of H. From then on, 
the magnetic field and the helicity can decay only at a much slower rate due to the bulk resistivity 
of the plasma, [ref A.M. Dixon Astron. Astrophys. 225, 156-166 (1989)] 
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Appendix B 



Magnetic 



computations 



for probe 



calibration 



In this appendix we present the details of a variety of approaches that were explored 
to tackle the problem of magnetic measurements in a finite conducting well. (See Chapter 2 for 
application.) 

B.l Plane with circular hole 

The simplest is to regard the probe well as simply a hole in the vessel wall which we can 
reasonably approximate as flat. The magnetostatic field near a circular hole in a thin flat conducting 
plate has a tractable analytic solution expressible as an integral that can be numerically evaluated, 
[ref jackson]The magnetic field H has the form H = — V 4>m where the magnetostatic potential <\>m 
is given by 



-H 0 y + & 



(x) if z > 0 



4>M = < 



_$1 



if z < 0 
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This form guarantees that H = — V </>m will approach a constant vector of magnitude H 0 
in the y-direction, as you move away from the hole in the upper half space, while the field goes 
to zero in the lower half plane. Thus, describes the localized variation away from complete 

uniformity. The solution of V 2 (& 1 =0 given the boundary conditions of a hole of radius a can be 
found by Bessel series expansion in cylindrical coordinates (p, <fi, z) 

9H n 2 r°° 

$!(£) = — °— / dkj 1 {ka)e- k \ z \j l {kp)sm<f) (B.l) 
P l Jo 

where j\{x) = — j s the spherical Bessel function of order 1, and J\{x) is just an ordinary 
Bessel function of the first kind, order 1. 

This solution is attractive to work with because it is fully 3-dimensional, so the total flux 
through a coil oriented arbitrarily in space can be accurately calculated. However, the simplicity of 
the boundary geometry yields some quantitatively different effects than the actual, more complex, 
probe well geometry. The root cause of this discrepancy is that in the case of the simple hole, the 
field is allowed to expand out infinitely into to void beneath the plane, instead of just the finite 
sized well that we actually have. This means that more total flux enters the hole than would enter 
the opening of a well. Since more flux enters, the field beneath the plane is stronger for the hole 
geometry than the well (i.e., the field below the hole is diminished less than it would be inside the 
well). Then the calibration factor j^y- due to field diminishment caused by expansion, would be 
somewhat less in the hole geometry than for a finite well. 

Secondly, each field line that enters the hole will go deeper than it will in the case of a 
finite well. Thus for a given probe position, the resulting effective depth will be noticeably larger 
for the simple hole than it would be for the finite well. 

Since we are primarily seeking an accurate estimate of the effective depth of the probe in 
order to compare measured field components to those predicted by theory, we need to do a better 
job of including the effects of the finite well geometry. So, let us consider a few more possible ways 
to approach the full magnetostatic problem for a finite well. 
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Figure B.l: Analytic result for magnetic field near a hole in a conducting plane. Here the field 
vectors are represented by a color-wheel type coloring of the complex plane, which corresponds to 
the (y, z) directions in physical space. Red indicates a vector pointing directly in the -y direction, 
while other colors represent vectors at other directions as indicated by the color legend in top left. 
Vector magnitude is represented by a color's brightness, with black for the zero vector. This color 
code avoids the "pick-up-stix" effect that can occur with arrow-field representations. A sampling of 
field lines are also plotted to guide the eye. 

B.2 Analytic approach to full problem 

The most ambitious route would be to perform a similar Bessel series expansion of the 
potential as above. However, there is a significant complication that two distinct regions of space 
would need to be analyzed separately, and then a unique solution to the full problem could be found 
by matching series coefficients at the interface of the two regions. 

The natural choice would be to set one region as the cylindrical volume inside the well, 
with conducting boundary conditions on the side and base, and the interface matching condition 
on the top surface. Then the second region would be the upper half space above the plane, with 
conducting boundary at the z = 0 plane, everywhere except for the circular disk where the field needs 
to match the other region. To implement this, first some analytic simplifications to the general Bessel 
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series would need to be made that take into account the different symmetries of the two regions. 
Then you would try to see how matching the coefficient at the interface would constrain the form 
of the solution. However, at this point, progress toward a purely analytic solution begins to fall 
into that grey unknown between difficult and impossible, and a numerical solution of the system 
of coupled integral Fourier-Bessel equations becomes the only tractable option. Then considering 
that the primary aim of my work was not to spend weeks or months to perfectly solve an advanced 
3-D magnetostatic problem by numerical evaluation of Fourier-Bessel equations, I decided to find a 
simpler, quicker approach. 

B.3 Conformal mapping approximation 

Since the diameter of the probe coil is relatively small compared with the diameter of the 
port hole, a 2-D approximation can be made that should be quite accurate for the fields near the 
central plane of the probe port well. Field expansion effects should be more faithfully recreated by 
a 2-D solution of the right cross-section geometry than a 3-D solution of the wrong cross-section 
geometry. 

Essentially, the problem in 2-D is equivalent to a field above a conducting plate with an 
infinitely long trench cut into it, with a rectangular cross-section. If the trench runs in the x- 
direction, then the solution of the field will be independent of x, it will only vary in the 2-D (y,z) 
plane. This 2-D solution would correspond to the x = 0 plane in the full 3-D problem. 

One standard analytic method for dealing with this 2-D problem is the use of conformal 
mapping, [smythe] [F.M. Henderson, Elliptic functions with complex arguments, Univ. Michigan 
Press 1960, pg 10]. In general it is possible to transform the upper half of the complex plane into 
the interior of a polygon, by a conformal map according to a straight-forward formula. 

For the case at hand, the resulting polygon would be the "upside-down top hat" shape of 
the vacuum-conductor boundary. It would have four vertices, and interior angles of {^f , f , f , }■ 
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The original untransformed plane is the z-plane, and the transformed plane is the w-plane; the 
mapping function is w — w(z), which is formulated implicitly through a differential relation: 



dw (z — c) 2 {z + c) 2 
dz (z-i)*(z + i)% 



G 



(B.2) 




z - plane 




w - plane 



Figure B.2: Conformal mapping method for 2-D approximation 

Here, G is a complex constant that accounts for overall scale and orientation of the resulting 
polygon, and c is a positive real number, where the points z — ±c get mapped to the vertices at 
the opening of the well w — ±a + bi, which are labeled points D and A, respectively (see Fig 5). To 
solve for w(z) we need to take the indefinite path integral, where the integration dummy variable u 
takes a path from m = -oo to u = z, which is route- independent as long as no branch cuts (along 
Im z = 0) are crossed. The integral is 



— oo 



'!-(£) 



u(z)=Gl x j- r - Y du = Gcl \J^^fdu (B.3) 



This integral is not possible to evaluate "analytically" in terms of elementary functions, and 
does not appear in any standard table of integrals. This is because it corresponds to a combination 
of elliptic integrals of the second kind, usually notated by E(fc, z), here taking a complex argument. 

Some formal progress can be made from this realization, 
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w(z) = Gc[E(l/c,z)-E(l/c,oc)} (B.4) 

and the resulting dimensions of the well radius = a, depth = b, are determined from G and c 
according to 



a = Gc 



£(l/c)-E(l/c,oo)], b = Gc[k (^i - (i/c) 2 ) -E(yji- (l/c) 2 )] (B.5) 



where K (k) and E(k) are complete elliptic integrals of the first and second kind, respectively. The 
shape of the well (i.e., the ratio b/a) is then independent of G and has the form 

a E(l/c) -E(l/c,oo) 

It was not too difficult to implement a direct numerical evaluation of the path integral 
(5.3), and some qualitatively good initial results were obtained. 

Once you have the mathematical mechanics set up, it is very easy to find the shape of 
field lines in the vicinity of the well. You just sequentially evaluate w p = w(z p ) N times, where 
z p = (xq + p ■ Ax) + iy for point index taking values of p <E {0, 1,2,..., N}. This just maps the 
horizontal line Im z = y = constant to the corresponding field line. The resulting array of complex 
numbers w p can be plotted to visualize the field line in space. 

Furthermore, the effective depth can be found by tracking field lines away from the well 
and then finding how high above the wall they are after going several hole-radii away from the well. 
Field expansion ratios can be calculated by comparing the distance between two neighboring field 
lines, when they are far from the well, vs. when they are at the center of the well. 

Since the magnitude of the magnetic field is proportional to the density of flux lines, in this 
way we can get the ratio of the far field to its corresponding value of the distorted near field. To 
measure the average flux through the probe coil, we just choose one field line that passes through 
the upper tip of the coil, and a second that passes through the lower tip of the coil. If the distance 
between these two field lines far from the probe well is 5z and the height of the probe coil is Az co u 
then the calibration factor due to field expansion is C exp = Az s "° il ■ 
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The one stumbling block that prevented the full application of this method to solve the 
problem at hand is that, for reasons not yet understood, the ratio of b/a was not precisely controllable 
with the algorithms I have implemented. Essentially, only a single value of b/a ~ 1 was achievable, 
regardless of the choice of c and G. An intrinsic part of the problem is that my algorithm was 
divergent for points close to the real axis in the z-plane. This prevented me from making an accurate 
determination of what the real b/a ratio was for any particular instance of the mapping. There is 
some hope that after some further work with elliptic integrals, the increased insight and experience 
may lead to a resolution of the bugs in my code. But for the sake of time management, I abandoned 
this semi-analytic approach and switched to an even simpler fully numeric 2-D iterative relaxation 
method. 

B.4 Gauss-Seidel iterative relaxation 

We begin this final (successful) approach with a significant practical simplification. The 
magnetic field vectors themselves are not really important to find the effective depth and expansion 
factor, especially if we are only considering a 2-D case. Really we just need to know the shape of the 
flux surfaces on which the magnetic field lines reside. The simplifications arises from the fact that 
the magnetostatic potential <pM is what we call the harmonic dual of the electrostatic potential <j>E- 

The concept of harmonic duality is that if two scalar functions satisfy the Laplace equation 
(V 2 0 = 0, which is the definition of harmonic here) on the same boundary geometry, we say that 
the functions are dual to each other if one function satisfies Dirichlet boundary conditions = 
constant), while the second function satisfies Neumann boundary conditions = constant). The 
result is that equipotential surfaces of <f>i correspond to flux surfaces of 02, and vice versa. This 
is due in part to the fact that equipotential surfaces and the flux surfaces for the same potential 
function are locally orthogonal to each other. To go from a potential function to its harmonic dual, 
we just have to interchange the role of equipotentials and flux surfaces, and the resulting potential 
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will still satisfy the Laplace equation, but with the dual boundary condition. 

The magnetostatic potential satisfies the Neumann boundary condition at the surface of the 
conductor, while the electrostatic potential for the same conductor geometry satisfies the Dirichlet 
boundary condition. Therefore, we can find the flux surfaces of the magnetic field by actually 
computing the cquipotcntial surfaces of the dual electrostatic potential problem. This is an advantage 
because Dirchelet boundary conditions are intrinsically simpler in that you don't have to evaluate 
derivatives to apply them, and they lend themselves to numerical solution by iterative relaxation 
methods. 

I have implemented a simple Gauss-Seidel method using a 2-step iteration on two inter- 
leaved square meshes, [ref jackson 1.13 pg 47-50] 
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Figure B.3: Two interleaved meshes, one for each timestep, used for calculating magnetic field near 
probe port 

The iteration scheme divides the mesh into two lists of coordinates, corresponding to pair of 
interleaved meshes depicted in fig B.3. First, using list A of coordinates La = {{xq, yo), (xi,y%), . . .}, 
which correspond to the blue squares in fig B.3, we evaluate a new value for <j> that is equal to the 
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average of its neighboring values, including newly modified values on the same list. Then we repeat 
using list B (corresponding to the green squares in fig B.3). 

4 1 

Step A: <f>new(x,y) = ~((f>oid b) c + -^(finew a) s V(x,y)eL A 

4 1 

Step B : 4>new{x,y) = ~{4>old A) c + ~{4>new b) s V(x, y) € Lb 

where ((f)) c is the "cross" average of neighboring values in the opposite list, and (</>) s is the "square" 
average of neighboring values in the same list that is being currently updated, 

{4>{x,y)) c = i[<p(x + 5,y) + (p(x,y + 5) + <p(x-5,y) + <p(x,y-8)} 

(<t>(x,y)) a = | [<f>(x + S, y + 6) + <f>(x + S, y - 6) + 4>{x - S, y + 6) + <f>(x - S, y - 6)] 

where 6 is the pixel size of the mesh. 

The initial configuration at timestep 0 was such that <pM = 0 inside the well and along the 
walls, while a constant upward gradient in </>m existed in the upper half plane. As the iterations 
progress the field steadily expands into the well and relaxes into a solution of Laplace's equation. 
The program halts after a preset number of iterations (stop when i — N). We tried different values 
of N and recorded how the Laplacian error asymptotically leveled off toward a stable minimum. For 
a mesh of 200 x 200 pixels, after about N <~ 1500 iterations the largest error in the system was less 
than 10~ 2 (sec fig B.5). 

A quantitative analysis of the results are discussed in chaper 2. 
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Figure B.4: 2-D Relaxation simulation result of magnetic field near a conducting port well. Image 
on left uses a cyclic greyscale to visualize flux surfaces. Image on right represents the magnetic 
field vectors using a color-wheel type coloring of the complex plane, which corresponds to the (y, z) 
directions in physical space. Red indicates a vector pointing directly in the -y direction, while other 
colors represent vectors at other directions as indicated by the color legend in top left. Vector 
magnitude is represented by a color's brightness, with black for the zero vector. 
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Figure B.5: The log (base 10) of the value of the Laplacian of the magnetostatic potential near a 
conducting port well. Since the goal is V 2 0 = 0, this graph indicates the order of magnitude of the 
error in the calculation. 
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Figure B.6: Mechanical drawing of probe port at z = 142 cm showing measured dimensions. 
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Appendix C 

Method of reconstructing 
a moving pulse 

In the following discussion we present a general purpose method of analyzing a set of N time 
signals Fi(t), i S {1,2,..., N} that consist of local measurements of the same physical observable 
taken at the TV sequential diagnostic port locations of a linear accelerator device, with i being the 
port index (e.g. 57, 91, 142 for CTIX). We describe an algorithm that can reconstruct an accurate 
approximation of the true space and time dependence of the given observable F(z, t) in the regions 
between the points of measurement where diagnostic access is possible. This method has many 
conceptual parallels with the method of characteristics used in analyzing nonlinear wave motion, 
however we are looking at the problem from a vantage point that deals more with experimentally 
relevant issues of applied mathematics, and not at all with the issues involved in a theoretical analysis 
of some set of differential equations. We also explore the conditions needed for this approximation 
to hold, and we test the algorithm under a variety of circumstances. 
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C.l Sequential measurements of a moving pulse 

The motivation for developing this type of interpolation originates with the daily task of 
trying to interpret time-domain data from a set of spatially distributed diagnostic signals, on the 
fly as the data comes in during an experimental run. After some experience with the raw signals, 
a sort of mental-interchange of space and time coordinates is possible to visualize; resulting in an 
intuitive appreciation of the accelerator dynamics and an approximate understanding of how a given 
observable varies in space. Of course, it is often difficult to fully communicate this understanding to 
others who are less familiar with the nuances of the particular system under consideration. In this 
appendix we present the details of an algorithm that produces a quantitative result that makes the 
accelerator dynamics easy to interpret. 

Some sample data is presented below in fig. C.l. It shows the time dependence of the 
magnetic field measured by probes on the CTIX accelerator section. 
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Figure C.l: Example of time-signals from axially distributed probes. These are the three measure- 
ments of B z (t) taken within the CTIX accelerator at positions z = 57, 91, 142 cm. 
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This work has been primarily motivated by signals like this and the desire to use them to 
determine the true magnetic geometry of the compact toroid plasma. Essentially we would like to 
know what is the function B(z,t) that describes the magnetic field as a function of axial position and 
time (see chapter 2). Additionally a variety of other experimental observables are also measured 
at sequential positions and can benefit from this interpolation analysis. For example, these signals 
could be the average electron density visible light emission, or plasma probe current at several axial 
positions as a function of time. In chapter 6 we used this reconstruction method on a sequence of 
visible emission measurements to determine plasma luminosity as a function of space and time. 

It is interesting to note that a similar algorithm is used in video decompression schemes 
when trying to interpolate in time between actual frames of a video image, in order to reconstruct 
a better video at a higher frame rate.[ref,ref] The motion of distinct features within an image can 
be taken into account to predict where each feature will be at any intermediate time t, according to 
a quadratic model of its trajectory. In this problem the observables are vectors, and the details of 
any practical solution ultimately rely on some advanced image analysis. 

While the problem of one dimensional flow considered in our work is in some sense a lower- 
dimensional special case of the motion-compensated interpolation used in video analysis (where the 
analogy requires an interchange of the space and time coordinates used in the two problems) it is 
worth formulating this application on its own in a coherent and more directly usable way. 

The algorithm we have developed can be broken down into two distinct sub-processes. The 
goal of the first process is to estimate the fluid velocity field for all (z, t) in the system, by any means 
necessary. The method used to make this estimation will need to rely on specific physical properties 
of the dynamical system under consideration. What are the basic constraints on the velocity field? 
Can "fluid/ wave" elements pass through each other with several different velocities at the same 
point, such as would happen with a superposition of traveling waves in a dispersive medium? If this 
is the case, then an accurate reconstruction of the field observable and its wave phenomena would 
need to begin with a good characterization of all of different wave velocities present. 
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This kind of complication is far beyond what has been considered so far; we will only deal 
with the much simpler case of fluid elements that never interpenetrate along the axis of motion (i.e., a 
single valued and continuous velocity field). This seems to be sufficient to describe the experimental 
observations of magnetohydrodynamics on the CTIX accelerator, moreover, the assumption of single 
valued velocity is consistent with basic intuitive reasoning about the evolution of a magnetic field 
that is frozen into a moving fluid. 

The second process takes this estimated velocity field and uses it to create fluid element 
trajectories, and finally computes a unitary-weighted superposition of the adjacent real data to 
determine the interpolated value of the observable at intermediate positions along each trajectory. 

C.2 Estimation of flow field 

On the most fundamental level, everything will depend on some incarnation of a mapping 
function that defines a correspondence between points on one probe signal with points on the signals 
of neighboring probes. This kind of mapping function is essentially a way to label the fluid elements 
in some unique way so that they all can be tracked. This tracking information enters our awareness 
in terms of measurements of the time at which each fluid element passes by each probe location. 

The criteria that defines which pairs of points correspond to each other is really a very 
system-specific issue. To do a good job with this, some consideration of the physical dynamics will 
be needed. And each distinct physical system may need its own bag of tricks. Overall, there are two 
basic varieties of possible point-matching criteria, which can be classified as being either point-wise 
in nature, or having a dependence on the collective properties of some subset of the data points in 
the signals. 

We will now list a number of important examples of point-wise matching criteria. 

1) Peak-to-peak This is a discrete valued map that only tracks a single fluid element through 
the system. This criterion defines the fluid element of peak amplitude, and we measure the time it 
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crosses each successive probe position by finding the time at which each probe signal reached its 
maximum value. We can generalize this to track the fluid elements that ride along with any of the 
local extrema in the field. 

2) Level set labeling This is a continuous method that might be useful in certain non-dissipative 
systems, where we identify two points on adjacent probe signals if they have the same amplitude, 
and exist within corresponding "regions", which is an additional constraint needed to prevent the 
labeling from being many-to-one. It is often reasonable to partition the domain into regions that 
fall between successive local extrema in the signal. In this way we ensure that every point is labeled 
in a unique way. This type of labeling can be extended to dissipative systems if we can normalize 
all of the probe signal to give them the same range of values, yet still preserve the correct kinematic 
behavior of the system. 

3) Differential labeling The simplest instance of this type would be to make a correspondence 
between points of maximum (±) slope in the probe signals. Clearly, this could be extended in a 
variety of ways similar to those already mentioned. 

The second category of matching criteria, collective correspondence maps, define relations 
between groups of points in one probe signal to corresponding groups of points in the neighboring 
probe signals, often with no explicit relation between pairs of individual points, instead it is implicitly 
determined via the collective relation. There are two important examples relevant to our work. 

1) Correlation methods Ideally we want to find the time-of- flight of each distinct fluid element 
as it travels between adjacent probes. We can begin to do this by calculating the correlation function 
between adjacent probe signals. 



The global maximum of Cj(At) occurs when At results in the optimum translational correspondence 
between the two signals. 
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For a moving pulse that does not change shape very rapidly as it travels, this method 
has certain advantages over point-wise maps, such as the peak-to-peak method, which can produce 
erroneous results due to any moderately high frequency noise occurring near the peak of the pulse. 
Instead, the correlation averages out any noise and finds the best overall fit between the signals. 

For any real system in which acceleration is occurring, this is not enough to properly esti- 
mate the velocity field. One conceptually attractive approach is to make a first order approximation 
using the global correlation to match the peaks of the signals, then go back and make a second 
order correction, again using correlation integrals to calculate the optimal translation, but this time 
locally, over a sequence of subwindows that progress away from the peak toward (t = ±00). The 
result is a time-shift function Atij(t) that defines the correspondence between points in adjacent 
probe signals Fj, and Fj, according to Fi(t + Atij(t)) ~ Fj(t). The time-shift function depends on 
the first and second order correlations in a straight-forward way, Atjj(i) = AiW + At^(t). Overall, 
this method would find the optimal (non-uniform) stretch to make one function fit another. 

However, there are a number of difficult pitfalls involved in actually calculating this second- 
order correction. One problem is that any overall slope in the signals across the length of the 
subwindow will lead to the maximum value on a boundary point always being optimally correlated 
with the largest neighboring peak in the second signal, regardless of any real coincidence between 
the ups and downs of the two signals. To make this idea work requires some sophistication, possibly 
fitting smooth curves to the signals in the subwindows, and then only correlating the residuals from 
this fit. An added complication is the question of how to iteratively evaluate and then apply the 
local time-shift function, as the subwindow moves away from the peak. 

A more important issue with this method is that not all moving pulses will leave time- 
signals that can be used to reconstruct the original pulse via correlation methods. The primary 
pathological examples are dispersive pulses consisting of several distinct wave components, each 
traveling at a different phase velocity. Even the case of two wave components can lead to serious 
errors in estimation of wavelength of the individual components, and the overall shape of the pulse. 



167 



This error is simply due to the correlation method incorrectly matching wave crests. Basically, 
this method uses the overall group velocity of the pulse, when it should be finding distinct phase 
velocities. 
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Figure C.2: Example of waveform that causes problems with the correlation method (orange curve). 
A simple remedy exists (green curve), but this example shows the limitation of the correlation 
method. 
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The error in this reconstruction occurred because we were simply using the wrong velocity 
field. If we carefully match wave crests correctly by hand, and then simply use piecewise linear 
extrapolations to estimate the velocity, we do significantly better in our reconstruction. 

This goes back to my point about the importance of including your best understanding of 
a system's dynamics to make the velocity estimation. Correlation methods should not be used when 
they are in disagreement with the basic physics of the system. However, when multi-wave effects 
are unimportant, successfully implemented correlation methods can be robust and elegant. 

2) Kinematic modeling When some kinematic property of a dynamical system is well estab- 
lished, it can be used to directly estimate the velocity field. The correspondence between points in 
neighboring probe signals will be entirely implicit in this method, but could be found by calculating 
trajectories from the estimated velocity field. 

Implementation specific to CTIX Consideration of the acceleration dynamics leads to a rela- 
tively simple model of the velocity field 

C.3 Lagrangian interpolation algorithm 

This method can be simply described by the following general formula for the interpolated 
value of the observable F(z, t) at position z between the ith and (i + \)th real probes, which involves 
a unitary superposition of neighboring probe signals, evaluated at retarded and advanced times: 

F(z, t) = Wi (z, t)Fi{t + AU(z, t)) + (1 - Wi(z, t))F (t+1) (t + At i+1 {z, t)) (C.l) 

where Wi is a weighting factor (< 1) dependent on the fractional distance between the position z of 
and the positions of the real probes. Although Wi can in general depend on z and t, the simplest 
good implementation is actually independent of time, Wi(z) = (z — Zi)/(Z i+ i — Zj) where Zi is the 
position of the ith real probe, i G {1, 2, ...(N — 1)}. The function Fi(t) is the time signal of the 
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ith real probe, and Ati(z,t), Ati + \{z,t) are some delay functions that will need to be optimized 
for best agreement with the approximately measured flow dynamics inferred from the time of flight 
between sequential probes. For fluid that is moving in the positive z direction, a given fluid element 
first passes by the lower probe at Zi, then it passes through the value of z under consideration for 
interpolation, and lastly it passes the upper probe at Z i+ i. And so with the sign convention used 
in (C.l) we see that Att(z,t) will have negative values, while At i+ i(z,t) will have positive values. 

Examples of the application of this method to actual data from CTIX have been presented 
in chapters 2, 4, and 6. 

C.4 Trajectory calculations 

[figure, illustrate how flow field is used to find the timeshifts for each point on the grid] 

C.5 Validity of approximation 

First the good news. If carefully implemented, this interpolation is 100% accurate at all of 
the locations where the real measurements have been made. Of course, this must be so for any kind 
of interpolation, and formula (C.l) is no different. 

In the following subsections we examine two fundamental categories of errors in this al- 
gorithm. The first kind of error can occur in the estimate of the velocity field. This is the most 
important practical concern because if the velocity field is wrong then the final reconstruction will 
probably suffer from a variety of serious errors. The second kind of error can occur even if you have 
a perfect formulation of the velocity field, and we will explore the hypothetical worst case that is 
mathematical allowable, in which transient pulses appear and disappear in the physically real field 
variable, but do so in the regions between probe measurement and thereby go undetected. 
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C.5.1 First order velocity error 

Here we will examine in some detail the errors produced in a certain formulation of a point- 
wise velocity estimation method that has been particularly useful on CTIX data. We need to find 
the error in the velocity curve of the CT/pushing- field interface (CT back edge), since any error at 
this stage will have an effect on the entire reconstruction. The key idea is that there will be some 
time error At in the measurement of the time that the CT back passes each probe position, and this 
error is dependent on the random fluctuation of the signal level, and inversely proportional to the 
slope of the signal at the crossing time. At = -^AB Z . 

[figure, illustrate how error in Bz causes Delta t] 

First we need to estimate the signal fluctuation level. Here we will investigate the high 
frequency noise in the B z signal as being the source of any errors in computing the velocity of the 
CT as a function of time. We will define the signal fluctuation to be the maximum absolute value 
of the residual of the B z signal after subtracting off a smoothed version Sj(i) of the original signal. 



AB Z = max 



B z (t) - lim Si(t) 



However, the result of this is dependent on exactly how much smoothing is done. The parameters 
of the limit, & and £ s top, simply represent a stopping condition for the smoothing process, and will 
be defined shortly. To provide an objective criterion for sufficient smoothing, and thereby yield 
uniform results in this analysis, we will adopt the following procedure to implement the smoothing. 
The primary concern of this procedure is that it smooths enough, but not too much. The fluctuations 
that we identify as noise in our B z signals typically exhibit a rapid change from positive to negative 
slope over a short interval of time. And so the "true" signals (as we conceptualize them) that lie 
underneath the noise will not have this "jumpy" character to them. As a reasonable measure of this 
"jumpiness" factor we can consider the first time derivative of the partially smoothed signal Si, and 
then multiply it by a weighting factor determined by the duration that the signal is near that value 
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of time-slope. The result is a function of time 

jumptt) = 3—— r • ^li^j 

Duration(slope) at 

Here, the term Duration(slope) is simply the time between zero crossings of the time 
derivative, and it measures (in a piecewise fashion) the duration that the signal has either positive 
or negative slope. The reason for dividing by the duration of same-sign slope is that the physically 
real signal may have large slopes, but it will typically maintain this slope for some extended period 
of time; whereas noise will change the sign of its slope every few digitization steps. 

As the original signal is iteratively smoothed, the noise will be removed, and the maximum 
absolute value of jump(t) of the ith partially smoothed signal Si , will fall monotonically. Let us 
call this jumpiness parameter £ = max \jump(t)\. The smoothing process should stop when £ falls 
below some small predetermined value £ s t op , which represents the goal for smoothness of the noise- 
free signal. It can also be thought of as corresponding to the fraction of total noise that has been 
removed from the signal, (1 — £ s *op) ~ (removed noise) / (total noise). The reason for following this 
procedure, is that we can make a single subjective choice for a reasonable value of £, s to P , and then 
use the same value across the board when analyzing all of the magnetic data; yet each signal will 
yield a different value of fluctuation level AB Z depending on its noise properties. This is far more 
objective than estimating AB Z by eye every time you do this analysis. 

For example, here is the Bz 57, 91, and 142 signals for shot ?. With a value for the smoothing 
limit of ^ s top = 0.01 we will tabulate the estimated signal error AB Z , and a corresponding arrival- 
time error ±Ai for each probe. Lastly we list the relative error of the final velocity of the CT, 
derived from the maximum possible error considering all possible sign combinations for each of the 
arrival time errors. 

Note that if you had a set of three new arrival times that were all shifted in the same 
direction (cither + or -) by the same amount At, then the calculated final velocity with these 
"errors" would be exactly the same as the true final velocity. Only if there is some asymmetric 
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distribution of positive or negative errors, will the resulting error of the final velocity reach its 
maximum value. Here we list this worst case velocity error as an upper bound, 
[table] 

If we are concerned with the qualitative effects of this level of error, we will pleased to see 
that it is barely noticeable. 

[figure: graphs comparing two reconstructions with an 8% disagreement in final velocity] 

C.5.2 Constraints on the amplitude of undetectable transient pulses 

In order to estimate the reliability of this reconstruction method we need to carefully 
consider the worst possible errors that could exist in the space between consecutive probes. We will 
consider errors that are as large as the mathematics of the algorithm allow, with no concern initially 
for the physical limitations on such errors. The end result of this analysis is a constraint in the 
form of an uncertainty relation between the duration and spatial extent of any undetected transient 
pulse-like modulations of the field occurring in the unmeasured region between the probes. This will 
provide a solid upper limit to the total error of the reconstruction. 

Since the probes themselves yield extremely accurate and reproducible measurements of 
the field quantities in the immediate vicinity of the probe locations, there is very little error due 
to global fluctuations of the real field, since these would be detected simultaneously by multiple 
probes. Instead, the real cause for concern are transient fluctuations, that cause error because they 
are short enough in duration and spatial extent, and happen to occur deep enough into the empty 
space between probes so that they go undetected, and consequently the algorithm has no ability to 
include their existence in the reconstructed waveform. 

It is tractable and informative to consider transient pulses that have a Gaussian form in 
space and time (z,t) such as 




(C.2) 



173 




Figure C.3: Transient pulse at various times and its time-envelope as a function of axial position. 



here, A is the amplitude of the pulse, r is the characteristic rise/decay time of the pulse (i.e., its 
duration), L is the characteristic spatial extent, t p is the time at which the pulse reaches its peak 
value, and v is the pulse group velocity, which we will assume is approximately the average fluid 
velocity. To measure the damage this pulse could do to our reconstruction, we need the envelope of 
the pulse as a function of position z, which is the maximum value of the pulse at a fixed z, over all 
values of time. 

First we find, for each z, the value of t maximizes J-(z,t). For positive- valued Gaussians 
like we are considering, there is a finite maximum value but no minimum value for finite (z, t). The 
maximum occurs when dT jdt = 0, 
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And so the maximum occurs when 

t — t p v(z — vt) 

~ T 2 



= 0 



2 2 \ 2 2 

V T \ V T Z 



or t 1 + —pr = t P + 



L 2 J p L 2 

Now, to clarify the algebra, let us define a dimensionless parameter n — vt/L, and we see 
that for a fixed z the pulse reaches a maximum value at time 

2 

t 4- ^-2 

- WW - (1 + K 2) 

The envelope is given by T env (z) = J 7 (z 7 t max (zj), which simplifies to 

/ z — vtp \ 2 

Restated in terms of of r, (using Lyjl + n 2 ) = yj L 2 + v 2 t 2 )), this is 



r env {z) = Ae W^+^)J (C3) 

In order for this pulse to be undetected by neighboring probes located at za and wc 
would need the pulse envelope to be less than some noise floor c no i se 

°nv(.ZA) ^ ^noise and 3~ env{zB^} ^ ^noise (^'^) 

where e no ise is determined by the statistical fluctuations that are present in an ensemble of measured 
signals under equivalent experimental circumstances. It does not have to be restricted to just 
the electrical noise present in the probe diagnostics, but could include actual MHD modes in the 
plasma that constitute random deviations away from the average waveform of the probe signals. 
It is at least conceptually possible that occasionally these real fluctuations would occur at a larger 
than normal amplitude and in a transient fashion between the probes, and thus go undetected. In 
the worst case scenario the peak of the transient pulse would occur at exactly the half-way point 
between consecutive probes, thereby allowing the largest possible amplitude allowable by (C.4). 
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From this (unlikely) situation, we could make the most conservative estimate for the accuracy of 
this reconstruction method. The true accuracy should only be better than what we find here. The 



worst case is when vt z 



nondctection is 



\{za + zb), and so let AL pl . 0 f, e = 1{zb — Zmid)- Then condition for 



A 



And the maximum possible undetectable amplitude is constrained by 

A < £ n c 



< e r , 



AL pro b e 



(C.5) 



(C.6) 
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Figure C.4: Log (base 10) of the maximum possible amplitude of undetected transient pulses, plotted 
as a function of pulse duration and spatial extent. 



This maximum possible amplitude is plotted in fig (C.4) as a function of fluctuation dura- 
tion and spatial extent (r, L) using parameters typical for the CTIX experiment. Using this result 
we can make statements regarding an upper bound of the error of the reconstruction method, such 
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as: 



For an average fluid velocity of v = 20 cm/ [is we can say that in the region between 
za = 91 cm and zb = 142 cm, if there exists an undetected transient pulse that spans 
the space between the probes, and it exists for at least the transit time between probes 
[2L = 51 cm,4 L = 20.5 cm and 2t = (51/20)/zs, r = 1.275/xs], (sec fig C.3) then the 
amplitude of the pulse can be no larger than 

Am.ax 1.84493 • 6 nr)7 " sf > 



For pulses that are significantly shorter in duration and extent, the upper limit on amplitude 
becomes much larger than this example. In fact, formula (C.6) goes to oo in the limit of zero pulse 
length and duration. Ultimately however, this does not cause any serious problems if our goal is to 
reconstruct the average waveform of the field moving past a set of probes, where we care about a 
spatial resolution of several cm, and temporal resolution of a few [is. 

Also, there is a mitigating factor involving the probability of these hypothetical high- 
amplitude events. If we take the noise floor to be the expectation value of the natural fluctuations 
T that occur with the plasma, e noise = (J 7 ), then we can use Markov's inequality [ref] to bound the 
probability that a high-amplitude event will occur during any given shot. Let J 7 be a nonnegative 
random variable; then for any amplitude of fluctuation A > 0, the probability that T will exceed A 
is bounded according to 

P(T>A) < & 

For example, the probability that a fluctuation will exceed 100 times e noise , is at the very 
most 1/100. And so, although high amplitude errors are possible, they are constrained by basic 
probability laws to be correspondingly rare occurrences. Alternatively, if there is some physical 
reason for some finite upper bound on A, such as conservation of energy perhaps, then we can 
restate equation (C.6) in the form of an uncertainty relation 



A T 2 

L 2 + v 2 t 2 < probe . (C.7) 

We see that if A < e noise then L and r are unconstrained, which agrees with the fact 
that small amplitude fluctuations have already been included in our definition of e no ise- But as A 
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becomes larger, the region of spacetime with a significant reconstruction error must become smaller 
and smaller, and its occurrence becomes increasingly improbable. Lastly, for some fixed A > e noise} 
we see that there is a trade-off between the spatial extent and the duration of undetectable pulses, 
and that this relationship depends on the fluid velocity and the distance between consecutive probes. 
The best reconstruction occurs when consecutive probes are closely spaced, and there is a high flow 
velocity. 
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Appendix D 

Equations of incompressible 
resistive MHD 

The equations of magnetohydrodynamics (MHD) describe the evolution of the velocity 
field, v(x, t) and magnetic field, B(x, t) within a conducting fluid, and also how they couple to each 
other. We call this description ideal MHD if we consider case of infinite electrical conductivity, 
while it is called resistive MHD when we take into the effects of more practical, finite values of the 
conductivity. 

There exist certain symmetries between the two fields that is worth notice. Both v and B 
are solenoidal 

V-v = 0, V-B = 0 (D.l) 
Also, the definition of fluid vorticity closely parallels that of current density via Ampere's Law 

u = Vxv, j=VxB (D.2) 

The vorticity evolves in time according to a Navier-Stokes equation that includes the effect 
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of the Lorentz force term (j x B). The dissipative term can be generalized in a variety of ways. 

3,w-Vx(v xw+jxB) = (-1) v - 1 ij, v A v u> (D.3) 

Here \i v is the fluid viscosity and the exponent v determines the type of diffusion operator that is 
used, v = 1 gives the usual Navier-Stokes diffusion according to the Laplacian of w, while v = 2 
gives "hyperdiffusion" . 

The magnetic field also evolves in time according to a convection-diffusion equation: 

a ( B-Vx(vxB) = (-ir 1 ^A"B (D.4) 

Mo 

Where r\ v is the electrical resistivity of the fluid. Be aware that equations (D.3) and (D.4) 
are nonlinearly coupled through the convection terms. 

One subtle point, although the above equations give an accurate and deterministic descrip- 
tion of the evolution of the vector fields v(x, t) and B(x, t), there is not an absolute meaning to the 
evolution of any particular field line[ref, ref]. At any particular instant in time (ti), we can correctly 
trace out a field line that is locally parallel to the magnetic vector field at every point along its 
length. However, we can not uniquely identify the future version that field line at some later time 
(t 2 ). It could correspond to any possible field line that exists at time t 2 . This is true even if the 
fluid velocity field is known, and the resistivity is low enough to consider that the field is "frozen" 
into the fluid. The mathematical argument for this non-uniqueness goes a bit like this: 

Given that any magnetic field can be expressed in terms of two scalar functions a(x, t), 
and (3(x, t) by 

B(x,t) = V a x V/3 

then any field line at a fixed time is uniquely delineated by a pair of coordinates a — constant, 
0 = constant. In other words, any magnetic field line is the curve of intersection between some 
surface of constant a, and some surface of constant j3. If we apply this formulation for B to its 
evolution equation (D.4), and consider the perfectly conducting case (n — 0) to ensure that the field 
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is frozen-into the fluid, 

d t B- Vx (v x B) = 0 
we find that the fluid velocity cancels out and we have, 



Vx 



da . dB 



= o 



Due to the curl operation, this relation for the evolution of a(x, t), and (3{x,t) is only 
determined up to the gradient of a some arbitrary scalar potential <p(x,t), 

da _ dp 
at dt 

Thus, depending on your choice for <fi, an initial set of scalar fields a(x,ti), (3{x,t\) will 
evolve into any other possible scalar fields a(x,t2), (3(x,t2), with the intersections between level 
surfaces also changing in an arbitrary manner. In the case of (f> = 0, we find that the field lines do 
indeed flow with the fluid elements. But for any other <j>{x, t) with a non-zero gradient, this will not 
be the case and the field lines will have a relative motion to the fluid elements. 

Overall, the moral of this story is that the equations of MHD describe the vector fields 
v(a?, t) and B(af, t), but they say nothing about the evolution of magnetic field lines. And we must 
adjust our intuitions accordingly. 
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Appendix E 

Vector field visualization 
using a complex color code 

purpose of technique, description of method, black at center, white at infty, color deter- 
mined by vector direction 

convert vector to complex number 

choose "radius" of magnitude to fit brightness range to full extent of vector size range, 
formulas for rgb values 

vector field example showing problem with arrow representation, and advantage of color 

code 

formula for vector field [figure, pick up sticks] [figure, same vector field in color code] 



-4 -2 0 2 4 
Real 



Figure E.l: Complex color code 
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Appendix F 

Direction of future research 

1. Make improved measurement of the effective depth of accelerator region magnetic probes. 
Possibly make new probes whose position and orientation are controllable in a more accurate way. 
New probes could also include a third coil to measure B r . Effective depth could be calculated from 
actual measurements of a known sinusoidally varying field taken over a range of probe depths 

2. A major avenue of research that would advance our understanding in a variety of ways, 
would be to perfect Dopplcr measurements of ion velocity, ft would be of great theoretic value 
to known the velocity field down to a cm, sub — jj,s scale. If it was possible to simultaneously 
measure this at several locations fluctuations in velocity could be examined, in principle. A new 
high resolution VUV spectrometer is planned to be installed on CTIX in 2006. 

3. Conduct improved measurements of CT magnetic geometry, in a variety of situations, 
improve estimates of resisitvity. develop a real time resistivity computation that can be displayed 
during an experimental run. measure radial dependence of fields to compare with theory. Work out 
equilibrium with inclusion of Hall term invariants. 

4. Comparison with 3-D simulation. 

5. turbulence work 

6. improved versions of CTIX 



